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Analysis and Design of an Oscilloscope Deflection System With 
a Calculable Transfer Function 


D. M. Stonebraker* 


Institute for Basic Standards, National Bureau of Standards, Boulder, Colo. 80302 


(December 12, 1967) 


An oscilloscope deflector is described, which has a calculable transfer function. The deflector 
is analyzed, leading to its transfer function, in terms of the complex frequency variables. A practical 
strip line deflector, usable as a pulse standard, is designed; and its frequency response, sensitivity, 
bandwidth, rise time, step function, and impulse response are calculated. The predicted deflection 
is down to 70 percent of its d-c value at 2.82 GHz, while the step-response 10 to 90 percent rise time 
is 148 picoseconds. The effects of a drift space on the oscilloscopic display is also discussed. Results 
are compared with a well-known expression for the parallel plate deflector structure. 


Key Words: Deflector, drift space, oscilloscope, strip-line, standard, transfer function. 


1. Introduction 


The cathode-ray oscilloscope is the most common 
and most useful instrument for displaying voltage as 
a function of time. It directly provides a graphical 
presentation of a time function, y(t), which is related 
to the voltage applied to the input terminals of the 
oscilloscope by the oscilloscope transfer function, 
H(@), providing we assume the oscilloscope is a linear 
device. This concept is depicted by figure 1. 


OSCILLOSCOPE 





OUTPUT 


v(t) INPUT VOLTAGE DISPLACEMENT y(t) 


je(w) 





H(w) = |HWw)| e 














FIGURE 1. Linear system representation of oscilloscope. 
When using an oscilloscope to measure a sinu- 
soidal input voltage, we generally assume or require 
that the oscilloscope be calibrated so that the “true” 
amplitude of the sinusoid may be read directly from 
the display. This calibration can be made by measur- 
ing the oscilloscope amplitude response, |H(w)|, 
from d-c to a suitable high frequency and either 
adjusting the instrument to read correctly or provid- 
ing appropriate calibration data to correct the dis- 
play amplitude. If on the other hand, the input voltage 
is a pulse or a pulse train, with high harmonic con- 
tent, the oscilloscope must be calibrated to cor- 


*Present Address: 3361 Vista Drive, Boulder, Colo. 80302. 


‘If a standard pulse is available, it may also be used to calibrate the oscilloscope for 
pulse displays. 


rect for both amplitude and phase distortion.' This 
requires measuring both |H(w)| and g(@), i.e., 
the amplitude and phase response of the oscilloscope, 
from d-c to many gigahertz. This is a formidable 
task owing to the lack of voltage standards above 
4 GHz and the lack of a technique for measuring phase 
distortion. A suitable alternative is to calculate 
H(w) from knowledge of the physical characteris- 
tics of the oscilloscope and verify the calculations, 
insofar as possible, by sinusoidal measurements. 
Once H(qm) is found, an input pulse may be calcu- 
lated from a knowledge of the output deflection, 
y(t), by either Fourier series or integral methods. 


2. Design Considerations 


A system analysis of a typical oscilloscope, either 
real time or sampling, is complicated by the nature 
of the input probe and amplifiers and by the de- 
flector used. Also the requirement of linearity must 
be satisfied for a linear analysis. With these criteria 
in mind, it was decided to select for analysis and 
design an oscilloscope with a direct feed (no ampli- 
fiers), feed through (no probe), 50-0, real time (non- 
sampling) deflector system. It was further decided 
to avoid slow wave (traveling wave) structure for 
reasons to be discussed next. 

Owaki [1]? designed and analyzed a traveling wave 
cathode-ray tube and obtained a mathematical ex- 
pression relating spot displacement to the input 
voltage and geometry of this tube. Hollmann [2] 
analyzed an elementary parallel plate deflector sys- 
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tem and obtained the familiar deflector sensitivity 
relationship for sinusoidal inputs of the form? 


sin w7/2 
w7/2 (1) 


relative dynamic sensitivity = 


where w is the angular frequency of the input voltage 
and 7 is the time taken for the electron to pass be- 
tween the plates. Talbot [3] considered the same 
problem and applied convolution integral methods 
in order to find various pulse response character- 
istics of the elementary deflector. He applied the re- 
sults to a direct feed, feed through, slow wave helical 
deflector. The results of these three workers analyses 
could be used to write approximate scope transfer 
functions for their oscilloscopes. However, in the 
case of Owaki’s and Talbot’s deflectors, some un- 
certainties in the transfer function would be difficult 
to estimate. These uncertainties stem from the 
effects of coupling and field pertubations between 
loops of the slow wave structures and resulting 
specification of transit time. Because we want to 
avoid significant uncertainties of this type, slow 
wave structures have been eliminated from consider- 
ation. The elementary parallel plate structure could 
be used, but its sensitivity is extremely small and it 
is difficult to terminate in its characteristic imped- 
ance. Lee [4] showed that this structure can, how- 
ever, be used in a microoscillograph with sufficient 
sensitivity for pulse work. Nahman [5] gives an ex- 
cellent survey of pulse oscilloscopes with a very 
complete list of references. 

Keeping all the criteria mentioned above in mind, 
as well as Lee’s microtechniques, the vacuum- 
filled strip-line deflector shown in figure 2 was selected. 


7 ELECTRON BEAM 






































FIGURE 2. Strip-line deflector. 


The electron beam travels between the two con- 
ductors in the region where fields, as shown by 
Brooke [6], are almost uniform if the center con- 
ductor is made sufficiently wide. 

In order to make a reasonably simple analysis of 
the relationship between a voltage applied between 
the two conductors and the resultant deflection of 
the electron beam, the following design criteria and/or 
assumptions were made: 

a. The deflector characteristic impedance is 500 


?Figures in brackets indicate the literature references on page 

3 Hollmann’s result does not include the smal] displacement within the deflector space. 
It considers only the transverse velocity change within the deflector and the effects of the 
drift space. Also Hollmann assumes the electromagnetic field between the plates is spa- 
cially invariant and ignores the traveling wave nature of the voltage applied to the plates 
and the attendant reflections at the open ends of the plates. 


and can be terminated in 500 so that any standing 
waves are negligible. 

b. The electric fields in the region where the elec- 
tron beam travels are uniform, i.e., not a function of 
vertical position. 

c. The fringe fields at the ends of the deflector are 
negligible. 

d. Energy transfer from the beam to the field is 
negligible, i.e., the beam density is small. 

e. Only the TEM mode is propagated down the 
deflector. 

f. The conductor and dielectric losses are negligible. 

g. The length of the deflector is sufficiently short and 
the magnetic forces sufficiently small so the effect 
of the magnetic fields on the electron displacement 
and on the longitudinal velocity may be neglected. 
Hutter [7] and Spangenberg [8] discuss the magnetic 
effects involved in this assumption. 

h. The electron beam radius is assumed to be small 
compared with the distance between the deflector 
conductors and compared with the displacement of 
the electrons due to the voltage applied to the deflector. 


3. Analysis of the Deflection System 


We desire to obtain a transfer function that relates 
the input voltage to the resultant deflection of the 
electron beam at the end of the deflector. Since we 
wish to find the entire effect of the deflection system, 
we will consider initially only the displacement in 
the deflector and take up the effects of any possible 
drift space later. We note that the electrons travel 
at a velocity less than the velocity of propagation of 
the electromagnetic wave so that an individual elec- 
tron slips behind the electromagnetic wave as the 
electron and the wave travel down the deflector. As 
a result, the system might be referred to as a slip- 
ping traveling wave structure. The slip will neces- 
sarily have to be accounted for in the final expression. 

A modified phasor approach will be used to find the 
oscilloscope transfer function. This involves finding 
the sinusoidal steady state output of the deflector 


y(t) =|H(w)|Vm cos [wt+ ¢(@) } (2) 
caused by the sinusoidal input voltage 
v(t) =Vm cos wt. (3) 


Once (2) is found, the system transfer function can 
be written directly since it is 


H(w) = |H(@) |e. (4) 


The output and the input are related through two 
physical phenomena, which occur simultaneously. 
These are the force and the resultant acceleration on 
the electrons due to the electromagnetic field and the 
effect of the traveling wave nature of the field. First, 
considering the force due to the electric field and using 
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FiGURE 3. Cross-sectional representation of deflector. 


the Lorentz force expression and Newton’s Second 
Law, we get 


Fy=eE (5) 


Teme (6) 


m 


where F; is the force in the y direction due to the elec- 
tric field E in the y direction, and y is the displacement 
of the electron in the y direction; e and m are the charge 
and mass of the electron. The coordinate system used 
is shown in figure 3. 


From the solution of the wave equation, we know that 
the electric field is a function of the position variable 
x and time ¢ and is of the functional form 


x 
v (:-= 
E(x, t)= q i 


for a uniform field, where d is the distance between 
the conductors, v is the input voltage and vp is the 
velocity of the electromagnetic wave, which is, through- 
out the discussion to follow, the velocity of light in a 
vacuum. The x displacement of the electron is a func- 
tion of the velocity in the x direction, vz. If we assume 
the electron velocity is constant, 


7 dx # 8 
2-2 = constant. (8) 


Then, 
x=vz(t—to) (9) 
where fo is the time an electron enters the deflector 


at x=0. If we substitute (9) into (7) and (7) into (6), 
we get 


d’*y e | 
t 


dt? md” 


Vr 
—— (ws) |: 
p 


If we specify the input voltage as a sinusoid, 
v(t) =V» cos ot, 


then (10) becomes 


d*y Ur 
We KV m cos w |e—22 (ea) | (12) 


'p 
where K=e/md. The initial conditions at x=0 are 


dy _ 


t= ts y=0, 7, =0. (13) 


By integrating (12) twice and satisfying the initial con- 
ditions, we find 


KV» be 
y(t) = Do? | CoS wtp — (t—to)wo sin wto 


— cos w (or+ to a] (14) 


Up — Vr 
Up 


o= 


(15) 


Sigma is a measure of the slip of the electron with 
respect to the electromagnetic field. Commonly, elec- 
tron velocities are described by a ratio B=vz/vp. The 
slip parameter o is seen to be c= 1—£. Examination 
of (14) shows the electrons follow a straight-line path 
with a sinusoid superposed, oscillating at an angular 
frequency w'=ow, the slip frequency. If we could 
make B approach one by increasing vz, the slip fre- 
quency would approach zero, and we would have a 
true traveling wave structure as described by Owaki 


[1] and Talbot [3]. 
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Equation (14) is plotted in figure 4. It shows the 
normalized displacement 


if . 
KV nlaro? COS Wlo — (t—t)wo sin @lo 
m 


v 
—cos w (or+ to e=) 
Up 


plotted versus w(t—to). Each curve represents the 
path an electron follows after entering the deflector 
at the indicated value of wto. Since a practical deflector 
would only be, at most, a few centimeters long, i.e., 
fractional nanosecond transit time, the oscillatory 
character of the electrons paths would normally not 
be observed except for extremely large w. 
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FIGURE 4. Electron path within deflector versus w(t—to) for sinu- 
soidal input. 


Transit time of the electron, 7, is given by 


tT— 
V;’ 


thus t=t)»+7 at the end of the deflector. If we solve 
(14) for displacement at the end of the deflector, we get 


( KV» [ ’ 
(to) = ~(LCOS Wlo— TWO SIN wt 
y 0) oo” 0 0 


—cos (wo7+ wto) }. (16) 
Equation (16) shows the displacement at the end of 
the deflector is a sinusoidal function of the time the 
electrons enter the deflector, t. Since we desire the 
electron displacement as a function of changing to, 
let us set 


to=t' (17) 


where t’ is a variable time at the end of the deflector. 


Replacing t by t’ yields the sinusoidal steady state 
response expressed as follows: 


KV», ’ e ’ 
—— [cos wt’ — wre sin wt 


a aah 


—cos (wt'+ ata) }. 


(18) 


If we write (18) in phase and magnitude form, we 
will have the deflector output in the proper form for 
finding H(w) as specified in (2) and (4), providing we 
properly account for the transit time delay between 
input and output. The identity 


cos (wt’ + ata) =cos wt’ cos wra— sin wt’ sin wro 
(19) 


allows us to write (18) as follows: 
y(t’) as [(1—cos wro) cos wt’ 


+(sin wra— Two) sin at’ | (20) 


which may be written as 


y(t!) = VA Bos (wt’+¢), 


wo” 
where 


A=(1—cos wro) 


B=(sin wto— Two) 


g=—tan" ri 


Equation (21) provides the necessary information for 
writing the oscilloscope transfer function except for 
the problem of accounting for the delay. Without the 
delay factor, we get from (21) 


H'(w) 2 V A? + B? e#, 


wo" 


(22) 


The delay is provided by introducing the delay factor 
e~J" giving the complete oscilloscope function 


K =, 
H(w) = VA? + B* elle-09, 


(23) 
where A, B, and ¢ are specified in (21). 
4. S-Domain Transfer Function 


While it is true that (23) is the desired oscilloscope 
transfer function, it offers little insight into character 
of the system. The S-domain (Laplace) transfer func- 
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tion, which can be derived from (23), is more compact 
and is in a more familiar Laplace transform. If we write 
(23) in rectangular form, with the delay term left in 
exponential form, it becomes 


H(w) = ee" (4A —jB) 


wo” 

K'e~ie" me 
=—.— [(1—cos oro) —j(sin ora — Two) J, 
2 


(24) 


@ 


K 
where K’ >, Writing the trigonometric terms of (24) 


in exponential form, we get: 


K'e-Je7 ( ejoro + e jutro 
i) 
w? 2 


(== sain e jure 


—j Fj _ rw) | (25) 


H(@) = 


If we set w=s/j and simplify, (25) becomes 


H(s)=K'e*| 3 


a 





H(s) =K' | (26) 


e7 sr l-o) — Tose~** — e~ 7 
s? 


Equation (26) provides immediate information on 
the response of the deflector to various inputs. For 
example, if the input is a step function, the output is 
the sum of two parabolic terms and one linear term. 
One of the parabolic terms and the linear term are 
both delayed by 7 seconds after the start of the step 
input; the other parabolic term is delayed (1—oa)r 
seconds after the step input. This demonstrates a re- 
sult which physical reasoning indicates, i.e., that the 
output response starts at 


which is the transit time of the electromagnetic wave 
and not t=7, the transit time of an electron. 

If we set //v»=7', an alternate form of the s-domain 
transfer function may be written, giving 





2 


e 
H(s)=K'| 


8’ — rorse—87 — <<) 


; (27) 


5. The Impulse and Step Function Response 


The impulse response of the deflector can be found 
by either finding the inverse Fourier transform of (23) 
or by finding the inverse Laplace transform of (27). 


The latter method leads directly to 
h(t)=K' [((t—7')u(t —7') — ro u(t — tr) —(t— 7) u(t — 7) 
(28) 


where u(t—7) is the unit step function delayed 7+ 
seconds. 

In the previous section, we discussed qualitatively 
the system step function response. Now, we would 
like to write an expression for the response to the step 
input voltage, v(t)=Vu(t), by taking the inverse 
Laplace transform of 


Y(s) =V(s)H(s) (29) 


V 
V(s) am 


and H(s) is as specified in (26), which yields 


y(t) =Kv | 


u(t—7') —ta(t—T)u(t—7) 


-£o we-2)| (31) 


as the step function response. Figure 5 shows a sketch 
of the step function response. 


KVr2 
2 





-7')? 
ae Ly u(t-r') 








FIGURE 5. Step function response of deflector. 


Equation (31) shows that the display at the end 
of the deflector (due to a step function of voltage 
applied at the deflector input) starts at t=7’, i.e., 
the transit time of the electromagnetic wave, and 
increases parabolically until t=7, the electron transit 
time. After t=7, the displacement becomes the sum 
of a positive going parabola, a negative going parabola 
and a negative going ramp. This sum term is simply 
KV7?/2, which is the steady state response and is, 
as would be expected, the deflection due to a d-c 
input voltage, V. The slip parameter o=1—vz/vp 
appears in the transit part of the response and cancels 
out in the steady state part. This is as would be ex- 
pected since the concept of slip applies to a traveling 
wave phenomenon and not to a steady d-c voltage 
applied to the deflector. 
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6. Frequency Response Function 


The frequency response function for the deflec- 
tor is simply H(w) as expressed in (23) where the 
amplitude response is 


K 


ow" 





|H(w)|= 


V (1-—cos wra)?+ (sin wro — wt)? 
(32) 


and the phase response is 


@TO — SiN W710 
1—¢os wto 





¢(w) =— r+ tan“! (33) 


The linear phase term w7 simply delays all components 
of the signal an amount 7 and causes no phase dis- 
tortion; therefore, it is omitted in the plot shown in 
figure 6, where the 


normalized amplitude = |H(w)|/|H(0)| (34) 


and the phase g(w) are plotted versus 


normalized angular frequency = wro. (35) 
The bandwidth of the deflector may be defined as 
the frequency where the normalized amplitude drops 
to 0.707, which occurs at wro = 1.127, which gives 


— @.707 _ 1.127 _ 0.56, 

Sao 2x 2ntTe To (36) 

This result shows that the bandwidth is inversely 
proportional to o7 where ' 


ee 3 (1 —#)=1 (2). 
Ur Up Vr Up 
which indicates that the bandwidth is proportional 
to the electron velocity, vz, and inversely propor- 
tional to the length of the deflector, /. In the next 
section, we see that exactly the inverse is true for 
deflector sensitivity. 

Actually, the concept of bandwidth is of little con- 
cern in this discussion since with closed expressions 
for |H(w)| and g(w), we may determine exactly the 
contribution of harmonic components to as high a 
frequency as we please. Notwithstanding, since 
bandwidth and its counterpart, rise time, are the usual, 
but approximate, figures of merit used to describe 
oscilloscopes, we will give them suitable emphasis. 


7. Design of a Practical Deflector 


In designing a practical deflector of the type shown 
in figure 3 with a calculable transfer function, we must 
consider the length, /, conductor spacing, d, and the 
velocity of the electrons, vz. Since a d-c voltage 
, gives the maximum deflection, as shown in figure 6, 
it will be used as the basis for design. The d-c deflec- 
tion is the steady state value of (31), which is 


y (ve) = Aire, 


(37) 


where V is the amplitude of the d-c input voltage or 


KV I? eV I? 


WO) =F = Sd of (38) 


where m is the relativistic electron mass, mo/V 1 — B?. 
The sensitivity then, is 


Deflector Sensitivity = 





we) _o Vi-# meter, (39) 


2modv*. volt 

The symbol mp is the rest mass of the electron, and 
B is the ratio of the velocity of the electron to the 
velocity of light, which is the same as the B discussed 
in the paragraph following (15) for’ an air-filled strip- 
line. We must consider the effect of relativity since 
we want to minimize the slip, i.e., use relatively large 
electron velocities, in order to increase the bandwidth. 
This, of course, decreases the sensitivity since (39) 
shows that the sensitivity is directly proportional to 
the length squared and inversely proportional to the 
electron velocity squared if the effect of B? is ignored. 
It also shows that a compromise is necessary between 
length and conductor spacing in order to optimize 
the sensitivity and avoid the problem of electrons 
striking the wall of the deflector before reaching the 
end. Some practical upper limit of vz and some usable 
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FIGURE 6. Amplitude and phase response of deflector. 
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———— f= 2.125" 
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FIGURE 7. Slipping traveling wave deflector (r = 429 ps for B= 0.42). 


lower limit to the deflector sensitivity are the govern- 
ing criteria in the design, although the other oscillo- 
scope criteria such as light-output intensity, writing 
speed, and size of screen display are also of primary 
importance. 

An example of a usable deflector, where Lee’s [4] 
microoscillographic techniques must be used,’ is 
shown in figure 7 which is designed for an accelerating 
potential of 50 kV, i.e., B=0.42. The structure is de- 
signed with a characteristic impedance of 50 1 based 
on Metcalf’s [9] design curves of characteristic im 
pedance for the rectangular strip lines. The calculated 
deflector sensitivity from (39) is 9.64 10-§ m/V or 
about 0.01 mm/V. The lowest TM mode cutoff wave- 
length is approximately twice the ground plane spac- 
ing according to Howe [10] 


Acorm = 2(0.15) = (0.30) in= 0.662 cm, 
which is equal to a cutoff frequency of 39.3 GHz. 


The lowest TE mode cutoff wavelength equals approxi- 
mately the mean circumferential path 


=0.53 in= 1.346 cm 


which is equal to a cutoff frequency of 22.3 GHz. 
The bandwidth of this design from (36) is 


0.56 _ 0.56 
fa =e = 929 X 10-2 X0.48 





= 2.82 GHz. 


It is interesting in this case, to calculate, using the 
derived transfer function, the actual rise time of the 
deflector to a step input. This is easily done using (31). 
Equation (31) leads to the oscilloscope display shown 
in figure 5, which has a 10 to 90 percent rise time 


tr=OT V0.9 —OT V0.1 = 0.63407 = 148 ps. 


8. Drift Space Considerations 


Thus far in our discussion, we have considered an 
oscilloscope with no drift space beyond the deflector. 
In general, some drift space will be involved; if we 
include a drift space, the oscilloscope transfer function, 
(27), must be modified accordingly. As shown in the 
appendix, the total oscilloscope transfer function in 
the S-domain is 


*Lee [4] describes the nature of the appropriate beam size and focusing techniques 
as well as suitable methods of recording for a microoscillographic deflector. 


7'K e78(t' +7") = e+"), 
S (A-8) 





Hr(S) =H (S)e-"s+— 
o 
where 7” is the transit time of the electrons in the drift 

space 


where L is the length of the drift space. The drift space 
has the effect of increasing the sensitivity. D-C sen- 
sitivity, for example, is increased by 2L/I < 100 per- 
cent. Thus, if L =/, the deflector sensitivity is increased 
by 200 percent. 

In the appendix the expression for the relative sen- 
sitivity of the slipping strip line deflector is developed 
for the drift space contribution of displacement only, 
i.e., the deflection within the deflector space is ignored. 
The result is shown in (A-—12). 


in 227 
[Y¥(w)|_ "2 
Y(0)| ow 
2 
It is interesting to compare this result with Holl- 


mann’s expression for parallel plate deflector given 
in (1). 





(A-12) 


. @T 
¥(w)|_ sim z 
[Y(0)| wr 
2 





(1) 


We see that the results are identical except for the 
slip factor o. 

Hollmann’s result is based on the assumption of a 
space invariant electromagnetic field within the de- 
flector while (A—-12) includes the effect of a slipping 
traveling wave. We see, since 


Up — Ur 
CS 
Up Y 


(15) 
that under Hollmann’s assumption, i.e., infinite wave 
velocity, the slip approaches one and (A-12) ap- 
proaches (1) as it should if this work and Hollmann’s 
results are to be compatible. 


9. Conclusion 


Knowledge of the oscilloscope transfer function 
makes it possible to predict the response charac- 
teristics of the oscilloscope to the extent that the 
physical parameters in the expression can be measured 
and to the extent that the assumptions employed lead 
to errors in the calculated transfer function. Uncer- 
tainties of these types can be held to a minimum by 
careful construction and measurement methods. Also, 
the accuracy of the transfer function can be approxi- 
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mately verified experimentally from sinusoidal meas- 
urements. A thorough error analysis would contribute 
to knowledge of the magnitude of uncertainties in the 
oscilloscope transfer function. 

Once the transfer function and its uncertainties 
are known, it can be used to calculate the actual time 
function of either a periodic or nonperiodic input pulse 
from a knowledge of the output display. An analog- 
to-digital conversion of the output can be used for 
fast and accurate machine calculations. An instrument 
of this type, which requires microoscillographic tech- 
niques and machine calculations, is not practical for 
normal laboratory measurements of pulses. It can be 
used as a standard oscilloscope for measuring pulse 
parameters such as rise time, overshoot, duration and 
so forth. It also may be used to determine the input 
voltage amplitude as a function of time to a stated 
accuracy. 

A significant result is found by comparing the rela- 
tive sensitivity of the slipping traveling-wave strip-line 
deflector (for the drift space only) with Hollmann’s 
[2] results for the parallel-plate structure. It is con- 
cluded that the strip-line results provide a more general 
expression for the parallel-plate deflector than 
Hollmann’s expression. 


The author acknowledges the many helpful sugges- 
tions, counsel, and encouragement received from 
A. R. Ondrejka, and the careful manuscript prepara- 
tion by Mrs. Toni Hooper. 
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11. Appendix 


The addition of a drift space whose length is L 


meters introduces an additional electron displace- 
ment, Y, such that 


Y=y'r" (A-1) 
where 7” is the time the electron is in the drift space, 
i.e., drift space transit time, and y’ is the vertical 
velocity of the electron at the end of the deflector. 
Assuming the x velocity is constant, 


(A-2) 


7 ak. 
Ur 


The total deflection, yz, is 


yr=Y+y. (A-3) 

Once the electron leaves the deflector, the dis- 
placement force due to the traveling electromagnetic 
field stops, i.e., y becomes constant. Thus to find 
yr, we need only to find the expression for Y. This 
requires that we obtain an expression for the vertical 
velocity at the end of the deflector. The vertical 
velocity within the deflector is found by integrating 
(12) and evaluating the arbitrary constants using 
(13). The vertical velocity at the end of the deflector, 
y’, is found in the same manner as (18) was found, 
giving ; 


1» _ KV im . > a ~! 
ps [sin (wt’+wor)—sin wt'] (A-4) 


where t'=t—(to+7) for t2=to+7. Equation (A-4) 
can be written in phase and magnitude form giving 
y =A VOTB cos (wt’—9)] (A-5) 

where 

C= (—1+ cos wor) 

D= sin wot 

©= Tan-' — 1+ cos wot 

sin WoT 


The transfer function due to the drift space is then 


Hy(w) =~ Ver+ De® (A-6) 
Wo 


which must be combined with (23) to give the total 
transfer function: 


Kr" 


Hr(w) = H(w)e-”" + —  V C2 + D2 eO+ 07" +7) 
wo 


(A-7) 
where e~” and e~%" provides the delay caused by 
the drift space, and deflector respectively. 


The S-domain transfer function may be found by 
setting w=S/j in (A-7) which can be shown to yield 


ar'K oo — e-(r+1")s 





Hr(S) = H(S)e-"* + — 
o S 


(A-8) 
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While (A-8) is the transfer function for which we 
are searching, it is interesting to compare this re- 
sult for the strip line deflector with the work of Holl- 
mann for the parallel plate deflector, as expressed 
in eq (1). The comparison can be made by consider- 
ing only the Y deflection since Hollmann’s equation 
includes only this part of the total deflection. The 
comparison is made by finding the relative sensitivity, 
as defined by Hollmann (which is |¥(@)|/|¥(0)|), 
for the strip line. 


The expression for |Y(w)| may be found from 
(A-6) since 


How) = 


Y(w) = Hp(w)V(). 


(A-9) 


Substituting the magnitude of (A-6) into (A-9) gives 


Vink e 
l\Y(w)| = a V(—1+ cos wot)? +(sin wor)? (A-10) 


7) 





since V (w)=V ne, i.e., the input voltage is Vm cos wt. 
It can be shown that 


Y(0) = KV prt" (A-11) 
which is the displacement due to a d-c input, Vm. The 
desired ratio is then found from (A-10) and (A-1]), 
giving 





[Y(w)| _ V(—1+ cos wor)? (sin wor)? 
|Y(0)| 





WOT 


9 [1 —cos wor 
lY(w)| _ 2 


lY(0)| 





WOT 


or finally 
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1. Introduction 


This paper describes a mathematical analysis for 
determining the value of the substitution error of a 
bolometer with a Wollaston-wire element (barretter). 
The analysis reflects all significant nonlinearities in 
the heat flow, including some not covered before, and 
includes all appreciable heat transport mechanisms 
simultaneously. The microwave current distribution 
has been measured in a scaled model and found to be 
approximately sinusoidal in shape with a current 
maximum in the middle of the wire. Input data are 
chosen so as to represent the operating bolometer 
characteristics as nearly as possible. 

The values of substitution error thus obtained, in 
conjunction with efficiency data obtained by micro- 
wave techniques [1],' will be very useful in extending 
power meter calibrations to frequency ranges where 
extremely accurate microcalorimeters are not avail- 
able. A method of measuring substitution error in- 
volves two distinct measurement techniques [1], 2], 
and since agreement is quite good between the value 
of substitution error thus measured and the value 
calculated from this analysis, the two provide a cross- 
check for each other. 


2. Background 


In the measurement of microwave and millimeter- 
wave power (hereafter referred to as rf power) by 
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substitution techniques, any substitution error is of 
great concern. A common measurement technique is 
to replace.a known amount of d-c power with an 
unknown amount of rf power in a bolometer. Un- 
fortunately, the different current distributions generate 
different temperature fields which give the bolometer 
element slightly different values of total resistance 
for equal amounts of power. For balanced-bridge 
methods the resistance is maintained constant, which 
causes some non-equivalence of power. In this case 
Carlin and Sucher [3] defined substitution error as 


E = W sus = W vt 

Wa 
where W,, is the rf power dissipated in the barretter 
element and Wsgyg is the substituted d-c power which 
yields the same value of barretter resistance. Although 
a bolometer can be used to measure relative powers 
without calibration, it must be calibrated to measure 
absolute power. An accurate method of calibration 
is by the use of a microcalorimeter [2], in which the 
ratio of Wgypz to the net rf power flowing across an 
arbitrary plane Wygz into the bolometer is measured. 
This quantity is called the effective efficiency, 
Ne=Wsypl/Wner. Another calibration method (for 
barretters only) is the impedance technique [1], in 
which the ratio of W_,to Wygr is measured. This quan- 
tity is called the efficiency, n= Wy4/Wyer. By replacing 
Wu, and W,, in the equation for E, an equivalent form 
is given by 
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This provides a way of measuring substitution error, 
but with the substantial limitation that in typical cases 
this necessitates subtracting two nearly equal numbers, 
and the relative uncertainty of the result is very large. 
Some results, comparisons, and comments will be 
given later. 

Another equation derived from the previous ones 
is given by 
Wer = SB : 

“  M1+E) 

This form provides a method of measuring the net 
power from Wgyg, n, and E. The substituted d-c power 
is measured relatively easily, eficiency can be meas- 
ured by the impedance technique (as mentioned 
previously), and now with the results of the analysis 
of this paper, substitution error can be calculated. 
For frequency ranges where microcalorimeters are 
not available, or for organizations to which they are 
not accessible, this provides an accurate means of 
calibrating barretters. 

Previous work on determination of substitution 
error has been done by Carlin and Sucher [3], Weber 
[4], and Bleaney [5]. However, none of the previous 
workers attempted to solve the total heat-flow problem 
or to take into account all of the nonlinearities. Carlin 
and Sucher have determined an upper bound of sub- 
stitution error from a combination of solutions to 
parts of the total problem. Weber has handled one 
nonlinearity in a solution to part of the problem. 


3. Description of the Barretter 


The actual barretter as shown in figure 1 can be 
represented by the model shown in figure 3. Partial 
justification for this is based on results of a model 
that used electric potential to represent temperature 
potential (fig. 2). These results showed that flush sup- 
ports in a model yield negligible difference in potential 
from other common support structures. This is im- 
portant, not only because flush supports give the 
simplest mathematical model, but also because there 
are large variations in the shapes and dimensions of 
support structures, even for the same models by the 
same manufacturer (fig. 4). 

The rest of the justification for using the model 
shown in figure 3 becomes clear in the discussion of 
operation .and different heat transport mechanisms 
that follows. 

The barretter is a small-diameter, platinum wire of 
radius a and length / (fig. 3). It is supported at each 
end by a silver wire of larger diameter, and is normally 
immersed in air. Possible external, forced-convection 
currents are eliminated by a plastic shield of radius } 
which is nominally transparent to the electromagnetic 
fields. 

For power measurements, this fine wire is usually 
heated with a d-c current to a particular value of 
resistance in the absence of rf fields, and then the d-c 
current is reduced to maintain the previous value of 
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FIGURE 4. Various common support structures. 





resistance when the rf field of frequency f is applied. 
The substituted d-c power is (/? —/2)R, where J, is the 
d-c current without rf power, /> is the d-c current with 
rf power, and R is the value of resistance to which the 
barretter is biased. The substitution error is caused 
by the difference in temperature fields between the 
two heating conditions. The solution of these two tem- 
perature fields will then require the additional data 
concerning thermal properties of the wire and the 
surrounding fluid, usually air. The thermal conduc- 
tivity of the wire, with its temperature coefficient, 
may be represented by ki=kio+hi:7. Similarly the 
thermal conductivity of the fluid may be represented 
by ko= kx + koT. The resistivity of the wire and the 
corresponding temperature coefficient may be repre- 
sented by p= pot pil. 

The previously mentioned parameters could all 
be input data; however, as is mentioned later, it is 
more judicious and accurate to arrive at some of these 
parameters from other measurements taken under 
operating conditions. 


All heat-transfer mechanisms that can possibly act 
in a barretter will have to be considered and evaluated. 
Of conduction, convection and radiation, only conduc- 
tion turns out to be significant. Radiation is negligible 
since surface area of the wire is so very small. Power 
radiated is given by 


p=oSF(Ti—T%), 


where p is power in watts, 
o is the Stefan-Boltzman constant, 
S is the surface area of the fine wire, 
F is a surface emissivity factor and the 7’s are 
absolute temperatures in °K. 


The diameter of the wire is of the order of magnitude 
10-4 cm, and typical calculations show the radiated 
power to be less than 0.01 percent of the total power, 
even if the wire were at its melting temperature. 
Convection may also be shown to be negligible (0.01%) 
by evaluation of the Grashof number which is some- 
what proportional to convective heat transport. This 
may also be confirmed by varying the orientation while 
measuring d-c power. The Grashof number is 


+ Bop? L*( T;,—T)) 
7 





‘Yar 


where B is the cubical expansion coefficient, 
g is the gravitational constant, 
p is the density, 
L is the width of the gas space inside the plastic 
convection shield, 
n is viscosity, and 
T’s are the absolute temperatures. 


The Grashof number thus calculated for typical 
bolometers is less than 10; if it is less than 1000, 
convection may be considered negligible. 

The only remaining heat transport mechanism, 


Sliding dc Ammeter 
“-~Bolometer 


FIGURE 5. 8X 16 in Scale model of waveguide bolom- 
eter for measuring rf current distribution in 
element. 


conduction, can be through both the fine wire and the 
air (or other fluid). The thermal conductivities of both 
media are functions of temperature over the rather 
wide temperature range for normal operation. 


It is essential to consider the current distribution 
in the wire. The d-c current distribution is considered 
uniform both axially and radially. The rf current is 
essentially uniform radially, even to frequencies as 
high as 100 GHz (< 10% variation), primarily due to the 
smallness of the wire. The axial rf current distribution 
is of greater uncertainty, and has been investigated 
extensively with a scale model (fig. 5). Details are 
discussed in the following three paragraphs. 


The measured values indicate a current maximum 
near the middle of the wire and in general a sinusoidal 


_or parabolic shape. There is a frequency dependence 


on the amount of curvature. There is great variation 
in current magnitude but only slight variation in shape 
with respect to relative position of the shorting end- 
plate and with respect to position of tuning screws. 


Some variation in the location of the current maxi- 
mum may be introduced if the shorting end plate is 
tilted only a few degrees from perpendicular. The 
shape and size of the wire support structure cause 
variations in the current shape, but to a large extent, 
only insofar as the wire length itself is varied. An 
extreme case not encountered in most bolometers is 
when there is no support structure — the wire supported 
only by the wide walls. In this case the current dis- 
tribution is nominally uniform, with great susceptibility 
to asymmetry introduced by varying angles in the 
shorting end plate. 


The measured values seem to follow a sinusoidal 
distribution with wavelength greater than free-space 
wavelength and with some asymmetry (< 30° for com- 
mon support structures). The sinusoidal distribution 
has been used by previous workers [3, 4, 5], but effects 
of asymmetry have not been previously considered. 
Further investigation is being considered, but for 
present calculations rf current is assumed to be sinus- 
oidal with a free-space wavelength. 
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There are several pitfalls in determining three of 
the input parameters a, po, and py. 

First, the radius a or diameter d cannot be 
measured directly with a microscope without 
bringing the objective lens destructively close to 
the barretter if sufficient magnification is used. 

Second, the small diameter of the wire causes 
po to vary from the bulk properties of the metal 
[6] since the wire diameter approaches the mean 
free path of an electron; p; does not vary this 
way [7, 8]. 

Third, a set of three equations with a, po, p: 
unknown is homogeneous and the determinant 
vanishes; it is possible to solve only for two of the 
unknowns in terms of the third. 

The second circumstance indicates that if a value 
for one of these three unknowns must be assumed, 
assuming a value for p; causes the least uncertainty. 
The equation 


Ria? — pol = piTil, i=1, 2 


can then be solved for a and po if two sets of R and T 
are given; / can be measured accurately. One of the 
R and T sets is for room temperature while the other 
is for operating temperature. The average operating 
temperature is calculated by the computer program, 
as mentioned later. 


4. Analysis 


In this section, we give a description of the analysis 
used to obtain the substitution error and of the approxi- 
mations involved. Details are given in appendix A, to 
which reference is frequently made. 

Referring to the model shown in figure 3, the basic 
mathematical problem is to solve two distinct cases 
of the heat flow equation 


V (kot hT)V T=S(r, x, T) (1) 


for the temperature field 7(r, x), where (ko, ki) are 
conductivity coefficients (different for r<a, r> a) 
and S is the local power dissipation. (A—2, A—3) 

In the first case, the temperature distribution 7; is 
generated by the source S, (A—1) due to dissipation of 
the d-c current 4, flowing uniformly in the wire with 
resistivity linearly dependent on 7. 

In the second case, the temperature distribution 
Ty is generated by the source S, (A—4) due to dissi- 
pation of a d-c current .4., as above, and also due to 
the time-average dissipation of an rf current (A—9a): 


F(x, t)= Fro2z(x) cos wt, (2) 


where zx), the current profile, is given, and the 
amplitude “4x9 is an unknown constant. 


Boundary conditions on the temperature field 7(r, x) 
in each case are: 


T, (0, x)=T(b, x)=0, 


(0<x</l); 


T(r, 0)= T(r, )=0, 
{T(r, x, Fs Ar, x)} 


continuous across r=a, (0 <x </). 

The rf current amplitude Zo is to be determined 
by requiring that, given .7; and #2, the total wire re- 
sistance Ryo, be the same in each case (A—20): 


l l 
Rrot= [ dx(po + piT(O, x)). (4) 


Because of the nonlinearities, the problem posed 
above is not tractable, and several simplifying assump- 
tions are introduced which lead to a linear system of 
equations involving a single parameter, the mean 
axial temperature: 


ae l 
7\=1 [ dxT\0, 2). (5) 


The linear system is solved iteratively and converges 
rapidly. 

Before making the simplifying assumptions, it is 
useful to introduce the “‘quasi-temperatures” (A—5): 


G(r, x)= T(r, x) (1 re hy T\(r, »)), 
2 ko 


H(r, x)=Ty(r, x) (1 re hy T(r, »)), 
2 ko 
g(r, x)= H(r, x) —G(r, x). (6) 


These quantities simplify the heat flow equations, 
since 
V . (ko + kiT)VT, = koV?G, 
(7) 
V (kot ki TyWVTy= koV7H. 


The interface conditions at r=a are, of course, more 
complicated: 


{hot  - (1 ++ (Z G) ) 8) 


are continuous across r=a, (0 <x </), with identical 
conditions on H. Here 7(x) is the inverse function 


(A-7b): 
— 
G (+7(2 c) )=n, (9) 
ko 


The three basic simplifying assumptions that are 
made are the following: 


First, the difference function g is assumed small: 


la(r, x)| < |G(r, x)|. (10) 
This approximation appears to be adequately justified 
by the fact that computed values indicate |g| < 0.15|G| 
at worst, while the axial average g=0. Neglecting 
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O(\g|*), we obtain equations of the form (A—-8a, A-11): 
V26 + A,G[1 + ¢1(G)]= Ae, 
V2g + Asg[1 + g2(G)] = (x) 


: [po > piG(1 + ¢3(G))], (11) 


where the A; are constants (different for r >a, r< a), 


ki, 
the functions gi(G) depend on 7 e C), and for cases 
0 


studied, we have again |gi(G)| < 0.15. 

The function ¢€(x) depends on the difference in current 
flows in the d-c and rf cases. The second (nonlinear) 
interface condition is (A-14) 


IG. + ga(G)| 2 =0, 


(12) 
|g(1 + ¢5(G))|2: =9, 


with g(G) as above. 
The second simplifying assumption is that for 
0 <r<a, we may put 


giG)= pi(GO, x)) (13) 


using the axial average, while for r >a, we may put 


gilG) = p(G*) (14) 


where G* satisfies the interface condition in the mean 
G*(1 + glG*)) = G1 + gx). (15) 


This approximation makes use of the expected flatness 
of the temperature distribution in x, and is clearly 
better for large //b. Numerical tests indicated the 
results were very insensitive to the choice of G; a5 per- 
cent change in c gave a 0.06 percent change in E ina 
typical case. 

The third simplifying assumption, used indirectly 
above, is that we may put 


G(r, x)= GO, x) (16) 
in the source term for g, (r<a). This is readily justi- 
fied by the extremely small radial temperature drop 
in the wire. 

With these simplifications, the linear equations with 
parameter C© in (g, G) (A-12, A-13, A-17, A-18) are 
readily solved in terms of cylinder functions. With 
the same simplifying assumptions, we find the balanced 
bridge requirement (A—21): 


g(0, x)=0 (17) 
which determines the unknown rf current amplitude 
Fro implicit in €(x). oe 
After iteration to determine the parameter G, the 
required power levels Wgyg and W,,, the substitution 
error E (A-27), and the temperature distributions 
T,(0, x) and AT = T,,(0, x) — T,(0, x) (A-34a, A—34b) are 


readily found by numerical integrations involving 
G(0, x) and the rf current profile z(x). 

Details of the computer programs are available from 
the authors. Further details for the sinusoidal rf 
current distribution are given in (A-37, A—-38, A-39). 


5. Uncertainties 


With given input data, the analysis and computation 
described in the previous section is believed to yield 
a substitution error which is correct for the cases 
studied to within a few percent, though it has not 
been possible to construct an error bound. 

An evaluation of the total differential reflecting/the 
uncertainties of all the input data would seem to offer 
a limit of uncertainty to the calculated value of substi- 
tution error, so that one concerned with applying this 
technique would have sufficient justification for doing 
so. This figure is nominally 20 to 25 percent, which is 
quite adequate for the typically small values of sub- 
stitution error. 

Unfortunately the shape of the rf current distribu- 
tion (assumed sinusoidal with a free space wavelength 
and with a current maximum in the middle) is known 
no better than it can be measured with a scale model. 
Uncertainties due to variation of wavelength and shift 
in position of current maximum are partially reflected 
in this total differential; phase shift variations cause 
only 0.2-0.3 percent decrease in substitution error, 
and are included in the above uncertainty, but the 
uncertainty as to wavelength could cause 10 to 40 
percent change either way and is in addition to the 
above uncertainty. Even this is just a best-judgment 
estimate. 

This degradation could be substantially reduced if 
a sound theoretical analysis of the current distribution 
in the barretter were available. Since none is, further 
investigation is certainly in order. 

As to the input data, R.ojq and Ryo may be measured 
quite accurately, while only two of the four quantities 
Po, 4, Reoig and Ryo are independent. First, the follow- 
ing relation must hold: 


_ pol 
Reoia .” ry 
Ta 


(18) 


A second relation results as follows: Given the param- 
eters a and x;={Reoia, pr, ky, 6, 1, 42, F3, f, tan og} 
the temperature field 7(x, 0) is determined 1n the wire, 
as is the heated wire resistance R 


R= we (po + pil (x, 0)) 
Ta 


=R(a; xi) 


and a substitution error is determined 


E=4(a; xj). 
An iteration is made to find a* such that 


R(a*; xi)= Rrot 
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which is the required second relation. Then also 


E=&(a*; xi) (22) 
is the proper substitution error. 

These forms are also used to determine the effect 
on E of small uncertainties in the data x;. From small 
changes in (xj, a), we compute the partial derivative 
approximations: 

0 AOR da& OR 
ae og eR (23) 
da da Oxi Ox; 


With the bridge balanced 


oR aR 
bRpot = — 6a+ i 6x;=0 


while 


aR 
0€ 0€ OXi 


OXi da OR 
da 


gives the error in E due to small input error 6x;. 


6. Results, Comparisons, and Evaluation 


Values of substitution error have been calculated 
for several common barretters and where possible 
comparisons are made with some measured values 
(see table 1). Although the agreement is reassuring, 
the uncertainty in the measurement makes the com- 
parison inconclusive in itself. 

A number of variations of input data lead to several 
interesting observations concerning the calculated 
values of substitution error. 

Based on calculated results which use sinusoidal 
current distribution with free-space wavelength, 
substitution error varies approximately as the square 
of the frequency. 

The effect of asymmetry in the current distribution 
(induced by asymmetries in the physical structure) is 
to slightly reduce the substitution error, typically of 
the order of magnitude of 0.1 to 0.3 percent. 

The use of helium, a gas with higher thermal con- 
ductivity than air, greatly increases the power han- 
dling capacity of the barretter. It also increases the 
heat “coupling” to the surroundings and thus reduces 
the substitution error approximately 10 percent. The 
same effect is obtained (but without increasing the 
power handling capacity) by increasing the bias cur- 


rent (e.g., so that the resistance increases from 200 
to 267 2). 

A change in rf power level causes virtually no change 
in substitution error. 

The nonlinearity due to the temperature dependence 
of thé thermal conductivity of the air causes an ap- 
proximate 15 percent increase in substitution error 
over what would have occurred if that nonlinearity 
had not been present. The temperature dependence 
of the thermal conductivity of the metallic element 
causes only slight variation (~ 1%) in substitution 
error as only about 5 percent of the heat escapes the 
element by metallic conduction. 


TABLE 1. Comparison of measured and calculated values of 


substitution error 








Measured 
substitu- 
tion error 


i. 


NBS 
barretter 
unit 


Effective 
efficiency 
(ne) 


Calculated 
substitu- 
tion error 


Efficiency 
(n) 


Freq 
(GHz) 





0.9920 + 0.002 
.9709 + 0.002 























It is also of interest to compare results with those 
of previous workers, in particular Carlin and Sucher 
[3]. Although they obtained an upper bound for sub- 
stitution error, and this upper bound is typically 
20-25 percent higher than our calculated results, a 
closer examination of their work points to a flaw in 
one of their conclusions, namely that nonlinearities 
have little effect in these cases. The margin of safety 
in their upper limit comes from the fact that they solve 
two related linear problems, calculate the substitution 
error for each case and add the results together to 
obtain a safe upper limit. It happens that both of these 
related problems resemble the real problem, their 
equation (15), quite closely, and the results of their 
equation (15) should be a good solution to the linear 
approximation to the real problem. The corresponding 
results of equation (15) are about 40 percent lower than 
our results in most typical cases; this indicates the 
magnitude of the error made in assuming that non- 
linearities can be neglected. 


Appendix A. Analysis 


A wire of radius a and length / terminates on me- 
tallic surfaces at z=0 and z=/, the whole enclosed 
in a shield of inner radius 6, as shown in figure 3. The 
wire is heated by current sources uniform in 
r, (0<rs<a), and the temperature 7 above ambient 
is assumed to vanish on the shield r=b, (0<x</) 
and on the end surfaces z= (0, /), (O<r<b). 

Let 7, be the temperature field in the d-c-only case, 
in which the source strength S, in the wire due to d-c 
current .F; is 


S,= F?(pot+piT) (A-1) 
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where (Oo, p1) = (po0/A?, 
sectional area of the wire. 
The heat flow vector is 


p:/A?) and A is the cross- 


qi =— (kiot kinT,) VT, (A-2) 


and governing equation for heat flow is 


V-qa=S (A-3) 
with subscript i referring to O<r<a and a<r<b 
(fig. 1). 

Let 7, be the temperature field in the rf-substitution 
case. The source strength S,(x) is taken to be a time 
average due to a d-c current #2 and an rf current 


SF r(x, t): 


Su = ([%2+ Fr (x, t)]}?)(potpiTy).  (A-4) 


Equations analogous to (A-2, A-3) hold for the heat 
flow in this case. 


The fields (7, 7) and the heat flow vectors (q;, qu) 
are to be continuous across r=a(0<x</). 
Let us introduce new variables 


1 
=e 7,(1+5 nm) 


H(r, x) = Tu(1+5 rT) 


g(r, x) =H(r, x) —G(r, x), 


ri =kilkio, (A-6) 


so that the definitions differ in the regions r 2 a. The 
relations in (A—5) may be inverted to give, e.g.: 


=G(1+7(riG)), (A-7a) 





] 
a 
T(x) = i 


] | 1 | 
: 4°93 x| 


In terms of the Laplacian 


(A-7b) 


H 


and in the region r>a, (i=2) Laplace’s equation 
holds exactly. This is the critical region for heat loss. 


292-251 O-68—2 


We shall suppose the form of the rf-wave is given, 
but its amplitude 79 is determined (for given .7;, 72) 
to balance the bridge after substitution 


F p(x, t) =F roz(x) cos wt. (A-9a) 


Then taking time average 


([%2+ Froz(x) cos wt]?) = 4 +5 Foz? (x) 


Let ( ) denote the average on (0, /) and define 


ais) 


21 (x ems 2(x) 


zi(x) =0, (A-9b) 


so that 


I4+5 Fhoo(x) =F i +5 Fheo2” (x) (1 +21 (x). — (A-9e) 


In the long panna limit, z2(x) is a constant, 


z1(x)= 


0, so that if 73+5 5 F hoz (x) x) 


=JF2, the sources 


S; and Sy are identical, so are the fields 7;, 7, and 


there is no substitution error. We shall assume then 
that 


5 F oF (x) 2— 924 22, (A-9d) 


and thus that 


Su= [$+ (#— (+P) (A1+21(x))] (pot piTu), 


(A-9e) 


where /?/(.77—.72) and z;(x) are assumed small com- 
pared with unity. We must expect g(r, x) to be 
small also whenI?/(.4?—%3) and 21(x) are small, and 
we linearize (A-8b) accordingly in g: 


kino Lg=—8[47Bg+ BCg+ AC}, 


where 
A=potpiG(1+7(7,6)) 
B=pi(1+7(nG+nGr' (71,C)) 
C=R+ (2 —JI2)z1(x) + 0 (I?z,(x)). 





ia fo) feet | 


i ((%ot+F%r(x, t))?) [Pot piH(1+7(nH))) |, 





Now C is assumed small by virtue of /?/($2—.43) and 
z1(x) being small, so that 


BCg « $2Bg. 
Also since 
Pigs < pot PG, 
BCg < AC. 


Thus we can neglect the second source term. We 


shall also neglect O(/%z,;(x)) in C. Then we have 
kin Lg =— Si {e Fp, (1 +7(r4G) 
+ rmGr' (mG) ) . [pot piGi +7(mG))]- 
[12+ (75 — F921 (x) ]} 
We shall deal only with G(A-8) and g(A-11) to 


which must be appended boundary conditions: 


aC _ag_ 


r=0 
or or 


0,(0<x 


r=b: G=g=0, (0<x<l); 


x=0, 1: G=g=0, (0<r<b); (A-12) 


while the interface continuity conditions become 


aG |@* 

> ki —— =0, 
ar a- 
Kio 2£|" =0, 


io 
or a” 





for continuity of heat flow, and 
r=a:|G(1+7(rG))\¢7=0, 


lg(1+7(riG) + riGr' (riG)) jar = 0, 
(A-14) 


where we have again used |g|<|G| in the second 
condition in (A—14). 

Now even in the region of the wire where the tem- 
peratures are highest, 7(7r:C) ~ 0.1, so that the non- 
linearity is small. Furthermore, because of the high 
conductivity of the wire compared to the air, the 
radial temperature gradient in the wire is negligible, 
and 


G(r, x)= G(0, x) (A-15) 
inside the wire. 

The equations are still nonlinear in a way prevent- 
ing the easy use of linear methods of solution. Extend- 
ing the assumption that t(riG) is only a small correc- 


(A-11) 


tion at worst, we assume that the average axial value 
of G may be used to evaluate 7(riG): 


r(riG) = 7(rG) (A-16) 


Then we shall use the differential equations 
inca SF ¥pi (1+ 71) 8G =— F Fpodis 


king t+ F tpi (1+ 72) 8g 
=— [P+ (4?—43)z(x)] 
: [pot priG(x, 0)(1 +71) 8: 


where 7; = 7(r:G) , T= 1 +Gr' (r:G). 
On r= a, we have the interface condition: 


G(a-, x) =G(at, x)o; 
r=a 


(a~, x) =g(at, x)or, 


_1+7(nG*) 
sill 1+7; 


_14+7(nG*) + rG*r' (RC*) 
7” 1+ 72 
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where G* is the number which gives continuity of 
temperature at r= a* in the mean by (A-7a): 


G*(1+7(mG*))=G(1+7(nG)).  (A-19) 

The system of equations we solve then, (A-17), 
subject to conditions (A-18, A-13, A—12), are linear 
when G is given. From the solution, a new G is com- 
puted, and the solution iterated to make G stationary. 
With a reasonable guess, say only 25 percent off, five 
iterations brought G to within 1 percent. 

No attempt has been made to bound the uncertainty 
of the substitution error due to the assumptions made. 


Possibly a better choice for the constant G than that 
given in (A—16) might be made. If a 5 percent increase 
in G is used to compute the (7;, oj), the substitution 
error changes by 0.06 percent in a typical case. The 
approximation is clearly better for large //b, for which 
the temperature profile is flat for;most of the wire 
length. . 

The constant / in the source of (A—17) for g is deter- 
mined by the requirement that the total bolometer 
element resistance be unchanged when the d-c current 
YF, is replaced by the current ((.42+.%x(x, t))). 
Equating these resistances: 


[, atm + piG(o, x1)(1+ 7(riG(O, x)))] 
0 


= I, dx[pot+ piH(0, x) + r(r,HO, x)))| (A-20) 
0 





Expanding H to lst order in g, and again assuming 
G(0, x) is flat enough to be replaced by its average in 
7(riG), we find 


l 
I. dxg(0, x) =0; (A-21) 


for the balanced bridge this relation determines 
I in (A-17). 

When no rf power is introduced, the d-c power 
developed in the bolometer element is 


l 
Woc = tra? J? I. dx [po + pi G (0, x) +7(r1G))] (A-22) 
¢ 


When the d-c current is reduced to .%2 and rf current 
is added to rebalance the bridge, the total power 
developed is 


l 
Wor = ra? | 


0 


dx [43+ Fj z?(x)] 


* [po+ pil 0, x)(1+7(r1H))) 


= Woc2t+ We 


(A-23) 


where by virtue of (A—20); 


Woc2= 0 Woe, C= 2/5. (A-24) 


The substituted power is defined as 
Wsun=Woa — Woe = (1-2? W ver (A-25) 


The change in power developed AW (with A, B, C 
as before (A—10), and using the same approximations 
and the bridge balance relation (A—21)) is found to be: 


AW = Wrox — Woe = Wer — Wr 


l 
= nat | dx(AC + £?Bg) 
0 


I 
= na | dx[pot piG(0, x)(1+71)] 
0 


[2 +92 —$22(x)].  (A-26) 


Then the substitution error is 


_ AW 
— Wsup_ __/ sus . 

i+ AW 
W sus 


a (A-27) 


Let x(x, r; €, n) be Green’s functions for the sys- 
tems (A-17) with the superscript (j= 1, 2) referring to 
(71, 01) for G - (t2, G2) for g. The Green’s function is 
symmetric in (x, €). Let: 


GI()=— (14+ 7)p 53 [ dé I “daenitte, Os 8/90. 
0 
(A-28) 


We define the numbers 


Ai, =[ dx(1+ GG 1a) |¢ 


G(x | 


By. =| dx(1 + G(x))z1(x) Lan i (A-29) 


The bridge balance condition (A—21) leads to 


P Bs 
~~ | A-30 
R-H As lle, 


In terms of the (A;, Bj), the power terms can be writ- 
ten as 


W sup = (1 — p*)7a*poF4A 1, (A-31) 


AW = ra*po.F?(1 — p?) rae A+ Bi}, (A-32) 


i? 
— $2 
so that the substitution error is 


A,B.—B,A2 E _ A,B vil 


ee = A,A> 


(A-33) 


The axial temperature distributions themselves are 
found from 


T(x) = *oW sin <i 


(A-—34a) 
AT (x) = Ty(x)— T;(x) 
(A-34b) 


where >* is a sum over odd n only, the coefficients 
gi) are given in Appendix B, and 


I 
yf2)= ni?) oA Asn + Baa| 


(2) 
—— (A-35) 


: (BsnA2 — AsnB2), 


a a) nee} |: 
beg : dx(1 + G™(x)) sin ] vr 


Our primary concern is with a sinusoidal current 


distribution 
= (1x 
z(x)= cos ad (J *) +9| 


from which we compute z;(x) by (A—9a). 


(A-36) 


(A-37) 
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Then 


2*(x)= cos? ¢y cos? ki (3-4) 


+ sin? ¢ sin? k/ (3 ~5)+ A(x), 
\ 


] (A-38) 


where 


A(x) = —* sin 2g sin 2kl (5-4) 


— 


does not contribute since it is anti-symmetric on (0, /) 
while G(x) are symmetric. 


We find, averaging and noting . A(x)=0. 


ain MY con K! 

a ; Phat: 9 Cos 2 

z-(x)= cos* i ie ea 
Pn oe kl 


me kl kl 
ale = one = 
l " 2 as ay (A-39) 
‘(Ht ——_———— 


2 kl 


We obtain z,(x) from (A-38, A-39) using (A—9b). 

In Appendix B, details of the calculation of the 
G(x) by Fourier series, description of a program used 
to compute the (4;, Bj) and the temperature distribu- 
tions are presented. 


Appendix B. Green's Function 


We are to determine the Green’s function 


Xx, r; €, n) satisfying the system 

kin. SX)? i Stpill sya! o(x =F) -Or— nN): 6i1 
c (B-1) 

subject to the conditions 


he O<xl; 


Nx, r; €, n)=90, 
Xx, 13 E, M=O =(0,1),0<r<b, 


(B-2) 


ay) 
— (x, 0; €, n)=0, 


j 2 j) 
kio (x = k2o 2°) 
Or /r=a- Or /y=at, 


(XK) )p—a7 sail o(X))r=a*. 


A Fourier sine series for the solution is obtained by 
taking the finite Fourier transform. 


(B-3) 
(B-4) 


(B-5) 


l 
x(a, rs €, n={ dxx(x, r; €, ) sin ax, 
0 


Ci ie 
x(x, rz &, =F SY x"(an, 73 €, 0), 


n=0 


nT 
An = 


l 


which satisfies the conditions (B—-2) at x= (0, 1). We 
obtain the equations (dropping the superscript (j)): 


ly! 
kyo : i a + (42p1(1 + 7) — kioa”)x! 


=sin a& <2 


(r< a): (B-7a) 


" 
ko FX bypa?y? =0, (r> a). (B-7b) 


For r>a, the solution satisfying (A—2) is: 


xa, =D, | fotar) — fated kean)| (B-8) 
ola 


For 0 <r<_>, the solution satisfying (A—3) is: 


ya, r)= Dato (F c). (B-9) 


= p’—(al)’, 
B =F? pil(1+7)/ki. 


For n < r<a, the general solution is: 


x" (a, r)=DsJo (F«)+ DsYo (; e): (B-11) 


The four constants D; are readily found from con- 
ditions (B—4, B-5) at r=a, and from jump conditions 
at r=7: 


x"(a, n-)= x(a, n*)s 


We obtain for the axial value 


x0, r; €, 0) = ik > Fae (Ben)+¥ (2 .)| 


. mne . TNX 
‘sin —~ sin 


; sn = 18) 


1 (fe) orem (4) 


aekio l 
; * anlkoo (5 ) 
" & cx) ie 


« — 1o(a0n)Ko(ban) — 
— 1 \(aan)Ko( ban) =, 


R,=- 





To( ban)Ko(aan) 
To(ban)Ki(aen) 





(B-14) 





All that is required in the text is the integral of this 
on 7. We obtain: 


a ‘ oe l 
dinyx Hx, 0; €, n) =: x 
j - - é ”) 2k 1085 n=0 


— NX 
no) sin mae sin 2. 


l 
(B-15) 
op = 46L 


Tne? €ink 
J Jo (5 oy) + Ee Ji (5 ci) 





—\lI- 


(B-16) 


Note that the integral (B—15) is symmetric in (x, &). 
From (A-15), we obtain 


t a 
GY (x) =—pr42(14+7;) [ dé I dnnx'd)(x, 0; E, ) 


l a 
=—pA1+n) | ae | dnnx (E, 0; x, N) 
0 0 


(B-17) 


ES si cs 
* bi) sin>— 


n=1 


The Bessel functions occurring in ¢¥) are computed 
from polynomial approximations [9]. 
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Analysis of n-Degree Elliptical Elastic Rings of Non- 
uniform Cross Section 


R. A. Mitchell 
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A curved bar structural analysis is developed for an n-degree elliptical elastic ring with sinusoid- 
ally varying cross-section dimensions, where the n-degree ellipse is defined by 


(=) +(5) =1, n>20. 
a 6 


A sequential grid constrained optimization method is used to search for a minimum weight design in 
six dimensional shape parameter space. Numerical results indicate that rings of this shape can be 
designed to have considerably less weight and greater flexibility than comparable circular rings with 
uniform cross-section dimensions. 


Key Words: Curved bar analysis, elasticity, force transducer, n-degree elliptical ring, nonlinear 
programming, optimization, proving ring, structural design, theory. 


1. Introduction 


Elastic ring force transducers (such as proving rings, load rings, ring dynamometers, etc.), 
for a wide range of load carrying capacities have most frequently been made circular in shape with 
uniform rectangular cross section. This shape has imposed limitations on the usefulness of this 
type device, particularly for higher capacities (i.e., greater than 100,000 lbf), due to overall size, 
gross weight, limited deflection range, and size effect during heat treatment. For the same force 
resisting capacity, rings of more complex shapes can be designed to have greater deflection ranges, 
less weight, and smaller maximum dimensions. This paper describes methods of analysis of a 
broad geometric class of complex shaped rings which result in improved deflection and weight 
characteristics. 

Structural analysis and design optimization methods are developed for an n-degree elliptical 
elastic ring with sinusoidally varying cross-section dimensions. The ring is analyzed as a thick 
curved bar, assuming small deflections, plane strain, Winkler-Bach bending stress distribution, 
uniform normal stress distribution, and parabolic shear stress distribution (Fy. A sequential grid 
constrained optimization method is used to search for the minimum weight ring design that will 
satisfy the specific force capacity, deflection, and dimensional requirements of a force transducer. 


' Figures in brackets indicate the literature references on page 150. 





List of Symbols 


x, y rectangular coordinates of ellipse, 

r,@ polar coordinates of ellipse, 

degree of ellipse, 

semimajor and semiminor axes of ellipse, 

ring cross-section width, 

ring cross-section thickness normal to tangent of ellipse, 
average of maximum and minimum ring width, 
average of maximum and minimum ring thickness, 
dimensionless parameter of ring width variation, 
dimensionless parameter of ring thickness variation, 
arc length variable measured along ellipse, 
counterclockwise angle from radius vector to tangent line, 
total load applied to ring, 

= L/2, load applied to one quadrant of ring, 
moment, 

normal force, 

shearing fofce, 

normal stress, 

shearing stress, 

strain, 

curvature of ellipse, 

relative rotation of differential segment ends due to N, 
curved bar rotation, 

curved bar deflection, 

total ring deflection, 

strain energy of curved bar, 

strain energy per unit volume, 

Poisson’s ratio, 

modulus of elasticity, 

integral defined by eq (19), 

design maximum stress, 

scale factor, 

ring weight, 

material weight per unit volume, 

constraint function, 

constraint bound, 

variable of the optimization problem, 

objective function. 
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2. Structural Analysis 


An elastic ring force transducer is loaded by two forces acting in opposite directions along 
a diameter, as shown in figure 1. Corresponding ring deflections are measured and related to 
applied forces by a calibration factor.2 For a ring of the geometric class represented in figure 
1, the locus of centroids of ring cross sections is an n-degree ellipse defined by the equation 


cy att .. 
(*) +(%) =1, »a2.0, (1) 


? Wilson, Tate, and Borkowski [2] have described the uniform cross section, circular force transducer, its calibration, and performance under various conditions 
ot use. 
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FIGURE 1. 1,000,000 bf capacity n-degree elliptical ring with sinusoidally 
varying thickness. 


or in polar coordinate form 
r=ab(b" cos" 6+a" sin" @)-"/". (2) 


Several n-degree ellipses are plotted in figure 2. Ring rectangular cross-section width h and thick- 
ness ¢t are given by 


h=c(1—e cos 26), (3) 


and 
t=d(1—fcos 28), (4) 


in which c and d are positive parameters and e and f are parameters of absolute value less than 
unity. Using these equations, the shape of a ring within this geometric class can be specified by 
the seven parameters a, b, c, d, e, f, and n. If a=b, n=2, and e=f=0, the ring shape is circular 
with uniform rectangular cross section. 

One quadrant of the loaded ring is shown in figure 3. Due to symmetry there is no rotation at 
6=0 or at 9@=7/2, no shear force V at 9@=0, and no norinal force N at @= 7/2. Expressions for the 
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FIGURE 2. n-degree elliptical curves. Be OF 
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FicuRE 3. Loads acting on one quadrant of a ring. 





resultant moment M and forces N and V acting on a typical normal cross section are derived by 
statics using formulas from calculus [3] as follows: 


d6_ sin p_ 1 
ds r dr\2. (9) 
V+ (a) 


dr dé 





M=M,)2—Prcos 0=PM', 


ee d0\ _ py 
N=P sin (6+) =P (sino $+ reos 65°) =PN’, 


V=—P cos (0+ w=P(r sin on 6 o)=Pv". 
ds ds 


Differentiation of eq (2) gives 


dr__(b" cos"~” @ sin @—a" sin"~ @ cos 8) 
do” (b" cos" 0+ a" sin" 6) ‘ (10) 





The assumption that ring cross sections remain plane leads to the Winkler-Bach bending stress 
formula, which has been written by Seely and Smith [1] in the following form, valid for a straight 
or curved bar: 


fe me". } (11) 
1+ me a ee , 
a+ | Taye 


The variable y is the distance from the cross-section centroid to a point in the cross section, positive 
in the outward direction. The curvature k, the inverse of the radius of curvature, is given by [3] 


2 2 
r?+2 (3) _~,£! 


dn2]32 (12) 
24(4@r 
f +() 


Differentiation of eq (10) gives the required second derivative 





Ou 


k= 





a : 
—_— (n+1)ab[a" sin'"-” @ cos 0—b" cos"-” @ sin 0]2[a" sin" 6+ b" cos" 9]-2+1 


+a'"*” b[sin" 6@— (n—1) sin"-» 6 cos? 6][a” sin" 0+ b" cos" 9]-+1/ 


+ ab'"*" [cos" @— (n—1) cos» 6 sin? 6][a" sin” 0+ b" cos" O]-U+U”, ) (13) 


It is assumed that stress due to the normal force N is uniformly distributed over’ the cross 
section, that is 


ae. ? 14 
On ht’ | ( ) 
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It is assumed that shear stress is distributed parabolically over the rectangular cross section 


according to the formula 
V (3 y? 


Castigliano’s strain energy theorem [1] is used to determine the ring bending moment and 
deflection at 9@= 77/2. The theorem stated symbolically as applied to this problem is 


aU 
On = SS, 
a/2 dM aie 
au 
bz/2 ai aP’ 
in which U= total strain energy of the loaded bar, 
0,/2= bar rotation at 7/2, 
and 52/2 = bar deflection at 7/2 in direction of P. 
The strain energy per unit volume due to bending stress is 


| _U-v*)o%, 
um= 2 Cuteu= oF ’ 


in which €y = unit plane strain, 
v= Poisson’s ratio, 
and E= modulus of elasticity. 


Substituting eq (11) into (18) and integrating over a differential segment gives the corresponding 
strain energy of the segment 





— 2 " 2 es 2 
Bin Eis. Tes Pu ee Sats 
Eh F cee J - Eh 
t t Ql + ky) q ; 


The strain energy per anit volume due to the normal force N is 


aie ev=5 (Fr) (*)- N? 
tN 2 NEN 9 Nit} \AtE)— 2ER?t?” 


and the strain energy of a differential segment is 


dUy ds. 


~ QEht 


The strain energy per unit volume due to the shear force V is 





oe 


(l+v)2 (1+)? (3 y*\? 
— =: oe G- 7) 


and the strain energy of a differential segment is 


6 /l+v\ Vv? 
dUy =< )5 Ss: 
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The normal force N acting on a differential segment causes rotation of one end cross section 
relative to the other, due to the difference in arc length from inside radius to outside radius. The 
mean change in segment arc length is 


N | 
Eht ds = (;) B. 
in which B= rotation due to N. 


The bending moment M acting during this rotation contributes strain energy to the segment 
equal to 


kMNds 


dU yy = ~Eht 


Adding eqs (19), (21), (23), and (25) gives the total strain energy of a differential segment 


w= N? oa Pe ,, Mg MNk) 
c =Elaat3 + vp) tg) + | 


Applying eq (16) and the boundary condition of no rotation at @=7/2 gives 


aU Lif" yN aN 1 Vv wv eM aM kM aN kN aM 
OMre E Jo aac = O40 iG Sint) Site Be Mest Te Mew © 


"2 F(1—p2)gM kN] . 
|, +I ds=0. (27) 


Substitution of eqs (7) and (8) into (27) gives 





[ (l—v*)gr cos 0 i k sin 6 
| i a) ar—| 
jes — ds 
h 





Applying eq (17) to (26) gives 


bzj2= 


aU | [ E aN 141+vW av ,. &M aM kM AN kN > 1 
—_ ds 


‘PE ht aP 5h aP*"2-”“) 7 opt he apt ht oP 
P m/2 | 


12 
a : lA B+ = (1+v)V'?—(1—v*)gtr cos OM’ + kM'N' —kr cos an’ | ds. (29) 


The tangential stress on the inner surface is provisionally assumed to be the limiting stress 
that determines the force resisting capacity of the ring. The possibility of capacity being limited 
by some other component of stress or combined stress should be checked in any final design 
analysis. Substituting —¢/2 for y outside the integral in eq (11) and adding to eq (14) gives the 
tangential stress on the inner surface of the ring 


elie M't a4 
OTr= —_——_—— - Pr — ie 
ht y* ht 
h(2—k / 
( lam dy 
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For a circular ring of uniform cross section the maximum value of o7 occurs at @= 77/2. But for the 
more general geometry considered here the maximum may occur at 6=0 or at some point in 
between. 


The structural analysis thus far has been developed for one quadrant of the ring. For the total 
ring, figure 1, the applied force is 


L=2P, (31) 
and the deflection at the point of load application is 
A = 28z/2- 


3. Design Analysis 


The design analysis attempts to find the dimensions of the minimum weight ring that will 
satisfy the arbitrary force capacity, deflection, and dimensional requirements of a particular force 
transducer. The deflection range of the ring must be compatible with the deflection sensing device, 
and the inside dimensions of the ring must be sufficient to accommodate the sensor. Outside di- 
mensions of the ring must not interfere with its intended use. The thickness of the ring must permit 


the heat treatment required for the ring material. All trial solutions are scaled up or down to 
coincide with the force capacity requirement. 


3.1. Scale Factors 


The procedure adopted here for meeting the force capacity requirement is to scale all length 
dimensions in the r-@ plane by a factor A that will make the maximum stress (the maximum value 
of or for the entire ring) equal a prescribed maximum design stress S. The cross-section width h 


is not scaled. This gives the specialized equation for the maximum tangential stress on the inner 
surface of the scaled ring 





) M'k _ M't N' 


causiteiniates +— 
ht y? ht 
h(2—kt | d 
( ) (1 + ky) dl MAX. 


LaF 
2d” 


in which S= design maximum stress, 
L,= design capacity, 
and F=maximum value of bracketed function. 


This gives the scale factor 


LF 
A= "9S 


Scaled ring dimensions, \ subscripted, are 


and 





Scaled forces and moment are 


Scaled ring deflection is 


Ay =)A. 


Scaled ring weight is 
1/2 
n= 4y)? | htds = ?>W 
0 
in which y= material weight per unit volume. 


3.2. Thin Circular Ring 


By making the appropriate simplifying assumptions, the thick ring equations developed above 
can be reduced to relatively simple thin circular ring equations. The thin circular ring equations, 
although they may be in error by several percent, can be used to determine a reasonable range of 
ring proportions for beginning an optimum search procedure. The equations also give a useful 
indication of the comparative efficiencies of different ring materials. 

A thin ring is here defined as one in which the ratio of thickness to radius of curvature is so 
small that, with negligible error: (1) the stress distribution due to bending can be assumed linear 
over a cross section, and (2) deflection due to shear V and normal force N can be neglected. By 
assuming zero effective curvature and no deflection due to V or N, the following thin ring equations 


for moment, maximum bending stress, and deflection at the point of load application can be 
obtained from eqs (28) through (32): 


orm he 


- _ 24\ d— y)Lr 
A= (30 ) Fhe 


These equations, solved simultaneously, give 








a i 72EL?A 
(4 — 8rr2)(1 — v2) 3h?” 


tro ht? 


—— 
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The weight of a thin circular ring is 


24EyLA 





W = 2arhty=7 


@—8\1— Po ) 


Note that the weight, for a particular LA, is a function of the material properties E, v, a7, and y. 
Equation (48) suggests the importance of maximum permissible design stress in a comparison of 
different ring materials where weight is a significant factor. 


3.3. Optimization Method 


The optimization problem considered here is of the general nonlinear programming class. 


The goal of this type problem is to determine the set of values of the variables x;,. . ., xp which 
satisfies the constraint conditions 


BhRas- » + xt <, =, =}b, i=1,... (49) 


and minimizes the objective function 


z=f(x1,. . .» Xp). (50) 


In the present case the p variables are the six independent ring shape parameters a, c, d, e, 
f, and n. The structural analysis is first done for a ring of unit value of mean radius, that is 


a+b _ 
D =]. (51) 


Thus the value of the parameter a determines the value of b. The thickness parameter d is, for the 
unit size ring, the ratio of the mean of the extreme values of thickness to mean radius. After the 
structural equations are solved for the unit mean radius ring (using Simpson numerical integration), 
the results are scaled to full size using the scale factor A. The width parameter c is not scaled 


during the entire analysis, and the shape parameters e, f, and n are nondimensional and are not 
scaled. 


The m constraint functions g; and the corresponding constraint bounds 6; are as follows: 
(1) Ring deflection 
&=A\2 bi, 
(2) Maximum thickness of ring cross section 
8 =dd(1+|f]) <b, 
(3) Outside width of ring 
&=d(2a+d—df) < bs, 
(4) Outside height of ring 


g4=A(4—2a+d+df) < by, 
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(5) Maximum width of ring cross section 
g5=c(1+|e|) <b 
(6) Inside clear height of ring 


gs=(4—2a—d—df) > be. | (57) 


The constraint bounds b; are the limits on ring characteristics imposed by use, fabrication, or 
material requirements. 

The objective function z to be minimized is the ring weight Wy, eq (42). 

A sequential grid method was used with reasonable success to search for minimum weight 
solutions. The method requires several computer runs of moderate duration. Between computer 
runs changes are made in the input data based on the experience of previous runs. The “Elastic 
Ring Optimization Program” listed in the appendix was used for this search procedure. 

For the first computer run, one or more values of each of the six variables x; are chosen as 
input data. If two or more values of a variable are chosen, they are widely spaced, but within a 
reasonable range for that variable. By nested do loops the computer generates all combinations of 
the variables. Each combination is a solution point in the six dimensional shape parameter space. 
The ring weight and the six constraint function values are computedand printed for each solution 
point. The best solution point for a computer run is the one with the least weight that also falls 
within all six constraint bounds. 

Subsequent computer runs are made using as input data variable values distributed about the 
best known solution point. The range of values of any variable should include the value of that 
variable for the best known solution point. Any results of previous computer runs may be useful 
as a guide for selecting input data. The process is repeated until a solution point with satisfactorily 
small weight is found, or until results indicate that a satisfactory approximation of a minimum 
weight solution has been found. 

Computer time required for each solution point is approximately one-half second on a UNIVAC 
1108 computer. Obviously there are practical limits to the number of variable combinations one 
would want to generate in one computer run. 

A second computer program (the “Detailed Ring Analysis Program” listed in the appendix) 
can be used to compute scaled dimensions, forces, moments, and stresses for a particular set of 
input shape parameters. The program prints ry, t,, h, k, Ny, Vi, My, or, and the value of the integral 
in eq (11) for each incremental value of 0. These results can be used with eqs (11), (14), and (15) to 
determine whether ring capacity is limited by a stress condition other than the maximum value 
of or. 

The following numerical example demonstrates the optimization method. 


3.4. Numerical Example 


The problem is to determine, for a ring of 1,000,000 Ibf capacity and for the material constants 
and constraint bounds listed below: 
(1) The appreximate weight of a circular ring of uniform rectangular cross section. 
(2) The dimensions of an approximately minimum weight n-degree elliptical ring of constant 
cross section width and sinusoidally varying thickness. 
(3) The dimensions of an approximately minimum weight n- "degree elliptical ring of sinusoi- 
dally varying cross section width and thickness. 
The material constants are: S= 150,000 lb/in?, E = 30,000,000 Ib/in?, v»=0.3, and y=0,29 Ib/in?. 
The constraint bounds are: 6; = 0.25 in, b2= 4.5 in, b3 = 40.0 in, bs = 40.0 in, 6b; = 14.0 in, bg = 9.0 ins 


a. Solution 1 


Thin circular ring equations were used to compute a rough approximation of ring weight and 
dimensions. Substitution of L, E, y, v, S for or, and 6, for A into eq (48) gives an approximate ring 
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weight of 1365 lb. Substitution of the same values, along with 6; for A, into eqs (46) and (47) give 
t=3.66 in and r= 14.73 in. The ratio of thickness to radius is 0.248. The more accurate thick 
curved bar equations were used for all subsequent computations. 

Computer run No. 1 of the sequential grid search included nine uniform cross-section circular 
rings with widths of 10, 12, and 14 in and thickness to radius ratios of 0.20, 0.25, and 0.30. For 
these nine cases the ratios of weight to deflection differed by less than two percent, a result that 
might have been anticipated from eq (48). The two cases with characteristics nearest the constraint 
bounds had deflections of 0.234 and 0.277 in and weights of 1327 and 1587 lb respectively. Linear 


interpolation between these two cases gives an estimated weight of 1424 lb for a ring that would 
deflect 0.250 in. 


b. Solution 2 


A sequence of five computer runs were used to seek the dimensions of a minimum weight, 
constant width ring. The input shape parameter data and the resulting weight for the best point of 
each run are given in table 1. The best point parameter values for each run are underlined. The 
weight of 553 lb for the fifth run is about 39 percent of the estimated weight of 1424 lb for a uniform 
cross-section circular ring. The 553 lb ring was accepted as a satisfactory approximation of the 
minimum weight constant width ring. The ring shape is shown in figure 1. Critical ring dimensions 
and a deflection of 0.257 in all fall within the constraint bounds. 

Plots of or versus @ for the 553 lb ring and for a 1587 lb uniform circular ring (c= 14.0 in, 
d=0.25 in, A=0.277 in in computer run No. 1) are shown in figure 4. The high stress level over a 
greater portion of the 553 lb ring is an indication of more efficient use of material in this more 
complex ring shape. 


¢. Solution 3 


Computer runs 6 and 7 included many cases with variation of cross-section width. But for both 
runs the best points were constant width cases. The seventh run included a variable width case 
with parameter values c= 13.6 in, e= 0.025, and f= 0.45, and a resulting weight of 539.18 lb which 


TABLE 1. Sequential grid search data for numerical example 
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is only slightly greater than the best point weight. Additional searching might result in a slight 
weight reduction due to width variation. But the weight savings would probably not be enough 
to justify the additional complexity of shape. 


4. Summary 


A structural analysis method has been developed for an n-degree elliptical elastic ring with 
sinusoidally varying cross-section dimensions. A sequential grid constrained optimization method 
has been used to search for the minimum weight ring of this geometric class that would be suitable 
as the load supporting element of a force transducer. 

Limited numerical results were obtained for one-million pound force capacity n-degree ellip- 
tical rings with nonuniform cross-section dimensions. The results indicate that such a ring with 
constant cross-section width and sinusoidally varying thickness can be designed to weigh as little 
as 39 percent of the weight of a comparable uniform cross section, circular ring. Outside dimensions 
of the more complex ring are also significantly less than those of the uniform, circular ring. The 
results for an n-degree elliptical ring with sinusoidally varying cross-section width and thickness 
suggest that the weight savings due to variation of width may not be great enough to justify that 
additional complexity of shape. 

The Elastic Ring Optimization Program (given in appendix) can be used directly for design 
optimization by the sequential grid method described in this paper; or the program could be used 
in combination with some other optimization procedure. In either case, the Detailed Ring Analysis 
Program (see appendix) or another comprehensive stress analysis should be used for a detailed 
examination of stresses for the entire ring. 


The analysis was coded in FORTRAN V for the UNIVAC 1108 computer by R. M. Slesser 


of the National Bureau of Standards Center for Computer Sciences and Technology. 
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6. Appendix 


The two FORTRAN V computer programs used in the elastic ring design analysis are listed 
below. The Elastic Ring Optimization Program is used to search for the optimum set of shape 
parameters a, c, d, e, f, and n. The Detailed Ring Analysis Program is used to compute scaled 
dimensions, forces, moments, and stresses for a particular set of input shape parameters. Table 
A] lists the computer symbols used to represent mathematical symbols appearing in the equations 
of the text. Sample input data is included at the end of each program. 





TABLE Al. Corresponding symbols in text, computer programs, and 


computer output 








Programs 


Output 
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rT SSR SF BIT SHV we MHHSSEAVS 


EM ERPSBRY 
> 


MW 

2M 

Ar 

At 

OT 

6 for max. or 
WyIAy 

Integral in eq (11) 
ds 


dr 
do 
&r 
d? 
dé 


ds 
dr 





R(IJ) 
THETA(IJ) 
RHO(IJ) 
GI) 
PM(IJ) 
PN(IJ) 
PV(IJ) 
DMAX 
XLAM 
WLAM 
DELAM 
TCON4 
TCONS 
TCON6 
TCON7 
TCONS8 


ud 
XLN(I) 
XLV(I) 
XLM(I) 
RLAM 
TLAM 
SIGTL(D) 
JMAX 
RATIO 
YINT(IJ) 
DS 


DRDT(IJ) 
SDRDT(IJ) 
DTDS 


DRDS 





AOZAMr ANAS 


Tl 


LAMBDA 
WEIGHT 
Gl 
G2 
G3 
G4 
G5 
G6 


LAMN 
LAMV 
LAM**2*M 
LAMR 
LAMT 
SIGMA 

™ 
WT/DEFL 
INT 








Elastic Ring Optimization Program 


i RUN MITCHE? 16650202060 
wTl FOR MAINI*MAINI 
COMMON /INTEG/F (100) eAREA 
UViMENSION YINT(91) 
O{IMENSION PV(91) 
DIMENSION C1(20) C3120) 0C5(20) 6 C7(10) »C8(10) eXN(40) 0 B91) 0 T(91) > 
1 (91) eDRDT(91) eSDRDT(91) eRHO(91) 0G (91) 6PM(91) ePN(91) 
2 €2(20) e THETA(91) 
one READ(Se 10) NleN3SeNSeN7eNBeNNe IGO 
io FURMAT (2024) 
C IGOeNEed MEANS ANOTHER SET OF DATA FOLLOWS THIS SET. 
; ReAv(S 0©20)(C1(1) eI=10eN1) 
READ(S 0©20)(C3(1) oe I=19eN3) 
Re AD(5 0©20)(C5(1) eI=ieNS) 
REAU(S #20) (C7(1) eIS=ieN7) 
ReAU(S 0©20)(C8(1) eo L=19N8) 
ReAv(S 020) (XN(1) eI=1eNN) 
20 FORMAT (10E8.0) 
ReAD(S #©30) XNHUeXLLAMeSeE* GAMMA 
30 FURMAT (5615-0) 
C INSERT INITIAL wRITE CARDS HERE. 
C PRINT OUT ALL INPUT DATA AT START OF RUNe PROGRAM le 
WRITE(601000) NieN3eNSeN7eN&eNNeIGO 
1000 FORMAT (ISHLINTEGER PARAMETERS //(8I15)) 
wR1TE (601010) (Ci(I) oe I=1eN1) 
WRITE (601020) (C3(1) ep IT=1eN3) 
wkLTE (601030) (C5(1I) ep I=1+N5) 
WKITE (601040) (C7(1I) e IL=10N7) 
#RLTE (601050) (C8(1) eI=1eN8) 
WRITE (601060) (XN(I) ep I=L e NN) 
FORMAT ( 3HOC1 //(8E15¢5)) 
FORMAT (3HO0C3 //(8E15¢5)) 
FURMAT (3HOCS //(8E15e5)) 
FORMAT (3HOC7 //(8E15¢5)) 
FORMAT (3HOC8 //(8E15-5)) 
FORMAT (2HON 4/(8E15¢5)) 
wRITE (601070) 
FURMAT(AHO 11X*°3H NUP15Xe 7H LAMBDA? 18Xe1HSe19XelLHEe 14XeSHGAMMA ) 
wK1TE( 601080) XNHUeXLLAMeSe Ee GAMMA 
FORMAT (1HO*°5E20-8) 
JCOUNT=0 
PI =3-141592654 
Du 1 I=leN1 
Cel(I)=2.-0-Ci(I) 
DU 1 J=leNS 
vv 1 K=leN5 
DU i L=leN7 
Du 1 M=1leN8 
OU 1 N=1eNN 
A=Ci(I) 
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C vi 


121 
125 


122 
1235 


C RHO 


Elastic Ring Optimization Program (continued) 


D= CatlI) 
GOcS FROM 0 TO 90. 

Du 2 IJ =1°91 

JIl= [J-1 

THETA(IJ) =FLOAT( JI) *PI/180-0 

bIIUI= C3(U)#(1.-C7(L) #COS(2e*THETA(IU))) 

TILUJ= C5(K)4#(1-¢"-C8(M) *COS(2¢*THETA(IY) )) 

TERM = A **XN(N) *SIN(THETACIJ) ) #*®XN(N) 4D **XN(N)* 

COS(THETACIJ) ) **®XN(N) 

R(IUJ= A*D /TERM*¥* (1 6/XN(N) ) 

TPl= O **XN(N) *COS(THETACIJ) ) ** (XN(N) 120) *SIN(THETA(TIJ) ) 
TP2 = A **XN(N) *SIN(THETACIJ) ) *®*®(XN(N) =1.0) *COS(THETAC(IJ) ) 
DROT(IJI=S RI IJ) *«(TP1-TP2) /TERM 

IF (XNON) eEQ@e2e0eANDe THETACIU) eEQe000) GO TO 121 

TPS = A ** XN(N)* 
1 (SIN(THETAC IJ) ) #*#XN(N) = (XN(N) =100) *SIN(THETACIJ) ) ¥* 
2 (AN(NI“2e) * COS(THETA(IJ) )*#2 ) 

IF (XN(N) eEQe2e0eANDeIUcEQ-91) GO TO 122 

60 TO 125 

TPS = A*¥*2*(-1.0) 

CONTINUE 

TR4= D ** XN(N)* 

1 (CUS(THETACISJ) ) ®*®XNON) —(XN(N) 120) *COS(THETACIJ) ) * 
2 (XN(NI©20)* SIN(ITHETACIU) ) **2 ) 

GO TO 123 

TPY = D¥*2 * (-1.0) 

CONTINUE 

SOROT(IUI= A*D *(XN(N) 41.0) «(TP2=-TP1) *TERM** 

1 (m2eelLe/XN(N))) *®(TP2-TP1)+ A*D*¥TERM#* (910-01 0¢/XN(N)) ) &( TP3t+TPY) 
RHO(C IJ =(RI LJ) **#242.0*eDRDT (IU) **2-R( IU) *SORDT (IJ) )/ 

1 (R¢CIJ) **24DROT (IU) #2) #*%(32/26) 

IS NOW Ke 

IF (T( IU) *RHOC(IU) eGEe2-0) GO TO 1 

DO 1971 JMA=1e11 

Y=e-T(IJ)*-5 + (CTI(IU)*®.5) * (FLOAT(UMA=1)/5-0 ) 

F (UMA) S=Y¥**2/(1-0*RHO( IJ) *Y) 

CALL YSIMP(10-0/T(IJ)) 

YaNT( IJ) =AREA 

OLV= RHOC(IV)/T(IJ) 

DY 1972 JMA=1e11 

Ye"T(IU)*-5 + (TI(IU)*-5) *© (FLOAT(UMA=1)/5-0 ) 

SUM = (1-+ RHO(IJ)* Y) 

F (CIMA) =DIV*¥*¥242 -0*DIV*Y/ (SUM*¥YINT (IU) ) +Y¥**%2/ (SUM*YINT( IU) ) #2 
CALL YSIMP(10.-0/T(IuU)) 

G(IU)= AREA 

DTVS = le/ SQRT(RI IJ) #*24DRDT (IU) #*2) 

DRUS = DORDT(1IJ)*0TDS 

PV(IQI= RII) €SINITHETACIJ) ) *OTOS=COS( THETAC( IJ) ) *ORDS 
PN(IUI= SINCTHETA( IU) ) *DRDOStR( IU) *COS(THETA( IJ) ) *OTDS 
CONTINUE 
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C NOW 
C WE 


104 


C NOw 


105 


Elastic Ring Optmization Program (continued) 


WE WANT MPRIME(CALLED PM). 

DO SOME SIMPSON INTEGRATION. 

DO 101 JJ=1991 

F(UU)= RIJS) *COS( THETA (UU) ) *SQRT(R( JJ) ##24DROT (J) ##2) #6 (J) /B( Jd) 
1 *(1.¢0-XNHU**2) 

CONTINUE 

CALL SIMP(1.0) 

PAR1 =AREA 

DO 102 JJ=1°91 

F(JU)= SIN(THETA( Ju) ) *DROT (Ud) * RHO( JU) / (BUI) #T (Jd) ) 
CONTINUE 

CALL SIMP(1.0) 

PAR2 =AREA 

00 105) JJ=1"91 

F(Ju)= COS(THETA( JU) )#R( UI) * RHO( JJ) /(B( Jd) #T (Jd) ) 
CONTINUE 

CALL SIMP(1.0) 

PAR3= AREA 

DO 104 JJ= 1°91 

F(JJ)= SQ@RT(R( JJ) #*2 40RUT( UJ) *#2) *6( JU) /B( JU) 
1 * (1 ¢0-XNHU**2) 

CONTINUE 

CALL SIMP(1.0) 

PARY = AREA 

INTEGRATE TO GET W SUB LAMBDA. 

DO 105 JJ=1991 

F(UJI= SQRTIR( Jd) *#24+DRDT (UU) #*2) *B( UU) #T (UJ) #PI/180. 
CONTINUE 

CALL SIMP(1.0) 

PaRS = AREA 


C USE PARS TO COMPUTE w SUB LAMBDA LATER ON. 
C NOW GET PMe 


106 


DO 106 JJ=1791 


PM(UJ)= (PARI“PAR2=PAR3)/PAR4 =R(JJ)* COS(THETA(JJ)) 
CONTINUE 


C FORM NEW LOOP TO GET BMIN. 


JMAX=0 

DMAX =0-0 

DO 8811 IJ =1°91 

DIV = PM(IJ) *RHO(IJ)/(B( IU) *T(IJ)) 

DENOM = DIV= PM(IJ) *T( IU) / (Bb C1) ®(20¢"RHO( IJ) *T (IU) *VINT(IJ)) 
+ PN(LU)/(B(IU)* T(IJ)) 

IF (ABS(DENOM) «GT eOMAX) UMAX=IJ=1 

DMAX = AMAX1 (ABS(DENOM) +» ABS(DMAX)) 

CONTINUE 

BMIN = 2e0/DMAX 

XLAM XLLAM /(BMIN#*S) 

WLAM 4 ¢ *XLAM**#2*GAMMA*PARS 

CUMPUTE SOME VALUES TO COMPARE WITH CONSTAINTS. 

TCONG = XLAM*C5(K)*(1.+ABS(C8(M))) 


JRL. RES. SC. © 


154 





Elastic Ring Optimization Program (continued) 


TCONS = XLAM*(2-*C1(1)+T(1)) 
TCON6 = XLAM*(2.*C2(1)+T(91)) 
TCON7T = C3(U) *(1Le+ABS(C7(L))) 
TCONS = XLAM*(26*C2(1I)-T(91)) 
GO TO DOOPS IF CONSTRAINTS ARE NOT SATISFIED. 
DV 107 JJ=1°91 
Tm = (Le —XNHU**2) *G (UJ) *T (UJ) PM (Ud) &R CU) #COS (THETA (JU) ) 
tie: = RHO (JJ) #PM( JJ) #PN( JJ) 
DS= SQRT(R (JJ) **2+DRDT (UJ) **2) 
TS = RHO (JU) *PN( UJ) #R (JU) *COS (THETA (UJ) ) 
F( JJ) =10/(B( J) *T (Ud) ) ® (PN (UU) #2412 6/56*(1¢+XNHU) * 
1 PV( GJ) **2-TM1+TM2=-TM3)*DS *PI/180. 
107 CONTINUE 
CALL SIMP(1.0) 
ASTAR=AREA 
Uc LAM=XLAM*ASTAR*BMIN*S/E 
RaT LO=WLAM/DELAM 
C COMPUTE SOME OUTPUT ITEMS. 
IF (MOD (UCOUNT*® 25) eNEe0) GO TO 7 
WRITE (604) 
ue FORMAT(IH1 1X *WT/DEFL'4X* WEIGHT '6X'G1*6X'G2'5X"G3'6xX'Gu4r6x'Gs? 
1 T7XAK*GO*3X*LAMBDAt4xtAt7TX*C*EX* Dt 7X"E*TX*EPTX*NOOXtTMt//) 
7 WRITE (605) RATIO*e WLAMe DELAM?e TCONG se TCONSe TCONGe TCON7# TCONS? XLAMe 
2 CLI) 0 C364) 0 CS(K) 0 C7(L) 0 CB8UM) & XNCN) » UMAX 
> FORMAT (1HO*?2F9e20F9e3012F 82350135) 
JCOUNT=JCOUNT 41 
1 CONTINUE 
IF (IGOeNE-0) GO TO 44 
C wE HAVE WOW FINISHED MAJOR LOOP. 
STOP 
ENO 
WTI FOR SINTe SINT 
SUBROUTINE SIMP (HHH) 
CUMMON /ZINTEG /F (100) AREA 
H = 1¢/(30¢*HHH) 
AnEA=0-0 
QuLD=0-0 
EVEN=0.0 
DY 21 I =209002 
EVEN =EVEN + F(T) 
DU 22 I =3089r2 
Obv =O0DD+ F(I) 
AKEA = H* (F(1)4F(91)44.*EVEN+ 2-*0DD) 
RETURN 
ENO 
@TI FOR YSUBs YSUB 
SUBROUTINE YSIMP (HHH) 
COMMON /INTEG /F (100) »AREA 
H = Le/(30e*HHH) 
AREA=0.0 
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i XaT 
1 
1-21 
12.7 
e305 
-.025 
oh 
6e 
e350 
@ EOF 
@ FIN 


Elastic Ring Optimization Program (continued) 


OvD=0-.0 

EVEN=0.0 

DO 21 I=201002 

EVEN =EVEN + F(I) 

DO 22 I=39992 

Oud =O0D0+ F(I) 

AREA = H¥*¥(F(1)4F(11)4+4e0*EVEN+2.0*0DD) 
RETURN 

ENO 

MaIN1 

5 i 6 3 0 


13. 13.3 14. 


0. e025 075 el 
05 


1000000. 150000. 30000000- 
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Detailed Ring Analysis Program 


id RUN MITCHEe? 16650702760 
@WiT- FUR MAIW2eMAIN2 
COMMON /INTEG/F (100) »AREA 
DIMEWSION XLN(91) eXLV(91) »XLM(91) 
DIMENSION YINT(91) 
DIMENSION PVv(91) 
DIMeNSION 6191) 
OIMENSION B(91)0T(91) eRI91) eDRDT(91) eSDRDT(91) eRHO(91) » 
1 PM(91)¢PN(91) eSIGTL (91) e THETA(91) 
10 READ(5220) CleC3eCSeC7eCBeXN 
20 FURMAT (10E6.0) 
IF (C1eEQ29999-) STOP 
REA (5°30) XNHUeXLLAMeSeEe GAMMA 
30 FURMAT (5E15.0) 
C INSERT INITIAL PRINT HEREe 
C PRINT OUT ALL INPUT DATA AT START OF RUNe PROGRAM 3. 
WRITE(621000) CieC3eCSeC7eCBeXN 
1000 FURMAT (. IH1le SX*#2HC1*8Xe2HC39 BX 2HC5e 8X0 2HC7 ee BXe2HCBe BXe LHNe// 
z oF 10.3) 
WRITE (601070) 
WRITE (601080) XNHUseXLLAMeSe Ee GAMMA 
C INSERT FORMATS 1070 AND 1080 HERE. 
1070 FURMAT(1HO 10X°3H NUe15Xe 7H LAMBDA? 18X°1HS01 9X0 LHE® 14XeSHGAMMA ) 
1080 FORMAT (1H0O+5E20-8) 
Pl =32141592654 
C2=2.0-Cl 
Cy = C3*C7 
Co = C5*C8 
A =Cl 
D=Cz2 
Dv 2 I =1°91 
JI = [el 
THETACI)SFLOAT(JUI) *PI/180. 
B(I) = C3 * (1e"C7 *COS(2-*THETA(I))) 
Til) = C5 * (1e"C8 *COS(2-*THETACI))) 
TERM = A ** XN *SIN(THETA(I))**XN#D *&*XN * 
COS(THETA(I) )** XN 
R(I)SA*D/( TERM **(1¢/XN) ) 
TPL = D **#XN* COS(THETA(I) ) **®(XN=120) *SIN(THETA(I)) 
TP2 = A **XN* SIN(THETACI) ) **(XN=1-0) *COS(THETA(TI) ) 
DRDT(II=SR(I) *(TP1=-TP2)/TERM 
IF (XNeEQe2eQeANDelI eEQel1) GO TO 12l 
TPS = A **#XN*(SIN(THETA(I) )*®*® XN@(XN=-1.0)% SIN (THETA(I)) ¥* 
1  (XN"2-0)*COS (THETA(1))¥**2 ) 
IF (UXNeEQ@e2eOeANDel e£Qe91) GO TO 122 
GY TO 125 
121 TRS = A*x*2*(-1.0) 
4125 CONTINUE 
TPS = D¥*XN *(COS(THETA(I) )*®*® XN=(XN"1-0)*% COS (THETACI)) ** 
1 (XN72-0)*SIN (THETA(I) )**2 ) 
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Detailed Ring Analysis Program (continued) 


Gu TO 123 

TP = D*¥*2 * (—1.0) 

CONTINUE 

SURVT(I)= A*eDd * (XN#1-0)* (TP2=TP1)* TERM** 
1 (#20 Le/XN)) *(TP2=TP1)+ A*D*TERM#*(—1.—(16/XN) ) *#(TP3+TPY) 
RHOCLT) = (RCI) **242.*DRDT(1) **2=R(1) *SDROT(I)) 
1 / (RCI) ®*24+DROT(I) **2) **#(320/2-0) 

DU 1971 J=l1eil 

Y=-T( 1) *.5+ TI) *e5*«(FLOAT(U-1)/5.0) 

Fi J) =¥**2/(1.0+RHO(I) *Y) 

CALK YSIMP(10.0/T(I)) 

YaNNTCIV=SAREA 

DIV = RHOC(I)/T(I) 

DU 1972 J=lell 

Y=-T(1) #254 TUL) *e5*(FLOAT(U-1)/5.0) 
SuM=(1-¢0+RHO(I) *Y) 

F (J) SDI V¥*24+2¢*DIV*Y/ (SUM*YINT(I) )4Y¥*#%2/(SUM*#YINT(I)) #*2 
CALL YSIMP(10-0/T(I)) 

Gil) = AREA 

UTOS =Le/SQRT(R(I)**#2 +DRDT(I) **2) 

DRUS = DRDT(I)* DTOS 

PV(L)= RCI) *®SIN(THETA(I) ) *DTOS=COS(THETA(I) ) *DRDS 
Pili) = SIN(THETA(I))*DRDS #R(I)*COS(THETA(I)) *DTDS 
COWT NUE 

WE INTEGRATE TO GET MPRIME (CALLED PM). 

DY 1v1 J =1°91 

F(J)=R(U) *COS(THETA( J) ) *SQRT(R( JU) #*24DRDT (J) *#2) 
1 *6(JU)/8(40) *(10¢0—XNHU**2) 

CONTINUE 

CALL SIMP(1.0) 

PaRL=AREA 

OG 102 J =1°91 

FiJ)= SIN(THETA(J) ) *DRDT (JU) *RHO(J)/(B( JU) * T(J) ) 
CONTINUE 

CALL SIMP(1.0) 

PAR2 =AREA 

DY 103 J=1°91 

FiJ)= R(J)*COS(THETA(J) ) *RHO( J) / (BCU) #&T(U)) 

CONT INUE 

CALL SIMP(1.0) 

PAR3= AREA 

00 104 J=1991 ; 

F (J) =SQRT(R( JU) **24DRDT(U) **2) *6(U)/B(U) # (1. 0=XNHUF*2) 
CONTINUE 

CaLL SIMP(1.0) 

PARGY = AREA 

woe INTEGRATE TO GET w SUB LAMBDA. USE PARS LATER ONe 
Du 1iuS J=1°91 

F(J)=SQRT (RI U)#*2 *tDROT(U)**2)% B(U) *T(U) *P1/180.0 
CONT ANUE 
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Detailed Ring Analysis Program (continued) 


CALL SIMP(1.0) 
PARS = AREA 


C NOW GET PM. 


106 


DO 106 J=1°91 
PM(J)= (PARL=PAR2=PAR3)/PAR4Y =R(JU) *COS(THETA(J) ) 
CONTINUE 


C FORM LOOP TO GET BMIN AND SIGTL. 


6811 


107 


C NOW 
C NOW 


40 


70 


DMAX = 020 

DO 8811 I=1°91 

DIV = PM(I)* RHO(I)/(B(I) *T(I)) 

SIGTL(I)S DIV = PM(I)*T(I)/(B(I) &(2e-RHO(I)* TCI) ) *YINT(I)) 
1 +PN(I)/(B(I) *T(I)) 

DMAX = AMAXL(ABS(SIGTL(I1) ) eABS(DMAX) ) 

CONTINUE 

BMIN = 2¢0/0MAX 

XLAM= XLLAM /(BMIN®S ) 

FACTOR=BMIN*S¥*05 

DO 6833 I =1°91 

SIGTL(I) =FACTOR*SIGTL(I) 

XLN( I) =XLAM*FACTOR*PN(I) 

XLV(1)=XLAM*#FACTOR*PV (I) 

XLM(1)=XLAM#*2*FACTOR#*PM(T) 

CONTINUE 

DO 107 JJ=1991 

TML = (1¢"XNHU**2) *G (JJ) *T (UU) *PM (JU) *&R (UU) *COS ( THETA (Ju) ) 
T2 = RHO (JJ) *PM( JJ) *PN (UU) 

DS= SQRT(R( JU) **2+DRDT (UJ) *#2) 

TMS = RHO (JJ) #PN( JJ) *R (UU) *COS (THETA (UJ) ) 

Fi JU HL e/ (BCU) *T (UU) ) #(PN( UU) ®*#24120/56%(106+XNHU) * 

1 PV( UJ) **2=TM1L+TM2-TM3)*DS *PI/180. 

CONTINUE 

CALL SIMP(1.0) 

ASTAR=AREA 

DELAM=XLAM*ASTAR*BMIN*S/E 

WLAM = 4e* XLAM **2*% GAMMA *PARS5 

COMPUTE SOME OUTPUT ITEMS. 

TCONG=XLAM*C5%(1-+ABS(C8) ) 

TCONS= XLAM * (26% A +#T7(1)) 

TCONW6= XLAM * (2e* D +T7(91)) 

TCON7T= C3 * (1-+ ABS(C7)) 
TCONS=XLAM*(2-0*C2=T(91)) 

COMPUTE STRESSeFOR EACH THETAe 

WE PRINT OUT DESIRED OUTPUT» THEN READ ANOTHER DATA SET. 
WI TE (6040) 

FORMAT(1HO 1X "WI/DEFL®4X "WEIGHT *6X°G196X"*G2°5X"G3"6X"*G4t6x'Gs? 
lL 7X*G6*3X*LAMDAtSX*At7TX*CPex* Dt 7X*E*7TXtFPTXING/S ) 
RATLO=wWLAM/DELAM 

WRITE (6070) RATIO#e WLAMe DELAMe TCONY es TCONS# TCONG? TCONT# TCONS? XLAMe 
1 CieC3eCSeC7eCBeXNn 

FORMAT (1HOe 2F9e20F9.3012F 8-3) 
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Detailed Ring Analysis Program (continued) 


WRITE (6080) 
80 FURMAT(LHO* THETA" 6OX*R*OX*LAMR*9OX*T*BX*LAMT*9X*LYOX*K*OX* INT? 
1 7X*LAMN®7X*LAMV "GX *LAM*¥*#24M"6X"*SIGMAt//) 
DO 64% I =1°91 
J =I-1 
RLAM=R (I) *XLAM 
TULAM=T( 1) *XLAM 
WRITE (6090) UeR(I) eRLAM@ TCI) e TLAMeB( I) pRHOC I) eVINT(CI) @XLNCI) 
Ll ALV(I) eXLM(I) eSIGTL(I) 
90 FORMAT(IH 15011E11-4) 
IF (MOD (495) eEQ@eOeANDeJeNEoO) WRITE(60110) 
64 CONTINUE 
FORMAT (1H ) 
GU TO 10 
STOP 
End 
@N FOR SINT*SINT 
SUBROUTINE SIMP (HHH) 
COMMON /INTEG /F (100) ¢AREA 
H = 1-0/(3-0*HHH) 
AKEA=0-0 
Qv0=0 e 0 
cEVEN=0-0 
DO 21 I =2°9002 
Even =EVEN + FCI) 
DO 22 I =3069r2 
Ouuw =ODD+ FCI) 
ARCA = H*® (F(1)4+F(91)44.*EVEN+ 2-*00D) 
Re TURN. 
ED 
@N FOR YSuUBe YSUB 
SUBROUTINE YSIMP (HHH) 
COMMON /INTEG /F (100) AREA 
H = 1¢/(3¢*HHH) 
AKEA=0.0 
Ouv=0-0 
EVEN=0 0 
Dv 21 I =201002 
EVenw =EVEN + F(T) 
DO 22 I =30 9e2 
Ovi =O0DD+ FCI) 
ARCA = H* (F(1)4F(11)44.8EVEN+ 2.*000) 
RETURN 
EnD 
@ XaT MAIN2 
1-22 13.8 e305 Oe 045 5.5 
e3 1000000. 150000. 30000000~- 
le 14, e248 O- Oo 26 
ed 1000000. 150000- 30000000- 
i OF 


JRL. RES. SEC. C 


(Paper 72C2-273) 





JOURNAL OF RESEARCH of the National Bureau of Standards —C. Engineering and Instrumentation 
Vol. 72C, No. 2, April-June 1968 


Applications of Coaxial Chokes to A-C Bridge Circuits 


D. N. Homan* 


Institute for Basic Standards, National Bureau of Standards, Boulder, Colo. 80302 
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One or more turns of coaxial cable on a magnetic core form a coaxial choke which may be used to 
suppress unwanted loop currents in a-c bridge circuits. A review of several applications of coaxial 
chokes to a-c bridge circuits is presented. These applications include suppressing ground loop currents, 
assuring unique definition in the measurement of small 3-terminal capacitors, reducing stray magnetic 
field and reducing voltage drops in leads caused by large ground capacitance current. There is also a 
discussion of techniques for greatly increasing the effectiveness of the coaxial choke. 


Key Words: Bridge, capacitance, coaxial cable, ground loop, inductance, magnetic core, stray 


magnetic field. 


1. Introduction 


Large improvements in stability have recently been 
made in commercial 3-terminal capacitance standards, 
and transformer-ratio-arm bridges of high precision 
have become available. The availability of such equip- 
ment may result in measurements of better precision 
in a greater number of laboratories. The accuracy of 
such measurements can frequently be improved by 
the use of coaxial chokes in appropriate parts of the 
bridge circuit [1, 2].! Coaxial chokes have been used 
by workers in other fields as well—sometimes under 
different names, such as “cable chokes” [3] and 
“baluns” [4]. This paper gives a review of several 
applications of coaxial chokes as applied to a-c 
bridge measurements 

Often a-c bridge circuits have components inter- 
connected with coaxial cables to prevent undesired 
coupling between different parts of the circuit. Fre- 
quently, the physical layout permits more than one 
return path for the current to a given component. 
Such multiple returns have the effect of producing 
currents in the outer conductors of the coaxial cable 
which are not equal and opposite to the corresponding 
inner conductor currents. This departure from equality 
can be called a “net current.” The elimination of the 
net current is a major problem in precision a-c bridge 
measurements, and one technique, that of using 
coaxial chokes, will be discussed here. 


*Radio Standards Engineering Division, HF Impedance Standards, NBS, Boulder, 
Colo. 80302. 


' Figures in brackets indicate the literature references at the end of this paper. 


A coaxial choke is formed by winding one or more 
turns of coaxial cable around a high permeability 
toroidal magnetic core, as illustrated in figure 1. A 
higher impedance to net current can be obtained either 
by increasing the number of turns through the core or 
by increasing the size (or permeability) of the core. In 
the former case, the increased cable length results in 
increased cable impedance which may be undesirable. 
In the latter case, the increased cost may be undesir- 
able. The self-inductance of either conductor of the 
coaxial cable is given approximately by 


4X 10-2 yA 
s 





where L is inductance in henrys, N is the number of 
turns, « the permeability, A the cross-sectional area 


FIGURE 1. Coaxial choke (coaxial cable wound on a high 


permeability core). 
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FIGURE 2. Capacitance measuring system showing ground loop. 


of the core in cm?, and S the mean core length in cm. 
A typical range of inductance is 10 to 40 mH at 1 kHz, 
for example. This inductance reduces only the net 
current in a coaxial cable; the equal “‘go and return” 
currents induce equal and opposite magnetic fluxes 
in the core. Thus, the contribution to the equivalent 
self-inductance from the high permeability core is zero 
for equal “go and return” current. 


2. Applications 
2.1. Ground Loop 


If in the capacitance bridge illustrated in figure 2 [5] 
the connection OD were not present, the detector 
would sense Vz, the bridge unbalance voltage. How- 
ever, if the case of the detector electrically contacts 
the case of the oscillator (connection OD), a ground 
loop is formed, accompanied by a possibility of an 
error in the measurement. The contact between oscil- 
lator and detector may be through a relay rack, or 
may result from both oscillator and detector being 
connected to a common power source (a-c line for 
example). Since z’ is generally not zero, a portion of the 
current i may be diverted through z, such that i, ap- 
pears in the detector ground lead, resulting in a voltage 
drop V. In this case V, is an error voltage, adding to 
the bridge unbalance voltage Vz. When the detector 
indicates null, V.=V,, thus giving an erroneous 
result. In addition, the current i, in the loop completed 
by OD will generate a magnetic field which may couple 
into other parts of the circuit, or the loop may have 
voltage induced in it by other sources of magnetic 
field. 

The error can be reduced by the use of low imped- 
ance leads and connectors from the bridge to the 
detector and generator, e.g., coaxial cable with low 
impedance in the shield. Also, the ground loop can be 
avoided by using an insulated, battery-operated detec- 
tor, or an isolation transformer can be used in an ap- 
propriate place in the circuit. However, in cases where 
it is not convenient or possible to use an insulated, 


é 


FIGURE 3. Partial circuit showing detector cable passed 
through a high permeability toroidal core. 


battery-operated detector, or an isolation transformer, 
a coaxial choke in the coaxial cable between the bridge 
and the detector may achieve the desired result. A 
coaxial choke can also serve as a convenient test for 
trouble. A change in the bridge balance accompanying 
insertion of a coaxial choke indicates that a problem 
previously existed. 

Figure 3 illustrates the use of the coaxial choke in 
the detector cable. It may be useful first to think 
of the toroid being around the connection OD in posi- 
tion Q. The single lead choke thus formed is a high 
impedance relative to z, so that i, in figure 3 is much 
less than i, in figure 2 (where no choke is used). 
Thus, the error voltage is considerably reduced; for 
example, if z=z’=0.019 and the impedance of the 


choke is 1000, then 2X 10-4 One can now 


2 
imagine pushing the magnetic toroidal core into posi- 
tion P around the detector cable; the two positions 
are electrically equivalent. 

It may be useful, as a rule of thumb, to consider the 
use of a coaxial choke in every unwanted loop in the 
measuring circuit. For N paths to some reference 
point it may be necessary to use N-] coaxial chokes. 


2.2. Measurement of a 3-Terminal Capacitor 


A 3-terminal capacitor connected to a bridge is 
illustrated in figure 4. Here C; is the capacitance to 
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FIGURE 4. 3-Terminal capacitor connected to a bridge with leads 


having inductances L,, Lz and Ls. 
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L3 
FIGURE 5a. 3-Terminal capacitor connected 
to a bridge showing the use of a coaxial choke. 


i] 


be measured; C, and C3 are capacitances to ground; 
L,, Lz, and Ls are lead inductances. It has been shown 
[6] that the effective capacitance seen by the bridge 
is given by 


Cote = GC. 1 _ W°C2Cs by + wC2L, + wC3L 
1 


+ wl; (Le + Ls) 


If C, is small, relative to C. and C3, then 


: . w?C2C3L 
Cote = C; E |, 


a can be 10-> if C,=1 pF, 


1 
C.=C;= 1000 pF, @>= 104, and L,=0.1 wH. At point 


For example, 


1 
A, iz is divided between jwl, and om Thus the 
voltage developed across L, by the current is drives 
current i; into C3, and to the detector, though it did 
not pass through C,, which is the capacitance being 
measured. The capacitor’s apparent value therefore 
will depend upon the circuit in which it is measured, 

and thus will not be defined uniquely. 
Let us assume that the leads from the active elec- 




















FIGURE 5b. Transformer bridge with coaxial 
choke. 








FIGURE 5c. Bridge split to 
show effect of the coaxial 
choke. 


trodes of a 3-terminal capacitor are attached to the 
inner members of coaxial fittings, and that the case 
of the capacitor and the outer members of the coaxial 
connectors form the third “‘terminal.’”’ Also assume 
that such a capacitor is connected to a bridge with 
two coaxial cables. The effect that may cause an 
error (discussed in the preceding paragraph) can be 
nearly eliminated by inserting a coaxial choke into 
a cable connected to the capacitor as illustrated in 
figure 5a. The choke is a high impedance to net cur- 
rent (i—i’). Thus most of the current in C. returns 
to ground through L;, rather than through C; and L;. 


Figure 5b shows a transformer-ratio-arm bridge as 
a practical example. Here C2, C3, C,, and C; are ground 
capacitances. Lead inductances are also shown. The 
result of using the coaxial choke may be thought of as 
breaking the connection BB’ in figure 5c and inserting 
the high impedance, Z,. The impedance of the ground 
loop is greatly increased by the core. 


Figure 6 illustrates how the small inductance, /,, 
of the case of the capacitor itself, in conjunction with 
ground capacitances C, and C3, may cause the capaci- 
tance being measured to differ from the value calcu- 
lated from the capacitor’s dimensions. This difference 
would remain even if a coaxial choke of infinite im- 
pedance were used; it would be considered an error 
only in the case of a calculable capacitor and has, for 
example, been evaluated for the NBS cross capacitor 
[6]. Mention of it is included here only to show that not 








FIGURE 6. Partial bridge circuit to illustrate |, 
a small inductance which is part of the capacitor, 
not the connecting leads. 
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FIGURE 7. Circuit using bypass leads and coaxial chokes. 


all errors of this type can be eliminated by simply 
using a coaxial choke. Current iz is divided at A’ 


between jal, and jal 


2.3. Stray Magnetic Field 


A net current in any given coaxial cable of a bridge 
circuit flows in other conductors of the circuit. Be- 
cause of the physical layout these other conductors 
have the effect of creating loops that are large relative 
to loops designed into the circuit. The magnetic field 
associated with the currents in such loops can cause 
trouble. A coaxial choke may be used to suppress the 
net current, thereby reducing the stray magnetic 
field. 

However, if the net current to be eliminated is in- 
jected through a very high impedance, such as a small 
capacitor, the coaxial choke may be a relatively small 
impedance to the net current and nearly ineffective. 
The coaxial choke is an effective impedance to net 
current only when the net current is a result of a cur- 
rent divided by lead impedances that are small rela- 
tive to that of the coaxial choke. 


2.4. Special Uses 


Figure 7 illustrates a technique used by Shields [7] 
to eliminate errors due to voltage drops in critical 
parts of his circuit. The voltage drops in the circuit 
were caused by the loading of capacitances C3 and C,. 
These are capacitances to ground from circuit shield- 
ing that was at a‘high potential above ground. The 
errors were eliminated by supplying current to the 
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FiGURE 8. Use of triaxial cable and two coaxial chokes to 
multiply the effectiveness of the coaxial choke. 











FIGURE 9a. “Triaxial choke” system with the outer 
shield completely enclosing the inner shield. 


capacitances by low-impedance “‘bypass leads,” and 
coaxial chokes were used in the coaxial cables as 
shown to restrict the current to the bypass leads in 
which voltage drops were not critical. 


2.5. Extra High Impedance Coaxial Chokes 


A coaxial choke can be modified to multiply its 
effectiveness by adding a shield and a magnetic core 
as shown in figure 8. Figure 9a is an equivalent circuit 
of figure 8. A calculation shows that such a “triaxial 


Zz ; 
choke” can reduce net current by Rp where z is the 


impedance of either shield and Z is the impedance of 
each choke as shown in figure 8. Since a single core 
coaxial choke reduces net current by approximately 


< : 
74 large improvement can be expected from the 
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FIGURE 9b. Equivalent circuit of a “triaxial choke” system where 
outer shield is connected short of the ends of the inner shield. 
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FIGURE 10. Manual high impedance coaxial choke. 








“triaxial choke.” To realize the maximum improve- 
ment from the “triaxial choke,” care must be taken 
to connect the outer shield to the ends of the inner 
shield. If the outer shield is connected to the inner 
shield at points short of the ends, an equivalent circuit 
shown in figure 9b can be used. A calculation based 
on this equivalent circuit shows that net current is 
zw zetZ 


reduced approximately by nt Zz where the 


assumptions made are z3;=2; =Z, joL, =jol 2 =Z, 
and z<Z. 


Two additional techniques [8, 9] for increasing the 
effectiveness of coaxial chokes also involve two cores, 
one a sensing core and the other a voltage injection 
core. One technique is to perform these functions 
manually as shown in figure 10. The detector senses 
flux in the core from net current in the coaxial cable. 
The generator output level can be adjusted such that 
the detector indicates null. The other technique is 
illustrated in figure 1] where the same functions are 
performed electronically. If ep=Zpi, eg=Gep, and 
er=e,p+eg, it can be shown that Z7-Z,(1+6G). The 
gain G of the amplifier and the number of turns through 
the cores can be adjusted to reduce net current by an 
amount which will make the error current negligible. 
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FIGURE 11]. Active high impedance coaxial choke. 
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Creep Rupture Properties of Ti-8Al-1Mo-1V Alloy 
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Creep tests were made at 600, 800, 1000, and 1200 °F (588, 699, 811, and 921 °K) on duplex-annealed 
Ti-8Al-1Mo-1V alloy specimens. Creep-test conditions were selected to produce creep rates ranging 
from 0.1 to 50,000 percent per 1,000 hours and rupture times ranging from 0.1 to about 10,000 hours. 
Conformance of the test data to creep theories was obtained over a limited range of stress, strain, creep 
rate and temperature. Creep behavior was influenced by microstructural changes during creep and by 
prior thermal-mechanical treatment. Engineering-design curves are included in this paper. 


Key Words: Creep, engineering design, high temperatures, phase changes, rupture, titanium alloy. 


1. Introduction 


During the past several years a comprehensive 
program has been in progress at the National Bureau 
of Standards to evaluate factors influencing the me- 
chanical behavior of metals at elevated temperatures. 
Concurrent with the presentation of engineering data, 
obtained in tensile and creep tests made on a variety 
of metals, the mechanical properties were interpreted 
in terms of modern theories of plastic deformation. The 
validity of most of the theories were confined to a lim- 
ited range of variables. As the status of current theory 
of creep behavior is well documented [1, 2],! an effort 
will be made in the present paper to omit the deriva- 
tions of the formulas and to apply only those concepts 
pertinent to the data obtained on a titanium alloy. The 
variables considered here are stress, temperature, 
creep rate, ductility, rupture time, microstructural 
changes, and prior strain history. Results are pre- 
sented for creep tests made on the alpha-beta duplex- 
annealed alloy, Ti-8Al-1Mo-1V, at temperatures of 
600, 800, 1000, and 1200 °F (588, 699, 811, and 921 °K). 


2. Material and Testing Procedure 


The titanium alloy, Ti-8Al-1Mo-1V, selected for 
this investigation is one of the structural materials 
intended for use in supersonic transport aircraft. 

The test specimens were prepared from bars ob- 
tained from the same heat as that used in a previous 
research program to evaluate the effects of notch 
geometry on the tensile behavior of this alloy [3]. 


' Figures in brackets indicate the literature references at the end of this paper. 


The bars were supplied by the manufacturer in the 
form of 1-in rounds in the hot-rolled and annealed 
condition. The chemical composition of the alloy is 
given in table 1. 


TABLE 1. Chemical composition (percentage by 


weight) of the alloy used* 








Titanium 
Aluminum 
Molybdenum 
Vanadium 


0022 





4 Analysis made at the National Bureau of Standards. 


The bars were cut into lengths, each about 6.5-in 
and annealed at 1750 °F (1226 °K) for 1 hr and air- 
cooled. This was followed by a stabilizing anneal for 
8 hr at 1050 °F (838 °K). The specimens were then 
air-cooled. This process, called duplex-annealing, 
is recommended for obtaining a good combination of 
properties at high temperatures. 

Each unnotched specimen had a 2-in gage length 
and a gage diameter of 0.357 in. The reduced section 
was finished by grinding and polishing the axial 
direction. 

Before testing, each creep specimen was heated in 
a tube furnace to the desired temperature and equili- 
brated at temperature for 16 hr. It was then tested at 
temperature under constant load. The creep-rupture 
machines used were of a multiple lever type with 
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motorized jacks to prevent shock loading. They had 
been calibrated at room temperature and the errors 
for the stresses used were less than 0.5 percent. The 
temperature of each of the specimens was controlled 
to an accuracy of +2 °F. Strain-time data were ob- 
tained from electric contact follow-up type extensom- 
eters attached to shoulders of the specimens. The 
sensitivity of extension measurement was approxi- 
mately 0.00005 in. Metallographic examinations were 
made on specimens after fracturing in creep to assess 
the influence of the test variables on changes in 
microstructure. 

In order to permit selection of creep stresses to be 
used in this investigation, short-time tensile tests were 
made in duplicate at the same temperatures as those 
later used in the creep tests. The average tensile 
values are shown in table 2. 


TABLE 2. Short time tensile properties of the duplex-annealed alloy 








Yield 
strength, 
0.2% 
offset 


Reduction 
of area 


Tensile Elongation 


-strength 





Ksi 
135 
87 
80 
65 
45 

















3. Results and Discussion 


3.1. Creep Test Data 


Extension-time data were plotted on a linear scale 
for all the specimens tested under various conditions 
(not shown). All the curves produced in this manner 
had the same general shape in that the first, second, 
and third stages of creep were clearly discernible. 
Loci of beginning of the third stage of creep, loci of 
fracture, and second stage creep rates were determined 
from these curves. However, no simple translation 
process could make the curves coincide for tests made 
at the same temperatures with different stresses or at 
different temperatures with equal stresses. J. 
deLacomb [4] has suggested that creep curves of this 
type could be described by the sum of two parabolic 
curves. Thus, 


=€)+at™+ bt" (1) 


€=creep strain 
€9 = instantaneous strain 
= time 
a and b= constants 
m<landn>l. 
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FIGURE 1. Extension-time relations for specimens tested at 1000 °F. 


Extension-time curves for specimens tested at 
1000 and 1200 °F in the present investigation are 
shown in figures 1 and 2. These data, as plotted on 
a log-log scale, indicate a qualitative conformance to eq 
(1) only over limited ranges of time if € is considered 
to be equal to zero. However, from a design engineer- 
ing standpoint, the 0.1 hr and also the subsequent 
extensions are readily observable. Moreover, the loci 
of the beginning of the third stage indicate aminimum 
at 1000 °F, whereas both a minimum and maximum 
are present at 1200 °F. The loci of fracture at each 
temperature are quite similar in shape. 


A number of formulas have been derived, empirically 
and theoretically, describing the dependence of 
second-stage creep rate or fracture time on creep 
stress [l, 2]. Generally, it has been found that at 


% 
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FIGURE 2. Extension-time relations for specimens tested at 1200 °F. 
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FIGURE 3. Effect of stress on second-stage creep rate of specimens 
tested at different temperatures. 


“low” stresses and “high” temperatures the stress 
dependence is given by 


€;= Ao" (2) 
or 
=A,o7-"1 (2a) 


where 


€; = second stage creep rate 

ty=rupture time 

o=creep stress 

n,n, A, B, A,, and B,= constants which 
are. dependent on temperature. 


€,=A exp (Ba) (3) 


and 


ty =A, exp (—B,o). (3a) 
However, at “high” stresses and “low” temperatures 
the stress dependence is given by the formula first 
proposed by Kanter [5]. Equations (3) and (3a) indi- 
cate that there should be a linear relation between 
stress and either log second stage creep rate or log 
rupture time. Data obtained on the Ti-8Al-1Mo-1V 
specimens are plotted as a function of stress in fig- 
ures 3 and 4. Although the equations indicate that in 
each of these sets the two curves should be mirror 
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FicuRE 4. Effect of stress on rupture time of specimens tested at 
different temperatures. 


images of each other, quantitative evaluation of the 
data showed that, at each temperature, the slopes 
of log t- versus €,; curves (not shown) were differ- 
ent at each temperature and at each stress level. 
This observation indicated that'the mechanism of 
deformation at these temperatures is different as 
the test parameters are changed. However, the 
curves in figures 3 and 4 do show that the strength 
of the alloy is more rate-sensitive at the higher tem- 
peratures than at the lower temperatures. Additionally, 
the rate-sensitivity at the different stress levels is 
an inconsistent index of the creep behavior of the 
metal. 

The variation of stress with temperature to produce 
engineering-design curves (A) rupture at various times 
and (B) various second-stage creep rates is shown in 
figure 5. The data are derived from figures 3 and 4. 
The shape of these curves indicates the dangers of 
predicting behavior of metals for long-time service 
from data obtained in short-time tests. In addition, 
the slopes of these stress-temperature curves indicate 
continuously changing mechanisms of deformation and 
fracture at the different temperatures. The similarity 
of families of curves of the type shown in figures 5A 
and 5B led Monkman and Grant [6] to propose a rela- 
tion of the form 


t-= Ke-™ 


log t-=m, log “—" log K 
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FIGURE 5. Variation of stress with temperature required to (A) 
cause rupture at various times and (B) produce various second 
stage creep rates. 
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FicuRE 6. Relation stress 


Tx = Temperature in degrees Rankine 
t= Rupture time in hours. 


t-= rupture time 

€;= second stage creep rate 
m,= slope of log €; versus log t, 
K= constant. 


However, a plot of the product (t,€,) versus tempera- 
ture did not yield straight lines for the data in this 
investigation. 

A number of methods have been used for extrapo- 
lation of creep data by use of parametric formulas 
[1]. A time-temperature parameter derived by Larson 
and Miller [7] and based on the general Arrhenius 
rate equation is of the form 


P,=f(o) =Tr(log t+ As) (5) 
where 
P,= parameter 
o = stress 
Ty =temperaiure in degrees Rankine 
t=rupture time 
A,= constant (generally 20), 


Conformance of the present data to this expression 
is shown in figure 6. Although a single curve can be 
roughly drawn through the experimental points, it is 
impossible to predict in advance the shape of this 
curve. Until some parametric expression can be 
derived which will include metallurgical changes 
affecting ductility at the different temperatures, 
these formulas will be considered useful only for 
cataloging engineering data. 

Most theoretical analyses of creep phenomena indi- 
cate that there should be a continuous change in rup- 
ture ductility values with changes in stress or rupture 
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FIGURE 7. Effect of test time on elongation and reduction of area 
values for specimens ruptured at different temperatures. 


time. Examples in the literature have shown that the 


validity of this concept is often affected by precipi- 
tation [8], recrystallization [9], grain boundary weaken- 
ing [10] and other metallurgical phenomena. The 
relation between rupture-time and elongation and 
reduction of area values of the titanium alloy specimens 
used in the present investigation is shown in figure 7. 
Anomalies in the curves may be attributed to recovery, 
internal cracking, and the dissolution of the beta 
phase during the long-time tests. 

The interrelationship of of the variables associated 
with creep testing [1| generally presupposes the ex- 
istence of an equation of state for creep of the general 
form 


o=f(é,, T. €) (6) 


o=fi(t,, T, €) (7) 


o= creep stress 

€; = creep rate 

T= absolute temperature 
€= strain 

ty= rupture time. 


These equations indicate that creep behavior is inde- 
pendent of the prior thermal-mechanical history of 
the metal. This is not borne out by figure 8, which 
shows the effect of prior history on the subsequent 
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FIGURE 8. Effect of prior strain history on extension-time relations 
for specimens tested at 1000 °F with a stress of 30,000 psi. 
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strain-time relations of specimens tested in this 
investigation at 1000 °F with a stress of 30,000 psi. 
The extension values are the total extensions obtained 
during testing under each condition; the stress value is 
based on the original cross-sectional areas; and the 
time is the time at 1000 °F and 30,000 psi only. At 
all temperatures within the limits of this investiga- 
tion, prior creep decreased the elongation at rupture. 
However, the rupture time was increased with increase 
in the temperature at which the specimens were de- 
formed. Although prior thermal history was an in- 
consistent index of the elongation values, it appears 
that increasing the prior extension at constant tempera- 
ture (curves D, C, and B) caused a decrease in the 
elongation values. Within the range of temperatures 
used these data prove conclusively that the creep 
behavior of this titanium alloy is strongly dependent on 
its thermal-mechanical history. 


3.2. Metallography 


Photographs and micrographs of selected speci- 
mens after rupture in creep are shown in figures 9 


through 13. 


Surface appearances of some of the creep specimens 


are shown in figure 9. The surface of a specimen 
prestrained in creep at 600 °F (fig. 9B) before creep 
testing at 1000 °F with a stress of 30,000 psi was 
approximately the same as that of a specimen com- 
pletely tested at 1000 °F with the same stress (fig. 9A). 
Extensive cracking was observed on the surface of 
the specimen pretested at 1200 °F (fig. 9C). However, 
the ductility of this specimen was no more seriously 
affected by these cracks than the one prestrained at 
600 °F. Surface cracks were also observed on a 
specimen creep-tested for a long time (fig. 9D) at 
1000 °F even though the reduction of area was greater 
for this specimen than for the one tested for a rela- 
tively short time (fig. 9A), Apparently, a different 
deformation mechanism contributed to increase duc- 
tility and counteracted the embrittling effects of 
surface cracking. This observation is also evident 
when the surface conditions of two specimens tested 
at 1200 °F for a short time (fig. 9E) and a long time 
(fig. 9F) are compared. 

The number, size and shape of cracks near the 
main fracture was affected by temperature, prior 
strain history, and rupture time as shown in figure 10. 
Elliptical holes, associated with low temperature 
deformation, are evident in figure 10A. As the tem- 
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FIGURE 9. Photographs of specimens after rupture in creep. X 5. 
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perature is raised, fewer internal cracks and holes 
were observed (fig. 10C) and even fewer were evident 
as the test time was increased (fig. 10D). Additionally, 
the flow and fracture mechanisms associated with 
prestraining at 600 °F (fig. 10B) contributed to the mi- 
crocracking and decreased ductility of this specimen. 

As shown previously (3), microstructure of the 
duplex-annealed alloy indicated the presence of equi- 
axed grains and lamellae of the beta phase in the 
alpha matrix. These microstructures were still evident 
in specimens tested to rupture at 600 °F (fig. 11A and 
B) but they were more evident in regions away from 
complete fracture (fig. 11B) than in regions near 
fracture (11A). As the test temperature was raised to 
1000 °F, less lamellar structure was evident as shown 
in figures 11C and D. The presence of beta is associated 


both with a strengthening effect and a decrease in 
ductility of the specimens. Moreover, cracking gen- 
erally occurred in the regions of the alpha grain bound- 
aries and was never observed in the regions having a 
lamellar structure. 

As shown in figure 12, cracks at alpha grain bound- 
aries, less distinct grain boundaries, and appreciable 
lessening of beta is more evident as the test tempera- 
ture was increased to 1200 °F. The increase in elonga- 
tion from 90 to 184 percent may be attributed to both 
recovery and the dissolution of beta into alpha during 
creep. 

Effect of prior strain history on microstructure is 
shown in figure 13. Grain boundaries are less distinct 
and the amount of intercrystalline cracking is de- 
creased as the prestraining temperature was increased. 
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FIGURE 10. Microstructures of specimens near axis near fracture. 
x 100. 


Longitudinal sections etched in a solution containing 3 ml HF, 6 ml HNO; and 100 ml H,0. 
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FiGuRE 1l. Effect of test temperature on microstructures of speci- 
mens after rupture in creep. X 500. 


Longitudinal sections, etched in a solution containing 3 ml HF, 6 ml HNOs, and 100 ml H,0. 
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Near fracture 
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These two factors are consistent indices of the rupture 
life at this temperature and stress, as it was previously 
shown (fig. 8) that the rupture life increased with 
increase in prestrain temperature. 


4. Summary 


Creep tests were made on dupiex-annealed Ti-8Al- 
1Mo-1V alloy specimens at temperatures of 600, 800, 
1000, and 1200 °F. 

Creep stresses were used to produce second stage 
creep rates ranging from 0.1 to 50,000 percent per 
1,000 hr and rupture times ranging from 0.1 to about 
10,000 hr. 

Conformance of the data to theoretical concepts 
is limited to small ranges of temperature and stress. 
This lack of agreement is attributed to the type of 
deformation, recovery, extent of cracking, phase trans- 
formation and prior thermal-mechanical history. 

Ductility at fracture generally increased with in- 
crease in temperature but varied inconsistently with 
test time. 

Rupture time increased with increase in prestrain- 
ing temperature and, at constant temperature, with 
decrease in the amount of prestraining. 

More cracks at the surface than in the interior were 


observed in specimens tested at high temperatures 
and slow creep rates, whereas there were fewer cracks 
in the interior of the specimens ruptured at the high 
temperatures. 

Increased ductility and low strength were associated 
with the dissolution of beta particles into the alpha at 
the high temperatures. 


The authors acknowledge the general assistance 
given by D. E. Harne throughout the investigation 
and the assistance given by J. Mackie with the 
metallography. 
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Comparison of incompressible flow and isothermal com- 
pressible flow formulae, J. Hord, Tech. Note 356 (Aug. 17, 1967), 
25 cents. 

Mass flow formulae for incompressible and “‘modified-incompress- 
ible” flow are compared with the isothermal compressible flow 
relation under the following conditions: The gas flow is steady, 
isothermal, and fully developed in a horizontal pipe of constant 
cross section with a prescribed static pressure drop (P;— P2). The 
comparative data are limited to static pressure ratios (P2/P,) > }, and 
subsonic isothermal flow. Laminar and turbulent flows are treated. 
Under the limitations of the comparison, modified-incompressible 
flow and isothermal gas flow relations are identical when 


fL(2D > ln (P;/P2). 


Graphical plots indicate the ee of approximation or error 
involved in using incompressible relations to solve comgnensne flow 
problems. Pressure losses due to end effects are briefly discussed. 


Key Words: Compressible flow, flow comparison, fluid flow, incom- 
pressible flow, mass: flow, pressure drop. 


Correlations for predicting leakage through closed valves, 
J. Hord, Tech. Note 355 (Aug. 1, 1967), 20 cents. 

Safety, convenience, and economy often demand the inference of 
leakage of a given fluid from known leakage of another test fluid 
across valve seats, welded and threaded fittings, seals, etc. The 
temperatures of the fluids also may be different. An example is the 
prediction of gaseous hydrogen leakage at 20 °K from test data 
obtained with nitrogen gas at 77 °K. Various flow formulae (molecu- 
lar, transition, and continuum) are examined, and two simple 
methods of correlating leakage for single-phase fluids are deduced. 
The correlations obtained (excluding transition fiow through long 
channels) indicate the leakage is inversely proportional to the square 
root of the density, or inversely proportional to the absolute viscosity, 
of the fluid. Thus, for gases, the leakage is directly proportional to the 
sonic velocity of the gas. The proper relationship must be estab- 
lished by experiment for each valve, fitting, etc. 


Key Words: Leakage through valves, leak rates, leak detection, 
predicting leak rates, valve seat leakage, leak flow rates, seal 
leakage. 


Designs for surveillance of the volt maintained by a small 
group of saturated standard cells, W. G. Eicke and J. M. 
Cameron, Tech. Note 430 (Oct. 9, 1967), 20 cents. 

This technical note describes a procedure for maintaining surveil- 
lance over a small group of saturated standard cells. The measure- 
ment process is briefly discussed and the principle of left-right 
balance as a means of eliminating certain systematic errors is 
developed. Specific designs and their analysis for intercomparing 
3, 4, 5 and 6 cells in a single temperature controlled environment are 
given. Procedures for setting up control charts on the appropriate 
parameters are given, and a technique is described for detecting 
certain types of systematic errors. 


Key Words: Control charts, experiment design, saturated standard 
cells, standard cells calibration, statistics, voltage standard. 


Realistic uncertainties and the mass measurement process. 
An illustrated review, P. E. Pontius and J. M. Cameron, Mono. 103 
(Aug. 15, 1967), 20 cents 

This paper gives a review of the concepts and operations involved in 
measuring the mass of an object. The importance of viewing measure- 
ment as a production process is emphasized and methods of eval- 
uating process parameters are presented. The use of one of the 
laboratory’s standards as an additional unknown in routine calibra- 
tion provides an accuracy check and, as time goes on, the basis for 
precision and accuracy statements. 


Key Words: Measurement, measurement process, uncertainty, 
mass measurement, precision, accuracy, statistical control. 


Standard Reference Materials: Catalog and price list of 
standard materials issued by the National Bureau of Stand- 
ards, Misc. Publ. 260, 1967 Edition (Sept. 15, 1967), 45 cents. 
(Supersedes NBS Misc. Publ. 260, 1965 Edition.) 

This catalog describes the various Standard Reference Materials 
issued by the National Bureau of Standards. These materials are 
used to calibrate measurement systems and provide scientific 
information that can be referred to a common base. A schedule of 
prices and quantities is included for each material, as well as direc- 
tions for ordering. Listed are the types and compositions of those 
chemical standards that are presently available. Announcements 
of new standard reference materials are made in ithe Federal Register, 
in scientific and trade journals, and in the Technical News Bulletin 
of the National Bureau of Standards. Changes affecting the current 
status of the various standards will be indicated by an insert sheet 
available quarterly from the Bureau. 


Key Words: Analysis, catalog, certificate, characterization, composi- 
tion, price list, property of material, purity of material, standard 
reference materials, standards. 


Standard x-ray diffraction powder patterns, H. E. Swanson, 
H. F. McMurdie, M. C. Morris, and E. H. Evans, Mono. 25, Sect. 5, 
(Aug. 31, 1967), 55 cents. 

Standard x-ray diffraction powder patterns are presented for 80 
substances. Thirty-nine of these patterns represent experimental 
data and 41 are calculated. The experimental x-ray powder diffrac- 
tion patterns are made with a Geiger counter x-ray diffractometer, 
using samples of high purity. All d-values were assigned Miller 
indices determined by comparison with theoretical interplanar 
spacings and from consideration of space group extinctions. The 
densities and lattice constants were calculated, and the refractive 
indices were measured whenever possible. The calculated x-ray 
powder diffraction patterns were obtained from published crystal 
structure data. The reported peak height intensities for calculated 
patterns were converted from integrated intensities. 

Reference intensity values based upon the strongest line of corun- 
dum (113) in a 50 percent weight mixture are given for 73 materials. 


Key Words: standard, x-ray diffraction, powder-patterns, crystal, 
structure, measurements, lattice, constants, reference-intensities. 


X-ray wavelength conversion tables and graphs for qualita- 
tive electron probe microanalysis, K. F. J. Heinrich and M. A. M. 
Giles, Tech. Note 406, (Sept. 25, 1967), 70 cents. 
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Tables and graphs have been constructed for electron probe x-ray 
spectromeiers equipped with crystal changers and calibrated for LiF 
crystals in angstroms. These tables indicate the equivalent reading 
for the lines which can be normally expected to appear in wavelength 
scans using the following crystals: LiF, EDDT, ADP, KAP, and Pb 
stearate (LSD). The readings were extracted from the tables of 
E. W. White, et al. The tables are arranged by the atomic numbers 
of the emitting elements; the graphs show the readings for observable 
lines as a function of atomic number for each of the aforementioned 
crystals. A suggested method for using the tables and graphs is 
included. 


Key Words: X-ray spectrometers, wavelengths, electron probe 
microanalyzer, x-ray emission. 


Field burnout tests of apartment dwelling units, D. Gross, 
Bldg. Sci. Series 10 (Sept. 29, £967), 25 cents. 

Results are reported of three burnout tests in an experimental test 
building, using a wood crib fuel load of 6 lb/ft?, representing combus- 
tible contents, and a structural design load of 40 lb/ft? applied to the 
floor or roof above the test room. Measurements were made of tem- 
perature, radiation, smoke, gas composition, and structural deflec- 
tion. A discussion of the fire performance of materials and methods 
of construction, and conclusions with regard to specific fire-protec- 
tive objects are presented. 


Key Words: Burnout test, fire performance, apartment dwelling, 
smoke, flame penetration, structural load, potential heat, fire load. 


“Life cycling” test on several strain gage pressure trans- 
ducers, P. S. Lederer, Tech. Note 434 (Oct. 1967), 20 cents. 

This publication reports the effects of many thousands of pressure 
cycles on the performance characteristics of several types of strain 
gage pressure transducers. The results obtained indicate permanent 
changes in zero pressure output and sensitivity. Most of these 
changes tend to occur during the first few thousands of cycles. The 
equipment and procedures used to obtain the data are described. 
Key Words: Life cycling, pressure cycling, life testing, pressure 
transducer, strain gage. 


Organic Chemistry; Air pollution studies characterization of 
chemical structures; synthesis of research materials; novel 
research materials; isotopic methods for analysis of carbo- 
hydrates; occurrence, preparation, and properties of 
naturally occurring monosaccharides (including 6-deoxy 
sugars); standard reference materials (organic) July 1966 
through June 1967, Ed. R. Schaffer, Tech. Note 427 (Oct. 1967), 
45 cents. 

This report describes work in progress in the Organic Chemistry 
Section of the Analytical Chemistry Division of the NBS Institute 
for Materials Research. The following research areas are discussed: 
Air Pollution Studies. Information has been obtained on the detec- 
tion of impurities in, and the purification of, commercial samples of 
certain polycyclic, aromatic hydrocarbons that are known to occur 
in polluted air. The oxidation of these hydrocarbons under conditions 
resembling those encountered in polluted areas has been studied, 
and some of the products have been identified. 

Characterization of Chemical Structures. The infrared absorption 
spectra of a group of 2-oxo-1,3-bis(phenylhydrazono) derivatives have 
been studied, and compared with the spectra of four tris(phenyl- 
hydrazones) and six related mono- and _bis-(phenylhydrazono) 
derivatives, in order to identify the important features of the bands 
arising from N=N, N-H (bending), and Ph-N groups. In addition, 
the phenylhydrazono—phenylazo tautomerism of 2-oxo-1,3-bis 
(phenylhydrazono) compounds has been elucidated. A novel aroma- 
tization of inositols in methyl sulfoxide—acetic anhydride has been 
discovered. 

Synthesis of Research Materials: A Disaccharide of D-Lyxofuranose. 
A compound formed in the high-vacuum distillation. of 2,3-0-iso- 
propylidene-a-D-lyxofuranose has been obtained crystalline by 
acetylation, crystallization of the diacetate, and saponification of 
the latter. From the results of oxidation, reduction, hydrolysis, and 
nuclear magnetic resonance studies, the compound is believed to be 
the di-o-isopropylidene acetal of 5-0-a-D-lyxofuranosyl-a-D-lyxofur- 
anose. This is the first known example of a thermal condensation 
providing a reducing disaccharide in which both sugar moieties are 
furanoid. 


Novel Research Materials. The reaction of sulfonic esters of alditols 
with methoxide has been examined, and a new dianhydrohexitol has 
been isolated. The literature on the infrared spectroscopy of carbo- 
hydrates has been critically evaluated and systematized. 

Isotopic Methods for Analysis of Carbohydrates. The various methods 
of isotope-dilution analysis, the isotope effects that may accompany 
studies with labeled molecules, and the methods for determining 
isotopic distribution have been assessed and assembled. 
Occurrence, Preparation, and Properties of Naturally Occurring 
Monosaccharides, Including 6-Deoxy Sugars. A survey has been 
made of all of the monosaccharides and 6-deoxy sugars that are 
known to occur naturally, and their preparation has been reviewed. 
In addition, their chemical and physical properties have been 
collated. 

Standard Reference Materials. The preparation of highly pure sam- 
ples of quebrachitol, Jevo-inositol, 1,2-0-isopropylidene-B-L-ido- 
furanose, and 2,3-0-isopropylidene-8-D-threo-pentulofuranose has 
been achieved, and the compounds have been characterized by 
melting point, specific optical rotation, and infrared absorption 
spectrum. These compounds are now available to the public as 
standard reference materials. 

In order to describe experimental procedures adequately, it has 
occasionally been necessary to identify commercial materials and 
equipment in this report. In no case does such identification imply 
recommendation or endorsement by the National Bureau of Stand- 
ards, nor does it imply that the material or equipment identified is 
necessarily the best available for the purpose. 


Key Words: Air pollution studies, oxidation of polycyclic aromatic 
hydrocarbons, characterization of chemical structures, phenyl- 
hydrazono—phenylazo tautomerism, disaccharide of D-lyxofuranose, 
infrared absorption methods for analysis of carbohydrates, naturally 
occurring monosaccharides and 6-deoxy sugars, standard reference 
materials, quebrachitol, Jevo-inositol, 1,2-0-isopropylidene-B-L- 
idofuranose, 2,3-0-isopropylidene-B-D-threo-pentulofuranose. 


Radiometric methods for measuring laser output,.D. A. 
McSparron, C. A. Douglas, and H. L. Badger, Tech. Note 418 (Nov. 
1967), 15 cents. 

The output of pulsed lasers may be measured with reasonable 
accuracy by modifications of established radiometric methods. The 
receiver, thermopile or phototube, is used ballistically. It is cali- 
brated by an incandescent-lamp standard of spectral irradiance. 
The laser beam is attenuated to make its irradiation on the receiver 
comparable to that from the calibrating source. Attenuation is 
accomplished by intercepting the entire laser beam on a diffusely 
reflecting surface. The attenuation factor is determined from 
direct measurement of the reflectance of this diffusing surface and 
the geometric constants of the system. The accuracy of the measure- 
ment of laser output is dependent upon the accuracy of calibration 
of the standard of spectral irradiance and the accuracy with which the 
spectral sensitivity of the receiver, the spectral transmittance of 
the band-pass filters, and the reflectance of the diffuser can be 
determined. Measurements made with different photoelectric re- 
ceivers agree to about 1.5% using a 1/4-joule pulsed ruby laser as a 
source. Total uncertainty is estimated to be about 5%. Measurements 
made using these radiometric methods have been compared with 
calorimetric measurements and a discrepancy of 9% was observed. 


Key Words: Laser. radiometry, laser energy measurement, radio- 
metric calibration, photoelectric photometry, optical attenuator. 


Second breakdown in semiconductor devices—A bibliog- 
raphy, H. A. Schaff, Tech. Note 431 (Oct. 1967), 30 cents. 

Almost 200 references with appropriate key words are listed which 
provide, it is believed, a comprehensive coverage of the literature of 
second breakdown in transistors and other semiconductor devices 
from 1958 through much of 1967. A representative list of earlier 
papers dealing with what appears to be second breakdown in point- 
contact and p-n junction diodes is also included. The indexes 
consist of an author index and an index to subject matter with 
reference tabulations and with key word assignments. 


Key Words: Bibliography, semiconductor devices, transistors, 
diodes, second breakdown, thermal breakdown, device reliability, 
failure modes, failure mechanisms. 
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X-ray wavelengths, J. A. Bearden, NSRDS-NBS 14 (Sept. 25, 
1967), 40 cents. 

Inconsistencies in accepted values (in x units) of x-ray reference lines 
have recently been demonstrated, although all are supposedly based 
on “good” calcite crystals. Factors supporting the selection of the 
W Ka. line as the X-Ray Wavelength Standard are critically dis- 
cussed. A review is given of the experimental measurements which 
are used to establish the wavelength of this line on an absolute 
angstrom basis. Its value is A W Ka; = (0.2090100+5 ppm) A. This 
may be used to define a new unit, denoted by A*, such that the 
W Ka, wavelength is exactly 0.2090100 A*; hence IA*= 1A+5 ppm. 
The wavelengths of the Ag Ka:, Mo Ka, Cu Kay, and the Cr Kaz 
have been established as secondary standards with probable error 
of approximately one part per million. Sixty-one additional x-ray lines 
have been used as reference values in a comprehensive review and 
reevaluation of more than 2700 emission and absorption wavelengths. 
The recommended wavelength values are listed in A* units together 
with probable errors; corresponding energies are given in keV. 
A second table lists the wavelengths in numerical order, and likewise 
includes their energies in keV. 


Key Words: Absorption, atomic, data compilation, emission, energy 
level, x-ray wavelength. 


Reevaluation of x-ray atomic energy levels, J. A. Bearden and 
A. F. Burr, NSRDS—NBS 14 (Sept. 25, 1967), 40 cents. 

All of the x-ray emission wavelengths have recently been reevaluated 
and placed on a consistent A* scale. For most elements these data 
give a highly overdetermined set of equations for energy level 
differences, which have been solved by least-squares adjustment for 
each case. This procedure makes “best” use of all x-ray wavelength 
data, and also permits calculation of the probable error for each 
energy difference. Photoelectron measurements of absolute energy 
levels are more precise than x-ray absorption edge data. These have 
been used to establish the absolute scale for eighty-one elements 
and, in many cases, to provide additional energy level difference 
data. The x-ray absorption wavelengths were used for eight elements 
and ionization measurements for two; the remaining five were interp- 
olated by a Moseley diagram involving the output values of energy 
levels from adjacent elements. Probable errors are listed on an 
absolute energy basis. In the original source of the present data, a 
table of energy levels in Rydberg units is given. Difference tables 
in volts, Rydbergs, and milli-A* wavelength units, with the respective 
probable errors, are also included there. 


Key Words: Absorption, atomic, data compilation, emission, energy 
levels, x-ray wavelength. 


A method of producing sturdy specimens of pressed powders 
for use in x-ray spectrochemical analysis, B. L. Bean, Appl. 
Spectry. 20, No. 3, 191-193 (May—June 1966). 

Sturdy specimens of compacted or pressed powders may be pro- 
duced for x-ray spectrochemical analysis as follows: A sleeve with 
sides about 0.16 cm (7/16 inch) thick and diameter 0.02 cm less than 
the inside diameter of the mold is placed in the mold. The powder 
is spread on the lower plunger of the mold. The sleeve is then re- 
moved and boric acid poured in the formed peripheral trench and 
on top of the powder. The specimen is then pressed. A specimen with 
strong edges and backing is prepared that may be reused many times 
because the boric acid edges support the pressure of the sample 
holder instead of the pressed sample powder. 


Key Words: Boric acid edges, pressed powders, sturdy specimens, 
x-ray spectrochemical analysis. 


Aromatic fluorocarbon polymers, W. J. Pummer and J. M. 
Antonucci (Proc. 152d American Chemical Society Meeting, New 
York, N.Y., Sept. 13, 1966), Polymer Preprint 7, No. 2, 1071-1076 
(Sept. 1966). 

Recent synthetic work has led to the preparation of a variety of 
fluoroaryl polymers. Most of the polymers synthesized are similar in 
structure to the hydrocarbon analogs except for the fluorine atoms. 
On this basis, the fluorinated polymers, presented in this paper, will 
be classified and discussed according to structural features. The first 
class of polymers contains the polyperfluorophenylenes which will 
include both the para and meta linked polymers. In the second class, 
the polymers contain an atom or group of atoms between the fluori- 
nated rings, such as oxygen, sulfur and oxydifluoromethylene groups. 


Some polymers in this category are the polyperfluorophenylene 
ethers and sulfides. The final group of polymers contains a pendant 
pentafluorophenyl ring attached to the main polymer chain such as 
the various fluorinated styrenes and vinyl ethers. 

In all cases, the method of polymerization, physical properties and 
thermal stability of the polymers will be discussed and whenever 
possible, comparisons will be made between the fluorinated and 
hydrocarbon polymers. 


Key Words: Fluorinated styrenes; fluoroaryl polymers; para, meta 
linked polymers; pendant pentafluorophenyl ring; polyperfluoro- 
phenylene ethers, sulfides; polyperfluorophenylenes; vinyl ethers. 


Calculation of currents of local galvanic cells, J. A. Simmons, 
S. R. Corviell, and F. Ogburn, J. Electrochem. Soc. 114, No. 8, 
782-787 (Aug. 1967). 

The current flow for a local electrolyte cell with rectangular anode 
and cathode is calculated as a function of polarization parameter, 
electrolyte thickness and conductivity, anode and cathode size, and 
zero current potential difference. The total current is obtained from 
the numerical computation (with error analysis) of a series solution of 
Laplace’s equation assuming a linear polarization relation at the 
electrode-electrolyte interface. For high values of the polarization 
parameter and sufficiently thick electrolytes an asymptotic formula 
obtains for the current approximately proportional to the anodic 
fraction—as was found by Waber et al. However, in general the cell 
geometry greatly influences the amount of current flow, especially 
for anode shapes near to linear rather than square. This effect can 
be greater than an order of magnitude. The calculated results are 
compared with experimental data obtained by Ogburn and Schlissel 
and are shown to be in good agreement. Corrosion currents in ran- 
domly pitted surfaces may be estimated from the above results. 


Key Words: Corrosion currents, galvanic cells, galvanic currents, 
Laplace’s equation, pitting corrosion. 


Crystallization of natural rubber, N. Bekkedahl, Rubber Chem. 
Technol. 40, No. 3, xxv—xlvii (June 1967). 

This review article presents a resume of the research work that has 
been done at the National Bureau of Standards during the past 
35 years in the field of crystallization of polymeric materials, with 
special emphasis on natural rubber. It contains only information on 
results of research work that has already been published. 


Key Words: Polymers, rubber, crystallization, melting, dilatometry, 
stark rubber, chain folding, heat engine, contractility. 


Determination of bound styrene in insoluble emulsion 
polymerised styrene-butadiene copolymers, M. A. Post, 
J. Appl. Chem. 17, 203-208 (July 1967). 

Bound styrene is determined in insoluble emulsion polymerized 
styrene-butadiene copolymers by measurement: of the 10.3 sm/ 
13.2 wm absorbance ratio using the baseline method. This ratio is 
obtained from the infrared spectrogram of the copolymer and the 
percent styrene is determined from a standard curve. The standard 
curve is constructed from absorbance ratios of 10.3 wm/13.2 wm 
derived from the infrared spectrograms of cast films of styrene- 
butadiene latexes over a range of percentage compositions. Bound 
styrene is determined in these latexes by nitration of the isolated and 
extracted copolymers. Styrene content based on nitration is deter- 
mined from standard curves relating styrene concentration to 
spectrophotometric absorbance measurements at wavelengths of 
285.0 nm, 273.8 nm and 265.0 nm. The source of the styrene for these 
curves is NBS standard rubber #1500. 


Key Words: Absorbance ratio method, bound styrene, infrared, 
insoluble styrene-butadiene copolymers, masonry paints. 


Electrical conductivity of high purity copper, J. J. Gniewek, 
J. C. Moulder, and R. H. Kropschot, Proc. LT—10 Conf., Moscow, 
USSR, Aug. 31-Sept. 9, 1966. 

The low temperature electrical conductivity of copper single crys- 
tals (Re73/R4> 30,000) and polycrystalline wire (Re73/R4> 10,000) 
has been measured between 4 °K and 77 °K. 

Between 14 and 30 °K the resistance is proportional to 7" where 
n=4.8+0.1. Size effects are observed using both eddy current 
and direct current techniques in samples as large as 2mm in diam- 
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eter. Assuming diffuse reflection of the electron from the surface 
direct current size effect measurements give a 


Poly = 0.53 X 10-"'Ocm? 
which is less than the free electron value of 0.66 X 10-"Qcm?. 


Key Words: Copper, electrical conductivity, low temperature, size 
effect. 


Epitaxial growth of iron on tungsten field emission points, 
A. J. Melmed, Surface Sci. 7, No. 3, 478-481 (July 1967). 

Iron crystals can be grown from the vapor phasé on tungsten field 
emission points. The crystals tend to nucleate at edges of W (011) 
planes and grow in the same orientation as the substrate, some- 
times with an apparent small parallel displacement. Inter-diffusion 
of W into the growing Fe crystals can occur and evidence for this 
is presented. 


Key Words: Epitaxy, field emission microscopy, crystal growth, 
vapor deposition, Fe-W alloying. 


Growth of large sodium chloride crystals from solution for 
color center studies, P. M. Gruzensky, J. Chem. Phys. 43, No. 11, 
3807 (Dec. 1, 1965). 

Sodium chloride single crystals with edge dimensions as large as 
2.5 cm were grown from aqueous solutions. Solutions were puri- 
fied by chemical precipitation of trace contaminants and _ filtra- 
tion through membrane filters. Crystal defects were apparent in 
in the initial stages of growth, but the outer portions of crystals 
were free of dislocations. F-band coloration by x rays at room 
temperature was not linear and the “first stage” was either absent 
or nearly so. Transition from extrinsic to intrinsic de conductivity 
occurred at 350 C with the extrinsic conductivity two orders of magni- 
tude lower than for “pure” melt-grown crystals. 


Key Words: Alkali halides, aqueous, color centers, chloride, con- 
ductivity, crystal growth, and sodium. 


Helium field-ion microscopy of hexagonal close-packed 
metals, A. J. Melmed (Proc. Symp. The Structure of Surfaces, 
Durham, N. C., Nov. 1966), Surface Sci. 8, No. 12, 191-205 (July- 
Aug. 1967). 

The atomic structure of curved hexagonal close-packed metal sur- 
faces is complicated compared to curved surfaces of body-centered 
cubic metals. He field-ion microscopy of Ru and Re field-evapo- 
rated surfaces discloses apparent differences in the surface struc- 
tures developed under similar experimental conditions. The effect 
of temperature and field strength on the structure developed by field- 
evaporation is analyzed, and the structures of the (100) and (001) plane 
regions in particular are discussed in relation to their respective 
field-ion images. 


Key Words: Field-ion microscopy, field evaporation, surface struc- 
ture and hexagonal close-packed metals. 


High pressure adsorption isotherms of neon, hydrogen, and 
helium at 76 °K, A. J. Kidnay and M. J. Hiza (Proc. 1966 Conf., 
Boulder, Colo., June 13-15, 1966), Book, Advances in Cryogenic 
Engineering 12, 730-740 (Plenum Press Inc., New York, N.Y., 1967). 
A desorption technique was used to measure the physical adsorption 
isotherms of hydrogen and helium on a synthetic zeolite and hydro- 
gen, helium, and neon on an activated cocoanut shell charcoal at 
76 °K in the pressure range of 1 to 95 atmospheres. 

The Polanyi adsorption potential theory, as modified by Dubinin, 
was used in an attempt to obtain a general correlation of these iso- 
therms. The results, while not completely general, should allow the 
prediction of the high pressure isotherms of neon, hydrogen, and 
helium at any temperature with an accuracy sufficient for most 
engineering calculations. 


Key Words: Adsorption potential, charcoal, helium, high-pressure 
physical adsorption, hydrogen, low temperature, neon, zeolite. 


History of the National Standard of electromotive force, 
W. J. Hamer (Proc. 22nd ISA Conf. and Exhibit, Chicago, IIl., 


Sept. 11-14, 1967), ISA Preprint No. M2-1-MESTIND-67, pp. 
1-10, (1967). 

The history of the National Standard of electromotive force is traced 
from the meeting of the International Electrical Congress in Chicago 
in 1893 to the present time. Included are discussions of the transfer 
of the “volt standard” from its old location in Washington, D.C., 
to its new location near Gaithersburg, Md., of the stability of stand- 
ard cells, of the application of Gibbs’ Phase Rule to standard cells, 
and the relation between the legal and absolute volts. 


Key Words: Volt standards, legal volt, absolute volt, Gibbs’ phase 
rule, history of volt. 


Indentation hardness of the fingernail, S. B. Newman and 
R. W. Young, J. Invest. Dermatol. 49, No. 1, 103-105 (Mar. 17, 1967). 
Indentation hardness using the Knoop indenter has been used as a 
quantitative criterion in physiological studies of the human nail 
plate. Data were collected and are presented to show that values 
reported in the literature are in serious error. A number of param- 
eters affecting Knoop numbers are examined in an effort to deter- 
mine the source of error. 


Key Words: Fingernail, Knoop indenter, hardness. 


Polymerization of fluoroolefins and related monomers, 
L. A. Wall (Proc. 152nd American Chemical Society Meeting, New 
York, N.Y., Sept. 13, 1966), Polymer Preprint 7, No. 2, 1112-1115 
(Sept. 1966). 

Although some fluorine-containing monomers polymerize readily, 
there are many which have not as yet been polymerized either by 
free radical or ionic catalysts. In order to investigate the reasons 
for the somewhat general tendency of these monomers to resist 
polymerization, we have initiated a systematic study using high 
pressures and y-rays to induce polymer formation. 

The technique has been successful for many fluoro monomers, 
fluoropropylenes, fluorodienes, perfluorostyrene and some fluoro- 
vinyl phenyl ethers. In a general way, the results suggest in many 
cases that slow propagation steps, easy transfer processes with 
monomer or impurities, and in some cases a tendency to thermally 
dimerize, interfere with the formation of high polymer. 


Key Words: Fluorine-containing monomers, fluorodienes, fluoro 
monomers, fluoropropylenes, fluorovinyl, free radical, gamma-rays, 
high pressures, ionic catalysts, perfluorostyrene, phenyl ethers. 


A computer simulation model of railroad freight transporta- 
tion systems, W. P. Allman, Jntern. Railway Cong. Assoc. Mo. 
Bull. Cybernetics and Electronics 4, No. 2, 45-57 (Feb. 1967). 
A railroad may be viewed as a network of nodes (yards) and links 
(railroad lines), with time-dependent demands for the movement of 
freight cars imposed upon it. Incident to car movement are the 
following key operating policy questions of railroad freight opera- 
tions: 
a) When and where should regularly-scheduled trains run? 
b) For each yard, what should the car classification (sorting) policy 
be? 
c) Over each link of a train’s route, what cars should be assigned 
to the train for hauling? 
These interdependent questions must be answered simultaneously, 
and in accordance with overall operating objectives of the railroad 
enterprise. Policies must be revised periodically when significant 
changes in demand traffic patterns occur. 
The paper describes a simulation model which permits experimen- 
tation with various alternative railroads freight operating policies 
at a total-network level. The model has been constructed with the 
SIMSCRIPT simulation programming language. Model inputs 
include time-dependent freight car origin-destination demand 
data, train routes and schedules, yard sorting and operation policies, 
and assignments of cars to trains for hauling. Freight cars are sorted 
at yards, and picked up and dropped off by trains which haul them 
thru the network. Model outputs include several railroad oper- 
ation performance measures such as origin-destination transit 
times, activity volumes, train lengths, delays incurred by cars 
at yard operations, and operating costs. 


Key Words: Railroad freight transportation systems, computer 


simulation model, 


transportation systems, network-simulation. 
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A precision noise-power comparator, M. G. Arthur, C. M. All- 
red, and M. K. Cannon, JEEE Trans. Instr. Meas. IM-13, No. 4, 
301-305 (Dec. 1964). 

This paper describes a prototype noise power comparator based 
upon a theory given by Allred. Operating at 3 MHz, it is a null- 
type instrument, the principal components of which are a reference 
CW voltage generator, a hybrid four-port, a dual-channel ampli- 
fier and band-pass filter, and an analog multiplier. Unlike other 
radiometers, no rapid switching of the noise power or reference 
voltage is performed. 

The instrument can compare noise powers having effective noise 
temperatures in the range from below liquid nitrogen tempera- 
ture to greater than 30,000 °K. Two noise generators having known 
spectral densities are used to calibrate the comparator. The accuracy 
of comparison is 1% at 77 °K and increases to 0.2% at 29,000 °K. 


Key Words: Analog multiplier, attenuator, band-pass filter, hybrid 
four-port, noise-power comparator, noise powers, noise tempera- 
tures, radiometer. 


Microprinting fills the bill, B. M. Fry, Business Graphics, 
pp. 21-25 (Jan. 1967). 

The Clearinghouse for Federal Scientific and Technical Informa- 
tion, part of NBS, is a focal point for the collection, announcement 
and sale of copies of unclassified U.S. Government R&D reports 
released for public sale by the sponsoring Federal agencies. Approx- 
imately 50,000 new reports and Government-sponsored transla- 
tions of foreign technical material are received each year. The 
article describes Clearinghouse reproduction services as part of 
the system to meet the demand for over 2,000,000 copies of docu- 
ments in microfiche and paper copy each year. The Clearinghouse 
reproduction plant makes use of the latest technology in micro- 
photography, electrostatic copying and offset printing. The article 
highlights the role of the reproduction plant as a laboratory in test- 
ing new equipment and ideas to improve service to customers. 


Key Words: Clearinghouse for Federal Scientific and Technical 
Information, microphoto, microprinting, printing, systems. 


Microwave standards and measurements in the USA, R. W. 
Beatty (Proc. XV General Assembly of URSI, Munich, Germany, 
Sept. 5-15, 1966), Book, Progress in Radio Science 1963-1966, Part 
1, pp. 194-260 (International Scientific Radio Union, Munich, 
Germany, 1966). 

A report is given on progress within the U.S.A. in the development 
of radio standards and measurement methods at frequencies from 
1 to 300 GHz (GHz= 10° cycles per second). The period covered is 
the three years 1963-1966. The following topics are included within 
the scope of this report: power; noise; reflection coefficient and 
impedance; phase shift; attenuation and loss; and field strength and 
antenna gain. 


Key Words: Antenna gain, attenuation, field strength, impedance, 
measurement, microwave, noise, phase shift, power, progress, 
reflection coefficient, review, standards, survey. 


Proposed recommended practice for description and selec- 
tion of conditions for photographing specimens, C. S. Mc- 
Camy, Photo. Sci. Eng. 10, No. 4, 185-191 (July—Aug. 1966). 

A “Proposed Recommended Practice for Description and Selection 
of Conditions for Photographing Specimens” has been prepared for 
the American Society for Testing and Materials so that the photo- 
graphic method of recording appearance can be made more repro- 
ducible from one laboratory to another, for various specimens, and 
for a given specimen at different times. The proposal defines terms 
and symbols and provides a systematic method of describing the 
arrangement of lights, camera, and subject; the characteristics of 
the illumination; the nature of the photographic process, and the 
viewing system. Conditions for photographing certain common 
forms of specimens are recommended. Conventions for exhibiting 
and publishing photographs of specimens include an orientation 
mark and a coordinate system for pictures. 


Key Words: Specimen photography, camera description, lighting 
description, appearance. 


Protection circuit insulated for high voltage, E. Niesen, Rev. 
Sci. Instr. 38, No. 5, 689 (May 1967). 


A circuit using standard components for floating an optical meter 
relay at high voltage is presented. By plugging in available optical 
meter relays, sensitivities in microamperes to amperes, millivolts 


-to kilovolts can be controlled for a high or low limit. 


Key Words: Optical meter relay, protection circuit. 


Size and power requirements of 4.2 °K refrigerators, T. R. 
Strobridge and D. E. Chelton (Proc. 1966 Conf., Cryogenic Engineer- 
ing, Boulder, Colo., June 13-15, 1966). Book, Advances in Cryogenic 
Engineering 12, 576-584 (Plenum Press Inc., New York, N.Y., 1967). 
A survey to determine the physical characteristics of 4.2 °K cyclic 
refrigetators is presented. Although the data acquired through the 
survey show considerable scatter, the expected trends can be de- 
tected and attempts have been made to establish guides by which 
the weight, volume, and power requirements for various capacity 
refrigerators may be estimated. No attempt is made to assess the 
reliability, the interval of maintenance free operation, or the capital 
costs. 


Key Words: 4.2 °K refrigerators, cyclic refrigerators, liquid helium, 
power requirements, size requirements. 


Slush hydrogen characteristics, D. B. Mann, C. F. Sindt, P. R. 
Ludtke, and D. B. Chelton, Proc. Conf. Long-Term Cryo-Propellant 
Storage in Space, National Aeronautics and Space Administration, 
Huntsville, Ala., Oct. 12-13, 1966. 

The utilization of slush hydrogen as an upgraded fuel is dependent 
upon knowledge of its physical and transport properties and char- 
acteristics. A NASA-sponsored program designed to acquire this 
information is in progress. 

Experimental measurements of solid particle size distribution, aging 
effects and terminal velocities of the solid particles in the liquid 
melt have been completed. Transport characteristics have been 
predicted using these measured parameters. These predictions are 
currently being submitted to experimental verification. 

Two 100 gallon dewars connected by 75 feet of vacuum-insulated 
transfer line comprise the experimental flow loop. Pressure sensing 
transducers and visual observation through glass ports allow ac- 
quisition of transport data over a broad range of flow rates. 
Handling and investigative techniques such as dewar-to-dewar 
transfer, straining of the solids through wire mesh screen, and 
vertical temperature profiles of the settled mixture are presented 
in the context of a continuing experimental and analytical research 
program. 


Key Words: Cryogenic, fluid transport, freezing, hydrogen-liquid 
solid, production. 


The role of NBS in motor vehicle satefy research, P. J. Brown, 
Proc. Regional Meeting Society Automotive Engineers, Gaithers- 
burg, Md., April 18, 1967, SAE 670204, p. 3 (June 1967). 

The relationship of the National Traffic Safety Agency and NBS 
and an outline of the research programs in tire systems, occupant 
restraint systems, and braking systems. The immediate goals and 
objectives as well as the long-range objectives will be stipulated. 


Key Words: Tire system research, occupant restraint system re- 
search, braking system, research information exchange, role of 
NBS motor vehicle safety. 


A modified Monte-Carlo quadrature. II, S. Haber, Math. 
Compt. 21, No. 99, 388-397 (July 1967). 

A modification of Simple Monte-Carlo quadrature is proposed, which 
uses very simple forms of stratified sampling and of the “method of 
antithetic variates.”” The new procedure is fully automatic —requir- 
ing no preliminary analysis of the integrand—and converges some- 
what faster than Simple Monte-Carlo. Results of experimental 
calculations are presented. 


Key Words: Mathematics, analysis, numerical analysis, integration, 
numerical integration, quadrature, Monte-Carlo, multiple integration. 


Current distribution and impedance of lossless conductor 
systems, R. L. Brooke and J. E. Cruz, JEEE Trans. Microwave 
Theory Tech. MTT-15, No. 6, 358-364 (June 1967). 
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A general method for determining the characteristic impedance of 
lossless, uniform, go and return systems is developed. The current 
distribution within the system is determined by means of a matrix 
equation programed for computer solution. Knowing the current 
distribution, the inductance per tnit length and characteristic 
impedance are obtained. The results obtained by applying this 
method to several rectangular coaxial systems are compared with 
the predictions of an approximate analytic expression. The reflection 
coefficient of a variable characteristic impedance coaxial line is 
measured on a “Time Domain Reflectometer” (TDR) and the results 
are'compared with both the matrix method and the approximate 
analytic expression. 


Key Words: Characteristic impedance, current distribution, general 
method, matrix solution, rectangular transmission line, time-domain 
reflectometer, variable impedance. 


The spheroidal method in satellite astronomy, J. P. Vinti, 
Book, Space Mathematics 5, Part 1, 119-129 (American Mathe- 
matical Society, Providence, R.1., 1966). . 

To calculate the drag-free orbit of an artificial satellite of an oblate 
planet, the author has devised a very accurate approximation for its 
gravitational potential, leading to separability of the problem in 
oblate spheroidal coordinates. The present paper sketches the 
procedures for finding this potential and for calculating the corre- 
sponding orbit. It then discusses the remaining perturbing potential 
and a method for accounting for its effect on the orbit. 


Key Words: Artificial satellite, drag-free orbit, gravitational potential, 
oblate planet, oblate spheroidal coordinates, satellite astronomy, 
spheroidal method. 


A systems concept of electromagnetic measurements in the 
U.S.A., G. E. Schafer, Proc. IEEE 55, No. 6, 775-778 (June 1967). 
Abstract is not required, but the following summary describes con- 
tents of manuscript. A brief description of a systems concept of the 
National Measurement System proposed by Dr. R. D. Huntoon, 
Director, Institute for Basic Standards, is given first. Activities of the 
Radio Standards Laboratory as a functional element of this system 
are then described. Some benefits of looking at RSL as an element 
of this system are given. 


Key Words: National measurement system, radio standards, elec- 
tromagnetic measurement. 


Liquid hydrogen engineering instrumentation, T. M. Flynn 
(Proc. Intern. Inst. Refrigeration, June 9-11, 1965), Book, Liquid 
Hydrogen, pp. 295-338 (Paris, France, 1965). 

Prudent liquid hydrogen engineering requires the measurement of 
both extensive and intensive properties ofthe cryogenic liquid. 
Transducers are required for liquid level (quantity), both point and 
continuous systems, and mass rate systems. In addition, there must 
be transducers of pressure, temperature, density, and occasionally, 
quality. This paper discusses some of the devices and practices 
currently used for the measurement of pressure,‘ temperature, flow 
rate, and liquid level in liquid hydrogen systems. 


Key Words: Cryogenics, instrumentation, liquid hydrogen, trans- 
ducers. 


Progress in the distribution of standard time and frequency, 
1963 through 1965, J. M. Richardson (Proc. XV General Assem- 
bly of URSI, Munich, Germany, Sept. 5-15, 1966 ), Book, Progress 
in Radio Science 1963-1966, Part 1, pp. 40-62 (International Sci- 
entific Radio Union, Munich, Germany, 1966). 

Progress in accurate long distance distribution of standard time 
and frequency, as reported in the literature from 1963 through 1965, 
is summarized. Techniques are by VLF, LF, and HF radio propaga- 
tion, by satellite relay, and by portable clocks. Effects on standard 
frequency transmissions of variations in VLF propagation with 
geophysical phenomena are quantitatively understood. VLF and 
LF transmissions have provided careful, long-term, statistical com- 
parison of remotely located atomic frequency standards. Precision 
of at least 2 parts in 10" for a 24-hour observation period is possible 
at 5000 km. The phase of some standard frequency transmitters 
is routinely steered by VLF from distances up to 5300 km. Global 
distribution of standard time by VLF to microsecond resolution has 
been shown feasible. The null beat between two neighboring VLF 


> 


carriers propagates stably enough to mark a particular VLF cycle, 
and the beat period can be long enough to enable ordinary time 
signals to mark a particular null beat. Intercontinental time syn- 
chronization by microwave pulses has been accomplished via Tel- 
star and Relay II satellites. Accuracy is observed to be at several 
microseconds. Portable cesium clocks have served as global trans- 
fer standards with degradation of timing accuracy of only about a 
microsecond per trip. Results by all the above methods are consist- 
ent with each other and with stated accuracies of atomic standards 
involved. 


Key Words: Time, frequency, time signals, standard frequency 
broadcasts, VLF propagation, portable clocks. 


Progress in the United States on electromagnetic standards 
and measurements at 30 kHz to 1 GHz, 1963 through 1965, 
M. C. Selby, JEEE Trans. Instr. Meas. IM-16, No. 2, 172-178 
(June 1967). 

A digest of highlights is presented on the most significant U.S. 
contributions to the measurement of attenuation, impedance, volt- 
age, current, field strength, thermal noise, and phase at 30 kHz to 
1 GHz. A total of approximately 30 contributions are digested in 
this section of the report. The following accomplishments are among 
them: a supersensitive detector for a complex-insertion-ratio meas- 
urement system having accuracies of about 0.0005 dB/10 dB at 
30 MHz; exact equations for mutual and self-inductance of various 
combinations of filaments, tapes, and bars; a modified Twin-T-Bridge 
for measuring resistances of 100 to 10,000 ohms to 15 MHz; a set 
of Q-factor standards for frequencies to 45 MHz based on accumu- 
lated data and experience over 5 years; a unique adjustable char- 
acteristic-impedance coaxial line measurement of Q’s greater than 
100,000 of cryogenic circuits at frequencies to 300 MHz; a novel 
Tee-junction to enable calibrations of voltmeters of any practicable 
input impedances with VSWR’s ranging from 1 to 200, to 1 GHz and 
higher; a miniaturized dipole-antenna field strength meter, employ- 
ing a semiconducting plastic transmission line, to measure complex 
near-zone fields of 0.1 to 1000 volts per meter, from 150 kHz to 30 
MHz; a prototype 3-MHz model of precision thermal noise-power 
comparators for an equivalent noise-temperature range of 75 to 
30,000 °K at accuracies of 0.2 to 1 percent. 


Key Words: High frequency measurements, progress in radio meas- 
urements, attenuation measurements, impedance measurements, 
voltage measurements, field strength measurements, thermal noise 
measurements, phase measurements. 


The conference on precision electr gnetic re- 

ments, H. M. Altschuler, Proc. IEEE 55, No. 6, 782 (June 1967). 

The character of the biennial Conference on Precision Electro- 

magnetic Measurements (CPEM) is described. The scope of the 

Conference is compared with that of the special issue of the Pro- 

ceedings of the IEEE on Measurements. Some of CPEM’s history 

is recounted and its increasing international coverage is noted. It 

is concluded that the challenge to CPEM grows with the expansion | 
of the fields it covers. 





Key Words: Conference, measurements, precision measurements, 
electromagnetic measurements, conference on precision electro- 
magnetic measurements. 


The nature of measurement, A. G. McNish, Book, Handbook of 
Industrial Metrology, Chapt. 1, pp. 2-12 (Prentice-Hall, Inc., Engle- 
wood Cliffs, N.J., 1967). 

General principles for measurement. Accuracy and precision are 
discussed. 


Key Words: Accuracy, measureme?t, precision. 


U.S.A. National Committee Report, Fifteenth URSI General 
Assembly, Munich, September 1°56: Commission 1, Radio 
Measurement Methods and Standards, Progress in Radio 
Measurement Methods and Standards, J. M. Richardson, Edi- 
tor, Radio Science 1, No. 11, 1333-1342 (Nov. 1966). 

Progress in radio measurement methods and standards within the 
United States during the Triennium 1963 through 1965 is reviewed 
for inclusion in the report of the U.S. National Committee of URSI 
to the 1966 General Assembly. Topics surveyed are high precision 
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atomic frequency standards, high precision quartz frequency 
standards, scientific aspects of Universal and/or Atomic time and 
frequency transmissions, standards and measurements at 30 kc/s 
to 1 Ge/s, standards and measurements at 1 to 300 Ge/s, precision 
coaxial connectors, swept-frequency techniques, measurements of 
electromagnetic properties of materials, and laser standards and 
measurements. The topics are treated within the length limitatio 

imposed by USNC/URSI. 


Key Words: Radio measurement standards, radio measurement 
methods, frequency standards, time standards, electromagnetic 
properties of materials, coaxial connectors, microwaves, lasers. 


Laser measurements of long distances, P. L. Bender, Proc. 
IEEE 55, No. 6, 1039-1045 (June 1967). 

The extension of radio frequency techniques to the optical region 
of the spectrum has made possible many types of measurements 
which were not feasible before. One area in which important scien- 
tific contributions can be expected during the next few years is 
the use of lasers to measure long distances with high accuracy. 
Three types of distance measurements which have been discussed 
in the literature and which are now being investigated actively are: 
a) interferometric measurements over distances of up to hundreds of 
meters through evacuated or sealed-off tubes; b) measurements 
with modulated laser beams over distances of perhaps one to 
fifty kilometers with corrections made for the atmospheric index 
of refraction along the path; and c) range measurements to arti- 
ficial satellites and to the moon using laser radar. Some of the pos- 
sible geophysical and geodetic applications of such measurements 
will also be discussed. 


Key Words: Laser measurements, radio frequency techniques, 
interferometric measurements, modulated laser beams, range 
measurements. 


Low temperature (295-4 K) mechanical properties of 
selected copper alloys, R. P. Reed and R. Mikesell, J. Mate- 
rials 2, No. 2, 370-392 (June 1967). 

Fifteen copper alloys, including brasses, bronzes, and commer- 
cially pure coppers were tested to determine their tensile, notch 
tensile, Youngs modulus, and impact properties. Tests were per- 
formed at 295, 195, 76, 20, and 4 °K. It was found that copper al- 
loys, in general, have improved room low temperature properties 
as compared to room ambient temperature properties. The excep- 
tion to this result is a cast Ni-Al bronze alloy, which is brittle at 
low temperatures. In 9 alloys the yield and tensile strengths at 
4 °K were less than at 20 °K. A discussion of this strength reduc- 
tion at low temperatures is presented. 


Key Words: Mechanical properties, low temperatures, copper, 
copper alloys. 


Meteorology, essential of space age, A. G. McNish, George 
Washington University Mag. 1, No. 4, 10-14 (1965). 
Understanding the physical world and the progress of civilization 
depends upon measurement. 


Key Words: Metrology, standards, units, education. 


Optical properties of germanium in the far ultraviolet, 
L. Marton and J. Toots, Phys. Rev. 160, No. 3, 602-606 (Aug. 15, 
1967). 

Reflectance measurements for Ge films in situ are reported for 
16 wavelengths between 490 and 1610 A, corresponding to photon 
energies of 25.2 and 7.7 eV. Measurements were taken after 3-min 
exposures to ambient pressures of 510-7 Torr. Reflectance 
values at a number of angles of incidence were used to compute 
the optical constants: the refraction index n, the extinction coefh- 
cient k, and the real and imaginary parts of the frequency-dependent 
dielectric constant €(w)=€:+i€2. Values are also presented for 
the characteristic electron energy-loss functions, i.e., the volume 
loss, —Ime-!, and the surface loss, —Im(e€+1)~-!. The quantity 
—Ime-' peaked at 16.1 eV and had a half-width of 3.9 eV. Electron 
oscillator strengths, summed between 8 and 25 eV, gave a hap 
value of 13.7 eV. The present data are in better agreement with 
the accepted plasma theory than are previous data. Some changes 
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of optical characteristics caused by exposure to ambient pressures 
and air are presented. 


Key Words: Characteristic electron energy loss, dielectric con- 
stant, far ultraviolet, germanium, optical constants, reflectance. 


Properties of liquid hydrogen, R. J. Corruccini (Proc. Intern. 
Inst. Refrigeration, June 911, 1965), Book, Liquid Hydrogen, 
pp 65-106 (Paris, France, 1965). 

A review is presented of the physics of liquid hydrogen, and exten- 
sive extracts are provided of the available data on those of its physi- 
cal properties that are important in technology. 


Key Words: Cryogenic, hydrogen, liquid, low-temperature, proper- 
ties. 


U.S. Department of Commerce and the flammability of 
clothing fabric, A. F. Robertson, Bull. N.Y. Acad. Med. Second 
Series 43, No. 8, 706-710 (Aug. 1967). 

The National Bureau of Standards has conducted studies of existing 
and modified test methods for measurement of flammable hazards 
associated with clothing fabrics. Recently, a cooperative program 
has been initiated with the Public Health Service to secure burn 
accident experience records which should yield information on 
the type of fabrics currently associated with clothing burn injuries. 
This program is just starting and it will be some years before valid 
information of the national picture becomes available. To date, 
however, there appears to be no technical evidence to show that a 
significant modification of the burn accident experience can be 
achieved without major changes in the flammable behavior of 


clothing fabrics. Thus, we hope industry will continue research on 


and development of permanently retardant-treated fabrics which 
will be economically, esthetically, and usefully acceptable to the 
public. It appears, in spite of the lack of direct evidence, that such 
fabrics will significantly reduce, but not eliminate, clothing fire 
injuries. 


Key Words: Burn accident experience records, clothing fabrics, 
flammability, measurement of flammable hazards, permanently 
retardant-treated fabrics. 


Perturbation theorems for waveguide junctions, with 
applications, D. M. Kerns and W. T. Grandy, Jr., IEEE Trans. 
Microwave Theory Tech. MTT-14, No. 2, 85-92 (Feb. 1966). 
Perturbation theorems are derived in the context of a theory of 
waveguide junctions. These theorems express changes in impedance 
or admittance matrix elements, due to changes in a waveguide 
junction, in terms of integrals over products of perturbed and 
unperturbed basis fields associated with the junction and with its 
adjoint. Media involved are required only to be linear. 
Concepts of first-order perturbation theory are discussed briefly, 
and the term “correct to the lowest order” is precisely defined. 
The need of explicit theorems telling when one may expect results 
actually correct to the lowest order is noted. 
Two problems are solved approximately by the perturbation ap- 
eer (1) Reflection at the junction of rectangular waveguide with 
lleted waveguide of the same main dimensions, and (2) the effect 
of finite conductivity of both obstacle and waveguide wall for half- 
round inductive obstacles in rectangular waveguide. 


Key Words: Correct to the lowest order, first-order perturbation 
theory, perturbation theorems, waveguide junctions. 


Progress in measurement of electromagnetic properties 
of material, 1963 to mid-1966, H. E. Bussey (Proc. XV General 
Assembly of URSI, Munich, Germany, Sept. 5-15, 1966), Book, 
Progress in Radio Science 1963-1966, Part 1, pp. 265-289 (Inter- 
national Scientific Radio Union, Munich, Germany, 1966). 
Progress on radio and microwave measurements of dielectric and 
magnetic properties of materials is reviewed mainly for the period 
1963 to mid-1966. Important materials properties, theoretical analysis 
of circuits containing materials, and representative measurement 
methods are reviewed. In addition, some earlier basic work is 
referenced. 


Key Words: Review, permittivity, dielectric constant, permeability, 
ferromagnetic resonance, bibliography. 
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Development of an adhesive restorative material, R. L. Bowen 
(Proc. 2nd Workship, Adhesive Restorative Dental Materials, 
University of Virginia, Charlottesville, Va., Dec. 8-9, 1965), Book, 
Adhesive Restorative Dental Materials, No. 2, Chapt. 4, p. 225 (U.S. 
Public Health Service Publ. 1494, Supt. Doc. Government Printing 
Office, Washington, D.C. 20402 $1.00, 1966). 

Recent findings, primarily from the literature, are reviewed on 
the subject of adhesive, restorative, dental materials, One mode of 
attack on the problem is described. In this, an inorganic-organic 
composite is used as the restorative material. The inorganic re- 
inforcing fillers serve primarily to reduce the material’s thermal 
expansion. A coupling agent improves the bonding between the 
phases. Bonding between the organic resin phase and the tooth 
surface is also improved by a coupling agent (surface-active 
comonomer). 


Key Words: Adhesive, bonding, dental material boneing, resin, 
restorative. f 


Identification of dyes in paper by extraction and chromato- 
graphic analysis, M. R. Wallace, L. T. Milliken, and S. D. Toner, 
TAPPI 50, No. 9, 121A-124A (Sept. 1967). 

A chromatographic procedure is described for the determination 
of a range of direct and acid dyes in dyed paper. It can be used for 
the identification of individual dyes in complex mixtures on unknown 
samples, and is simple and inexpensive enough to use in routine 
quality control testing of dyes and dyed paper. 

R; values and visual appearances serve immediately to identify 
most of the paper dyes studied. Comparisons with known references 
give good estimates of the quantities of each dye present. Individ- 
ual basic dyes could not be identified by this procedure, but the 
presence of a basic dye could be detected by its appearance in the 
solvent front. In a mixed furnish containing both acidic and basic 
dyes the acidic dyes could be determined without interference by 
the basic dyes. Ry values for the principal characteristic bands of 
26 dyes are listed. 


Key Words: Chemical analysis, chemical tests, chromatography, 
dyes. 


Modification of the fluoride activity electrode for micro- 
chemical analysis, R. A. Durst and J. K. Taylor, Anal. Chem. 39, 
No. 12, 1483-1485 (Oct. 1967). 

A fluoride activity electrode has been modified for operation in the 
inverted position by converting the inner reference solution to a 
gel. Operation in this position makes possible the use of the single- 
crystal fluoride-sensitive membrane as the sample “container” 
and permits the measurement of sample volumes as small as one 
drop (50 1). A calibration curve is obtained for solutions containing 
10-' M to 10-* M NaF in 0.1 M KNOs. The precision to which the 
fluoride concentration can be read from this calibration curve is 
calculated from the slope of the curve at the various concentrations 
and the standard deviation of the e.m.f. measurements. In the con- 
centration range from 10-! M to 10-3 M NaF, the calculated error 
is about 2%. At the 10-5 M NaF level, 10 nanograms of fluoride can 
be determined with an error of approximately 0.7 nanogram. For 
sample volumes of 50 41, the modified fluoride electrode was found 
to exhibit a usable response in the range from 10-! M to5 x 10-*° M 
fluoride. 

Key Words: Microchemical analysis; potentiometry; fluoride, 
analysis for; Specific Ion Electrode; fluoride electrode; electro- 
chemical analysis. 


Modified linear null-point potentiometry, R..A. Durst and J. K. 
Taylor, J. Chem. Educ. 39, No. 12, 1374-1379 (Oct. 1967). 

A modification of the technique of precision null-point potentiometry 
has been developed in which a linear titration curve is obtained 
when the logarithm of the amount of coulometrically generated 
titrant is plotted against the concentration cell potential. The 
euivalence point is evaluated graphically and by computer tech- 
niques. The analyses of standard 0.1 ml samples containing 0 to 13 
pg of silver are discussed in detail. This technique was found to 
be applicable down to one nanoequivalent (0.1 yg) of silver with 
an error of less than 5%. Below this value, a positive deviation 
occurs due to the spontaneous dissolution of the silver indicator 
electrodes. The interference effects of several cations are also 
discussed. 


Key Words: Null-point potentiometry, silver, microtitration, linear 
potentiometry, electroanalysis, potentiometry. 


Dead weight piston gauge for pressures to 26 kilobars, D. P. 
Johnson and P. L. M. Heydemann, Rev. Sci. Instr. 38, No. 9, 1294- 
1300 (Sept. 1967). 

The design and operation of a dead-weight piston gage for pressures 
to 26 kbar is described. The gage is of the controlled clearance type. 
The diameter of the gage piston is .2 cm. It is loaded directly with 
up to 840 kg of weights. A volume of about 1 cm? is provided inside 
the gage for the material under test. The use of this gage for transi- 
tion pressure determinations is briefly discussed. 


Key Words: High pressure, piston gage, transition pressure, pressure 
measurement. 


Lou Costrell of NBS speaks out on the revolution in instru- 
ment packaging, L. Costrell, EEE 15, No. 8, 78-84 (Aug. 1967). 
The interchangeability problem posed by the advent of modular 
instrumentation has been alleviated by the development of a Stand- 
ard Nuclear Module (NIM) system by the AEC Committee on 
Nuclear Instrument Modules. The Committee held its first meeting 
in March 1964. The completed specifications were issued in July 
1964 and laboratory utilization and industry exhibit of the NIM 
system began in November 1964. Acceptance by the laboratories 
and by the nuclear instrument industry has been so rapid and so 
overwhelming that by September 1966 more than 70% of the total 
modular nuclear instruments produced in the U.S. were in the 
NIM system. The percentage today is even higher and the impact 
is world wide. Since all necessary components are available, other 
fields and other industries can readily adopt the system. 


Key Words: Instrument, interchangeability, modules, nuclear, 
standard. 


C. S. Peirce on benefit-cost analysis of scientific activity, 
W. E. Cushen, Operations Res. 15, No. 4, p. 641 (July—Aug. 1967). 
A very early paper on benefit-cost analysis of research projects was 
published in the Appendix to the Coast Survey’s Annual Report 
1876. A special issue of Operations Research, devoted to applica- 
tions, is reprinting the original article, authored by Charles Sanders 
Peirce, who is better known as the father of pragmatism in Ameri- 
can philosophical circles. The paper by Cushen serves to introduce 
the reprint and to paint a brief scene for its use and interpretation. 


Key Words: Benefit-cost, research priorities, Peirce, philosophy, 
pragmatism, history, weights and measures. 


Editor and translator, H. J. Oser, Book, Functional Analysis 
and Numerical Mathematics, by L. Collatz, Academic Press Inc., 
New York, N.Y., 1966. 

Written primarily for graduate students in mathematics and physics, 
this book highlights the dramatic changes which have occurred in 
numerical mathematics during the past twenty years. The increased 
use of electronic computers and the tendency to develop more 
abstract methods are the primary factors in these rapid advances. 
The text is mainly concerned with those parts of functional analysis 
which have proven useful for numerical applications including 
Hilbert space theory, Banach spaces, metric spaces and pseudo- 
metric spaces, the theory of Frechet derivatives, and topological 
fixed-point theorems. Problems concerning iterative methods, 
differential and integral equations, and approximation theory are 
used to illustrate the applications. 


Key Words: Banach spaces, Frechet derivative theory, functional 
analysis, Hilbert space theory, metric spaces, topological fixed- 
point theorems. 


Coping with multiple reflections in precise measurement 
of the speed of light with a microwave Michelson interferom- 
eter, D. M. Kerns (Proc. URSI Symp. Electromagnetic Wave Theory, 
Delft, The Netherlands, September 1965), Book, Electromagnetic 
Wave Theory, pp. 245-249 (Pergamon Press Inc., New York, N.Y., 
1967). 

A microwave version of the Michelson interferometer is currently 
being used for a precise determination of the speed of light, c. In 
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preliminary experiments, multiple reflections between the antenna 
and the movable mirror of the system have been found to produce 
apparent variations in c, which, if uncorrected, would be an order. 
of magnitude greater than tolerable. This paper establishes a theorem 
that should enable substantial elimination of errors due to multiple 
reflections. It is also shown that the four-arm junction in the arrange- 
ment —usually visualized as an ideal magic T—need have no special 
properties and in particular need have none of the ideal properties 
that distinguish an ideal magic 7. 


Key Words: Microwave Michelson interferometer, multiple reflec- 
tions, speed of light. 


Low-frequency electrical calibrations at the National Bureau 
of Standards, F. L. Hermach (Proc. 22nd Annual ISA Conf. and 
Exhibit, Chicago, Ill., Sept. 11-14, 1967), ISA Preprint M5-1- 
MESTIND-67 (May 1967). 

Charts are presented to show the present range and accuracy of 
NBS calibrations of standards of resistance, capacitance, inductance, 
voltage, current, from direct current through 50 kHz. The chains 
of measurements by which these and other calibrations are related 
to the basic NBS standards of voltage and resistance are also shown. 


Key Words: Electrical calibrations, electrical standards, low- 
frequency calibrations, low-frequency standards. 


A note on the evaluation of interference photographs, J. R. 
Meyer-Arendt, Acta Histochem. 22, 58-61 (1965). 

Interference patterns should be evaluated with respect to two 
parameters, namely the spatial displacement and the contrast of 
the fringes. The contrast permits the deduction of the degree of 
wavefront deformation in one of the two bundles brought to inter- 
ference. A series of representative figures are shown in a plot. 


Key Words: Interference photographs, photographs, wavefront 
deformation. 


A simple, inexpensive wavefront shearing interferometer, 
J. B. Saunders, Appl. Opt. 6, No. 9, 1581-1583 (Sept. 1967). 

An efficient wavefront shearing interferometer is easy to produce 
from two small right-angle prisms. A cube beam divider becomes 
such an interferometer when two if its faces are aluminized. Any 
chosen angle of shear is obtained by a simple rotation of one prism 
relative to the other. Chromatic compensation is achieved by sliding 
one prism relative to the other along their common face. 


Key Words: Interferometer, wavefront shearing interferometer, 
prism interferometer. 


Measurement standards, A. G. McNish, IEEE Student J. 4, 
No. 2, 24-31 (Mar. 1966). 

All units of measurement, including the ampere of current and 
candela of light intensity—two base units included in the new 
International System of Units, are ultimately derived from four 
original prototype units—the meter of length, kilogram of mass, 
degree Kelvin of temperature, and second of time. To make all 
measurement units more effective tools, the standards that physically 
embody or define three of the prototypes—length, time, and tem- 
perature—have been changed in recent years. 


Key Words: Length, mass, temperature, time, units of measure- 
ments. 


The functions of Commission I of the International Scientific 
Radio Union, J. M. Richardson, Proc. IEEE 55, No. 6, 743-745 
(June 1967). 

Within the context of the objectives and characteristics of URSI, 
the functions of its Commission I on radio standards and measure- 
ments is given. It deals with the subject in the framework of careful 
survey, evaluation, and guidance of precision of standards and 
measurements; the attention is to the highest national levels of 
accuracy for the long term; the involvement is by people technically 
experienced and also in positions of enough responsibility to carry 
out the conclusions of the commission. The nature of its international 
and national activities, especially in the U.S., is given. 


Key Words: Radio, standards, measurements, International Scien- 
tific Radio Union (URSI). 


Ionosphere, T. E. VanZandt, Encyclopedia of Physics, Ed. R. M. 
Besancon, pp. 348-349 (Reinhold Publ. Co., New York, N.Y., 1966). 
The structure and physics of the ionosphere are briefly reviewed. 
The ionospheric regions are named. Then for each region are given: 
its height range, ionic composition, electron concentration, causes 
of ionization, and loss processes. Transport of ionization is men- 
tioned. Sporadic E and spread F are described. 


Key Words: Electron concentration, height range, ionic composition, 
ionosphere, ionospheric regions. 


Mechanistic aspects of hearing, E. L. R. Corliss, J. Acoust. Soc. 
Am. 41, No. 6, 1500-1516 (June 1967). 

The results of a number of apparently diverse experiments on human 
hearing can be summarized via a simple model. Two mechanisms 
are required: One element emits unit responses whose number is 
proportional to the instantaneous amplitude of a sound. The other 
is a selector mechanism, analogous to a resonant circuit, that 
integrates the responses. Correspondingly, just two parameters are 
involved: the threshold of hearing and the “Q”’ of the selector 
mechanism. The system of equations representing the model makes 
it possible to derive the necessary parameters from the results of 
a number of different types of experiments. The theory of the 
selector mechanism is based upon a recent study of the limits of 
performance of a system capable of storing oscillatory energy 
reversibly, but the remainder of the theoretical description makes 
use of ordinary communication theory. 

The discrimination ability of the normal ear is shown to be propor- 
tional to the one-fourth power of the level above threshold, and to 
relate closely to the sensation of loudness. Many features of the 
behavior of normal and impaired ears can be predicted by applica- 
tion of familiar network theory to the parameters derived from the 
model. The properties of the model may be used as criteria in the 
choice of experiments to develop explicit data on the way the ear 
performs its function. 


Key Words: Hearing, mechanistic model of hearing, communication 
theory, resolving power of the ear, ear behavior, auditory perception, 
auditory analysis. 


Second breakdown—A comprehensive review, H. A. Schafft, 
Proc. IEEE 55, No. 8, 1272-1288 (Aug. 1967). 

This paper is a comprehensive review of the published literature 
dealing with the phenomenon of second breakdown in semiconductor 
devices and the problems it creates in the design, fabrication, testing, 
and application of transistors. 


Key Words: Review, semiconductor devices, transistors, diodes, 
second breakdown, thermal breakdown, transistor reliability, circuit 
reliability, failure modes, failure mechanisms. 


Superconductive properties of ceramic mixed titanates, 
J. F. Schooley, H. P. R. Frederikse, W. R. Hosler, and E. R. Pfeiffer, 
Phys. Rev. 159, No. 2, 301-305 (July 10, 1967). 

The superconducting transition temperatures of a series of mixed 
barium-strontium and calcium-strontium titanate ceramic specimens 
have been determined. For several of these samples the low-field 
magnetization has also been investigated. The results are compared 
with those obtained from single-crystal SrTiO3. 


Key Words: Barium, ceramic mixed titanates, low field magnetiza- 
tion, strontium, superconducting transition temperatures, super- 
conductive properties. 


The radiative band pinch effect and temperature depend- 
ence of radiative recombination in GaAs, N. N. Winogradoff, 
K. Owen, and R. M. Curnutt, Int. J. Electronics 22, No. 3, 229-233 
(Mar. 29, 1967). 

Localized heating of GaAs by intense ruby laser flashes cause a 
reduction or “pinching” of the band gap in the illuminated region. 
The increase in wavelength of the recombination radiation from this 
region permits the study of the temperature dependence of the 
internal quantum efficiency without absorption in the colder, wide 
band gap regions. 
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The results show that the internal quantum efficiency decrease 
with an increase in temperature. 


Key Words: GaAs, radiative recombination, temperature depend- 
ence, quantum efficiency, ruby laser. 


Another method of synthesizing nonuniformly spaced an- 
tenna arrays, M. T. Ma, Trans. IEEE Ant. Prop. AP-13, No. 5, 
833-834 (Sept. 1965). 

A formulation for synthesizing an optimum nonuniformly spaced but 
symmetrical array by adjusting both the amplitude excitations and 
element spacings is given. This is accomplished by applying Haar’s 
theorem which is known in a branch of mathematics. It is shown 
that the solution obtained according to the method proposed here 
is optimum in the senses that, with respect to a chosen set of ele- 
ment spacings, 1) the maximum deviation between the synthesized 
and desired patterns is minimized, 2) the side lobes can be made 
equal and their level minimum for a specified beamwidth, 3) the 
side lobe level and beamwidth are not mutually constrained so that 
a solution better than the Dolph-Chebyshev array is possible, 4) a 
minimum number of elements required to synthesize a desired 
pattern can be determined, and 5) the solution is unique. 


Key Words: Amplitude excitations, antenna arrays, element spac- 
ings, nonuniformly spaced antenna arrays, synthesizing nonuniformly 
spaced antenna arrays. 


Antennas, H. V. Cottony, Encyclapedia of Physics, Ed. R. M. 
Besancon, pp. 25-28 (Reinhold Publ. Co., New York, N.Y., 1966). 
The functions and properties of antennas are outlined in a 1500 word 
article for use in an encyclopedia. The presentation is presented 
from a view point that an antenna is a network designed to facilitate 
the transfer of information energy from circuits to space and reverse. 
The relationship between an antenna and the contiguous electronic 
circuit is determined and fully described by the terminal impedances 
of the circuit and antenna. Because of the three-dimensional nature 
of propagating medium and the nature of electromagnetic wave, 
the performance of an antenna as a radiator is more complex. For 
full description, it requires the specifications of its radiation pattern 
in space, and the polarization of the radiated electromagnetic wave. 
Six illustrations of different ‘types of antennas are presented. 


Key Words: Antennas, circuits, contiguous electronic circuit, 
electronic circuit, information energy, terminal impedance. 


The development of an International atomic time sacle, 
J. A. Barnes, Proc. IEEE 55, No. 6, 822-826 (June 1967). 

The paper reviews briefly the methods of generating atomic time 
and the errors inherent in the resulting scales. An atomic clock 
consists of an atomic frequency standard and an “integrator” to 
accumulate the phase of the signal. Because of noise perturbing 
the instantaneous frequency, an ensemble of identical atomic 
clocks will show a distribution of (epoch) times which is unbounded 
as the system evolves in time. The recognition of this problem has 
important consequences in national and international coordination 
of time scales and the construction of average atomic time scales. 
Also of significance is the not completely resolved question of 
weighting of individual standards in the construction of average time 
scales. In spite of these difficulties it is pointed out that through 
coordination and proper data handling, most of the advantages of 
astronomical time scales can be realized by atomic time scales. A 
statement of some of the problems facing any attempts at coordina- 
tion is presented without any suggested solution. 


Key Words: Atomic time, coordination, atomic frequency standards. 
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